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PREFACE 


A request  for  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  to  conduct  an  investigation  of  wave  heights  on  the  Great  Lakes  was 
made  by  the  U.  S.  Army  Engineer  Division,  North  Central  (NCD),  in  a con- 
ference held  in  Chicago,  Illinois,  on  22  July  197*+.  Funds  were  autho- 
rized by  NCD  on  30  August  197 *+•  The  study  was  conducted  during  the 
period  from  September  197^  to  June  1976  in  the  Coastal  Branch,  Wave 
Dynamics  Division,  Hydraulics  Laboratory,  WES,  under  the  direction  of 
Mr.  H.  B.  Simmons,  Chief  of  the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin, 
Chief  of  the  Wave  Dynamics  Division,  and  Dr.  C.  L.  Vincent,  Chief  of  the 
Coastal  Branch. 

Drs.  D.  T.  Resio  and  C.  L.  Vincent  conducted  the  study  and  also 
prepared  the  report.  During  the  investigation,  Mrs.  Rebecca  Brooks 
and  Mr.  W.  D.  Corson  were  especially  helpful  in  performing  analytical 
and  programming  tasks.  Dr.  L.  H.  Blakey  and  Messrs.  N.  Arno  and  L. 
Hiipakka  of  NCD  exhibited  a keen  perception  of  the  priority  need  and 
practical  use  of  this  information  and  contributed  valuable  suggestions 
to  the  form  of  the  design  wave  information  required  on  the  Great  Lakes. 

A significant  proportion  of  the  numerical  computations  were  per- 
formed on  a CDC-7600  computer  at  the  U.  S.  Army  facility  in  Huntsville, 
Alabama.  Mr.  Fred  Bourgeois,  U.  S.  Army  Engineer  Division,  Huntsville, 
assisted  significantly  in  coordinating  visits  to  Huntsville  and  sched- 
uling computer  time. 

Directors  of  WES  during  the  conduct  of  the  study  and  the  prepara- 
tion of  this  report  were  COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE. 
Technical  Director  was  Mr.  F.  R.  Brown. 
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CONVERSION  FACTORS,  METRIC  (Si)  TO  U.  S.  CUSTOMARY  AND 
U.  S.  CUSTOMARY  TO  METRIC  (Si)  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in 

Multiply 

this  report 

By 

can  be  converted  as  follows : 

To  Obtain 

Metric 

(SI)  to  U.  S. 

Customary 

cent imetres 

0.3937 

inches 

metres 

3.28081* 

feet 

square  metres 

10.7638 

square  feet 

Celsius  degrees  or  Kelvins 

9/5 

Fahrenheit  degrees* 

U.  S.  Customary  to  Metric  (Si) 

feet 

0.30U8 

metres 

miles  (U.  S.  statute) 

1.60931*1* 

kilometres 

miles  (U.  S.  nautical) 

1.852 

kilometres 

miles  per  hour 

1.60931+1* 

kilometres  per  hour 

knots  (international) 

0.511*1*1*1+1* 

metres  per  second 

square  feet 

0.09290301* 

square  metres 

degrees  (angle) 

0.01 71*5329 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins** 

* To  obtain  Fahrenheit  (F)  temperature  readings  from  Celsius  (C)  read- 
ings, use  the  following  formula:  F = 9/5 (c)  + 32.  To  obtain  Fahren- 
heit readings  from  Kelvins  (k),  use:  F = 9/ 5 ( K - 273.15)  + 32. 

**  To  obtain  Celsius  (c)  temperature  readings  from  Fahrenheit  (F)  read- 
ings, use  the  following  formula:  C = 5/9 (F  - 32).  To  obtain  Kelvin 
(K)  readings  from  Fahrenheit  readings,  use:  K = 5/9(F  + 1*59-67). 
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DESIGN  WAVE  INFORMATION  FOR  THE  GREAT  LAKES 


, LAKE  HURON 

PART  I : INTRODUCTION 

1.  The  primary  purpose  of  this  study  is  to  provide  wave  informa- 
tion for  planned  sites  of  dredged  material  retaining  structures  on  the 
Great  Lakes.  Secondary  goals  are  the  provision  of  wind-field  informa- 
tion over  the  Great  Lakes  and  detailed  wave  information  along  U.  S. 
shorelines  of  the  Great  Lakes.  This  report  deals  with  wave  conditions 
on  Lake  Huron. 

2.  Past  studies  of  waves  on  Lake  Huron  have  generally  been  made 
on  a site-to-site  basis.  Very  few  data  are  available  from  wave  gages 
(Figure  l).  Visual  wave  observations  are  available  from  Coast  Guard 
stations  at  several  shore  locations  (Figure  2)  and  from  commercial  ships 
for  waves  offshore.  Data  from  wave  hindcasts  are  also  available  for 
several  locations  (Figure  3).  For  a specific  site,  inconsistencies 
exist  between  these  data  that  make  them  unreliable  for  estimates  of  re- 
currence intervals  (Figure  4).  The  sensitivity  of  economic  analyses  of 
optimal  design  criteria  to  even  small  variations  in  recurrence  intervals 
of  large  waves  clearly  demonstrates  an  acute  need  for  a reliable  source 
of  wave  data  as  well  as  an  evaluation  of  the  distribution  of  errors  in 
these  data. 

3.  The  sparseness  of  reliable  wave  information  is  critical  for 
many  areas  on  Lake  Huron.  Visual  observations  from  shore  cannot  be 
extrapolated  to  other  sites  along  a shore  due  to  the  unknown  effects  of 
shallow-water  transformations  and  shorelines  geometry.  Visual  obser- 
vations from  ships  provide  estimates  of  waves  offshore  but  are  dif- 
ficult to  treat  properly  in  a shallow-water  environment.  Although  a 
network  of  wave  gages  might  eventually  provide  good  data,  the  expense 
involved  makes  it  economically  infeasible  to  provide  detailed  coverage 
of  the  entire  Lake  Huron  shoreline.  Even  if  such  a wave  gage  network 
were  established,  there  would  be  a lag  time  of  many  years  before 
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Figure  1.  Location  of  wave  gage  site  on  Lake  Huron 
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MICHIGAN 


COAST  GUARD  STATIONS  PERIOD  OF  RECORD 

1 ST.  IGNACE  1970  - 1971 

2 MACKINAW  CITY  1914-1968 

3 THUNDER  BAY  1945  - 1971 

4 TAWAS  POINT  1944  - 1971 

5 HARBOR  BEACH  I960  - 1971 

6 PORT  HURON  1941  - 1971 

Figure  2.  Locations  of  visual  wave  observation  sites  on  Lake  Huron 


1  COCKBURN  ISLAND,  ONTARIO 


1965  - 1967 
1965  - 1967 
1965  - 1967 
1965  - 1967 
196  5 - 1967 


COLE  AND  HILFIKER 


2 ALPENA.  MICHIGAN 

3 HARBOR  BEACH,  MICHIGAN 

4 PORT  HURON,  MICHIGAN 

5 DOUGLAS  POINT,  ONTARIO 

Figure  3.  Locations  of  sites  on  Lake  Huron  hindcast 


SOURCE  OF  WAVE  DATA 


Ial-2  HINDCASTS  BY  GLENN 
lb)  HINDCASTS  BY  SAVILLE  3 
Id  HINDCASTS  BY  LEE  a 
Hal  SHIPS  OBSERVATIONS 
Ibl-2  SHIPS  OBSERVATIONS 
lcl*2  COAST  GUARD  (VISUAL-) 
Hal  WAVE  GAGE 


WAVE  HEIGHT,  FT 

Figure  4.  Chart  showing  range  of  estimates  of  design  wave  heights 
for  Cleveland.  Arabic  numbers  1 and  2 refer  to  full  population 
and  asymptotic  methods  of  extreme  analysis,  respectively 


sufficient  wave  information  was  available  for  design  purposes.  Hind- 
casted  waves  can  provide  data  along  the  entire  Lake  Superior  shoreline; 
however,  there  are  many  problems  that  must  be  overcome,  if  reliable 
hindcasts  sire  to  be  obtained.  The  treatment  of  some  of  these  and  the 
selection  of  an  appropriate  wave  model  will  be  discussed  in  detail  in 
Part  III,  along  with  comparisons  of  this  numerical  wave  model  to  the 
Sverdrup-Munk-Bretschneider  (SMB)  technique  described  in  the  Shore  Pro- 
tection Manual  (SPM).^ 

H.  In  spite  of  uncertainties  of  hindcast  methods,  the  immediate 
need  for  wave  data  on  Lake  Huron  makes  an  approach  via  hindcasting  the 
only  viable  alternative.  Historical  wind  fields  can  be  used  to  estab- 
lish a length  of  record  up  to  22  yr  for  Lake  Huron.  The  sparsity  of 
stations  and  large  lake  surface  area  make  any  effort  to  hindcast  beyond 
22  yr  dubious. 

5.  The  largest  wind-generated  waves  on  the  Great  Lakes  are  those 
produced  by  synoptic-scale  circulation  features  such  as  extratropical 
and  tropical  cyclones.  Smaller-scale  wind  phenomena,  such  as  individual 
thunderstorms  or  squall  lines,  do  not  maintain  high  winds  over  suffi- 
cient fetch  to  generate  waves  of  comparable  height.*  Consequently,  the 
study  of  extreme  wave  conditions  on  the  Great  Lakes  is  equivalent  to  the 
study  of  extreme  synoptic  events. 

6.  The  use  of  a hindcast  technique  to  estimate  wave  heights 
assumes  that  the  coupling  between  the  atmospheric  boundary  layer  and 
waves  generated  by  the  motion  in  this  boundary  layer  is  known.  Other 
factors,  such  as  the  interaction  between  waves  and  the  lake  bottom  and 
interactions  between  different  spectral  wave  components,  also  must  be 
considered.  If  all  these  effects  can  be  treated  adequately  and  if  the 

tmospheric  motion  over  a lake  is  known,  then  a reliable  estimate  of 
waves  can  be  made.  In  any  estimate  of  wave  height,  there  will  always 
be  some  error  whether  the  estimate  is  from  a wave  gage,  visual 

* This  refers  only  to  waves  with  periods  of  less  than  30  sec.  Long 
waves  can  be  generated  by  a resonant  interaction  with  a front  or 
squall  line,  but  this  is  essentially  a different  process  than  the 
growth  of  wave  spectra  covered  in  this  report. 


observations,  or  hindcast  techniques.  The  essential  question  that  must 
be  answered  for  engineering  and  planning  purposes  is  how  much  error  is 
there  in  these  methods  relative  to  the  requisite  accuracies  for  a given 
application.  In  order  to  facilitate  answering  this  question  for  poten- 
tial users  of  wave  information,  an  effort  will  be  made  to  include  con- 
trol bands  for  all  results  presented. 

7 • Figure  5 demonstrates  an  overview  of  the  approach  to  be  fol- 
lowed in  the  estimation  of  wave  information  for  the  design  of  dredged 
material  retaining  structures.  The  general  approach  is  valid  as  well 
for  any  application  in  which  information  is  needed,  particularly  in  the 
range  of  extreme  values. 

OVERVIEW  OF  WAVE  INFORMATION  PROJECT 
FOR  GREAT  LAKES 


WINDS 

o . LAND  STATIONS 

b.  SHIPS  OBSERVATIONS 

c.  WEATHERMAPS 

I 

STORM  WINDS 


t 

WAVE  MODEL 
(NUMERICAL) 


STORM  WAVES 

I 

WAVE  CLIMATE  PARAMETERS 

a.  RECURRENCE  INTERVALS 

b.  STORM  HYDROGRAPHS 

Figure  5.  Overview  of  hindcasting 
approach  to  design  wave  estimates 
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PART  II:  ESTIMATION  OF  WINDS 

8.  Lack  of  continuous  wind  observations  over  the  Great  Lakes 
forces  the  conclusion  that,  for  practical  purposes,  winds  over  a lake 
must  be  estimated  from  other  sources  of  data.  Time  and  space  scales  of 
wave  generation,  as  well  as  the  adaptability  to  computer  processing, 
must  be  considered  in  choosing  the  methodology  for  these  wind  esti- 
mates. The  only  three  sources  of  data  with  sufficient  length  of  record 
to  permit  calculation  of  extremes  for  the  Great  Lakes  are  (a)  pressure 
observations  at  adjacent  land  stations,  (b)  synoptic  weather  maps,  and 
(c)  wind  observations  at  adjacent  land  stations.  Recent  studies  by 
Cole  and  Hilfiker^  have  shown  that  a network  of  pressure  data  around  a 
lake  can  produce  reasonable  forecasts  of  winds  over  a lake.  However,  a 
study  by  Resio  and  Vincent0  demonstrated  that  the  use  of  suitably  trans- 
formed overland  winds  from  around  a lake  can  also  provide  an  accurate 
overlake  wind  field.  The  ease  of  data  manipulation  and  wind  field  inter- 
polation (Appendix  C)  in  this  technique  make  it  advantageous  to  consider 
the  use  of  these  wind  transformations  for  lake  hindcasts.  Figure  6 
shows  the  locations  of  the  major  sites  at  which  wind  data  have  been 
recorded  on  Lake  Huron. 

9.  Reports  1 (Lake  Erie)  and  2 (Lake  Ontario)  of  this  series  and 
WES  Miscellaneous  Paper  H-76-12^  established  the  relationship  between 
winds  over  land  and  over  the  lakes.  Theoretical  and  empirical  support 

of  the  technique  employed  is  detailed  in  this  Miscellaneous  Paper.  Thus, 

specific  examples  presented  in  this  report  may  be  taken  from  Lakes  Erie 

or  Ontario  with  the  data  from  Lake  Huron  serving  as  a supplement  to 

increase  the  general  confidence  in  earlier  results.  The  concept  is  to 
* 

accumulate  general  data  in  support  of  the  proposed  wind  and  wave  tech- 
niques wiLh  each  report  in  this  series  rather  than  to  treat  each  lake 
independently . 

7 

10.  Cole  reviewed  a number  of  techniques  used  for  estimating 
winds  over  the  Great  Lakes.  His  conclusions  indicated  that  none  of  the 
presently  used  techniques  correlated  well  with  winds  observed  over  a 
lake.  Table  1 gives  the  results  Cole  obtained  in  his  study.  Figure  7 
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1 

SAULT  STE.  MARIE,  MICHIGAN 

1948  - 1973 

2 

ALPENA,  MICHIGAN 

1959  - 197  4 

3 

OSCODA-WURTSMITH.  MICHIGAN 

1950  - 1970 

4 

SAGINAW,  MICHIGAN 

1948  - 1954 

5 

SEL FRIDGE,  MICHIGAN 

1949  - 1970 

6 

CENTRALI A,  ONTARIO 

1953  - 1966 

7 

WIARTON,  ONTARIO 

1953  - 1972 

8 

GORE  BAY,  ONTARIO 

1953  - 1972 

Figure  6.  Locations  of  major  sites  at  which  wind  data  have  been 

recorded  for  Lake  Huron 
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Figure  7*  Schematic  diagram  of  wind  techniques 
investigated  by  Cole? 

shows  the  analyses  used  in  each  of  the  methods  listed  in  Table  1. 

8 9 

The  Bretschneider  and  Jacobs  techniques  would  have  to  be  taken  from 
weather  maps  and  involve  exceedingly  large  amounts  of  hand  processing. 
Only  the  comparison  to  the  Richards  et  al.10  technique  is  relevant  to 
the  method  to  be  used  in  the  present  study.  This  comparison  indicates 
a very  poor  correlation  of  0.36  between  the  l6-m*  Richards  wind  and  winds 
observed  at  a tower  in  Lake  Michigan.  On  the  basis  of  such  an  inade- 
quate wind  representation,  the  use  of  land-wind  information  to  estimate 
over lake  winds  might  be  dismissed.  However,  a possible  explanation  of 
the  lack  of  correlation  in  this  study  will  be  offered  after  a more  de- 
tailed discussion  of  the  general  problem. 

11.  The  interaction  between  a marine  surface  and  the  overlying 
atmosphere  is  very  complex.  Only  at  some  large  distance  from  the  inter- 
face (somewhere  in  the  range  of  a kilometre)  can  the  effect  of  the  bound- 
ary be  neglected.  At  this  level,  winds  at  the  latitude  of  the  Great 
Lakes  can  be  treated  as  geostrophic  or,  if  isobar ic  curvature  is 


* A table  of  factors  for  converting  metric  (Si)  units  of  measurement  to 

U.  S.  customary  units  and  U.  S.  customary  units  to  metric  (Si)  units 
is  presented  on  page  3. 
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significant,  gradient.  As  previously  mentioned,  synoptic-scale  systems, 
such  as  extratropical  cyclones,  produce  the  dominant  large-wave  con- 
ditions on  the  Great  Lakes.  The  size  of  a typical  storm  of  this  type 
is  shown  in  Figure  8.  From  the  patterns  of  isobars  in  this  figure,  it 


NOVEMBER  19,  1879  3 PM 

Figure  8.  Isobaric  pattern  of  a typical  extratropical 
low  in  the  Great  Lakes  region 

is  evident  that  the  geostrophic  winds  above  the  Great  Lakes  are  not 
different  from  the  winds  over  adjacent  land  stations.  Hence,  the 
driving  mechanism  for  lower  level  winds  in  each  case  is  the  same. 

The  fundamental  differences  between  the  lower  level  winds  must  then 
be  primarily  a function  of  the  boundary  characteristics . This  means 
that,  if  the  transformations  from  geostrophic  winds  to  overland  winds 
and  geostrophic  winds  to  over lake  winds  were  both  known,  a relation- 
ship between  land  winds  and  lake  winds  could  be  defined.  In  practice, 
this  is  not  an  easy  task.  Most  attempts  to  relate  land  winds  to  lake 
winds  have  been  purely  empirical.  For  comparison  of  lake  winds  to 

ik 


where  U and  represent  wind  speed  over  water  and  over  land,  re- 

spectively. Then  Hunt,1'*  Lemire,1^  Richards Richards  et  al.  and 

l8 

Richards  and  Phillips  used  this  ratio  to  study  the  relationship  be- 
tween land  winds  and  lake  winds.  Table  2,  taken  from  Reference  18,  gives 
a summary  of  results  obtained  in  these  studies.  The  study  by  Richards 
et  al.  also  documented  the  effects  of  over-the-water  fetch,  air-water 
temperature  differences,  and  wind  velocity  on  this  ratio  (Figure  9). 

From  this  study,  it  is  apparent  that  the  effect  of  air-water  temperature 
differences  on  R is  decreased  with  increasing  wind  velocity,  which,  for 
the  investigation  of  extreme  wave  conditions,  can  be  very  significant. 

12.  The  systematic  behavior  of  the  wind  ratio  suggests  that  these 
results  are  not  spurious  but  actually  represent  a real  relationship  be- 
tween land  winds  and  lake  winds.  However,  the  studies  do  not  indicate 
to  what  degree  the  relationship  was  reproducible  in  a given  application, 
i.e.,  there  is  no  correlation  coefficient  calculated  between  observed 
lake  winds  and  winds  estimated  using  the  relationship.  An  evaluation 
of  the  error  in  the  wind  estimate  is  essential  to  the  determination  of 
errors  in  wave  hindcasts.  Thus,  additional  empirical  work  was  necessary 
to  complete  the  evaluation  of  the  relationship  between  lake  winds  and 
land  winds. 

l8 

13.  Richards  and  Phillips  points  out  that  the  wind  ratio  R 
"combines  the  influences  which  cause  winds  over  the  water  to  vary  from 
winds  on  the  land;  i.e., 

(i)  differences  in  frictional  effects  created  by  land  and 
water ; 

(ii)  changes  in  atmospheric  stability  created  by  air-water 
temperature  differences;  and 

* For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de- 
fined in  the  Notation  (Appendix  F). 
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(iii)  the  length  of  overwater  fetch." 

When  errors  in  estimated  winds  are  examined  relative  to  observed  winds 
over  a lake,  another  factor,  inhomogeneities  in  the  wind  field,  must  be 
considered.  Many  scales  of  motion  are  superimposed  in  the  planetary 
boundary  layer.  Any  observation  of  wind  speed  must  be  regarded  as  a 
statistic  containing  elements  from  each  of  these  scales  of  motion.  Thus, 
a comparison  between  two  winds  at  different  locations  must  be  carefully 
scrutinized  to  retain  only  the  scales  of  motion  that  are  being  investi- 
gated. In  the  present  study,  synoptic-scale  features  are  being  exam- 
ined. If  a single  land  station  is  correlated  with  the  ships  observa- 
tions, there  is  no  consideration  given  to  spatial  gradients  in  wind 
velocity.  Results  from  such  a correlation  analysis  can  be  quite  mis- 
leading since  they  contain  sources  of  variability  other  than  those 
attributable  to  inherent  differences  between  land  winds  and  lake  winds. 
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To  emphasize  this  point,  a correlation  coefficient  was  computed  between 
simultaneous  wind  observations  at  Toledo  and  Cleveland  airports  (using  a 
random  sample  of  50  observations  during  December  1972).  The  calculated 
correlation  was  0.85,  which  means  that  only  72  percent  of  the  variance 
at  Cleveland  could  be  explained  using  wind  data  from  Toledo.  On  close 
examination,  much  of  this  difference  appears  to  be  related  to  shifts  in 
the  time  of  peak  winds  during  the  passage  of  storms.  For  example,  a 
storm  moving  at  30  mph  from  west  to  east  along  a line  from  Toledo  to 
Cleveland  might  take  about  3 hr  to  advect  the  same  wind  conditions  from 
the  Toledo  area  to  the  Cleveland  area.  The  remainder  of  the  difference 
is  most  likely  due  to  mesoscale  and  smaller  variations  in  the  wind  field. 

lU.  It  should  be  noted  that  the  land  winds  for  airport  stations 
typically  represent  the  average  of  two  1-min  wind  observations  per  hour. 
Cole^  does  not  discuss  the  averaging  length  for  the  winds  observed  over 
the  lake,  which  together  with  the  neglect  of  spatial  gradients  in  the 
wind  field,  may  explain  why  the  correlation  between  lake  winds  and  land 
winds  was  so  low  in  his  study.  The  ships  observations  used  in  this 
study  are  approximate  winds  averaged  visually  by  an  observer  watching 
an  anemometer  for  a 1-min  period.  Deviations  from  such  an  observation 
schedule  might  be  expected  during  severe  weather,  but  for  most  obser- 
vations the  quality  control  appears  to  be  quite  high.  Figure  10  shows 
a comparison  between  several  sequences  of  winds  observed  at  Toledo  and 
Cleveland  and  simultaneous  observations  by  ships  on  the  lake  between 
longitude  8l.5  deg  and  the  western  end  of  Lake  Erie.  As  seen  in  this 
figure,  there  is  considerable  variance  between  the  ships  observations 
for  the  same  time;  yet,  if  an  average  wind  is  calculated  from  the  ships 
observations,  a consistency  is  apparent  between  the  land  winds  and  lake 
winds.  To  test  this  result  with  the  entire  data  sample,  regressions 
were  run  between  the  ships  observations  and  the  mean  wind  calculated 
from  Toledo  and  Cleveland.  These  regressions  are  stratified  by  the  num- 
ber of  ships  observations  per  3-hr  increment  of  time  (Table  3).  Clearly, 
there  is  an  increase  in  the  correlation  as  one  moves  from  one  observa- 
tion per  time  increment.  Hence,  the  actual  correlations  between  a 1-hr 
average  land  wind  and  a 1-hr  average  lake  wind  should  tend  to  be  well 
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Figure  10.  Wind  sequences  from  Toledo  and  Cleveland  (Lake  Erie) 
compared  to  winds  observed  over  the  lake  ships 


correlated  with  correlation  coefficients  of  0.8  or  better. 

15.  Extensive  wind  data  are  available  for  eight  stations  sur- 
rounding Lake  Huron:  Sault  Ste.  Marie,  Self ridge  AFB,  Alpena,  Oscada 
Wurtsmith  AFB,  Saginaw,  Wiarton,  Gore  Bay,  and  Centralia  (Figure  6). 
These  records  are  stored  on  magnetic  tapes  and  are  thus  easily  retriev- 
able for  the  complete  22  yr  of  hindcasts,  but  a minimum  of  four  stations 
was  maintained  for  all  hindcast  storms.  In  this  report,  the  basic  form 
for  ratio  R derived  in  Reference  6 is  expressed  as 

R = tp  4>  (2) 

v n 

where  i|»  is  a velocity-dependent  function  (Figure  A2  in  Appendix  A), 
and  is  a normalized  function  dependent  on  air-sea  temperature 

difference . 

16.  Table  U presents  the  root-mean-square  (rms)  errors  of  wind 
speeds  estimated  from  land  stations,  using  Equation  2,  compared  to  wind 
speeds  observed  by  ships.  As  seen  in  this  table,  the  expected  rms  error 
for  most  conditions  of  interest  to  hindcasts  is  under  5 knots. 

17.  Reports  1 and  2 indicated  a linear  relation  between  the  ratio 

of  anemometer-level  winds  to  geostrophic  winds  and  air-sea  temperature 
difference.  However,  analyses  documented  in  Reference  6 have  shown  that 
although  the  portion  of  curve  studied  in  the  earlier  reports  is  linear, 
the  complete  curve  is  of  the  form  shown  in  Figure  11.  The  relationship 
is  given  in  terms  of  the  normalized  function  in  this  figure.  In 

order  to  convert  V/G  , the  ratio  of  anemometer-level  wind  speed  to 

Q 

geostrophic  wind  speeds,  to  <J>n  , the  curve  from  Bretschneider  was 
normalized  by  dividing  by  the  value  of  V/G  for  neutral  stability. 

18.  Using  a two-layer  model  of  atmospheric  motion  in  the  plane- 

19 

tary  boundary  layer,  Cardone  showed  that  the  linear  portion  of  the 

variation  of  <t>  with  air-sea  temperature  difference  was  theoretically 

n 20 
justifiable.  Figure  12  presents  empirical  data  from  Bijvoet,  H.O. 

12  21 

Pub.  60U,  and  Johnson  that  also  support  the  slope  of  the  linear 
portion  of  the  curve  shown  in  Figure  11. 

19.  The  dependence  of  the  ratio  R on  velocity  is  clearly 
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AIR-SEA  TEMPERATURE  DIFFERENCE,  DEGREES  CENTIGRADE 

Figure  11.  Comparison  of  empirical  results  of  Richards  et  al.,10 
Bret schneider ,8  and  the  present  study.  The  data  of  Richards 
et  al.  and  Bretschneider  were  normalized  to  1 at  neutral 


and  air-sea  temperature  differences  as  determined  in  the  present 
study  and  published  empirical  relationships  of  the  ratio  of 
anemometer -level  wind  to  geostrophic  wind  speed.  A value  of  V/G 
= 0.1+5  at  Ta  - Ts  = 0 was  assumed  in  constructing  the  line 
for  the  present  study  (from  the  "b"  coefficients) 


established  in  several  studies  (Johnson,21  Richards  et  al.,1^  and  Resio 
and  Vincent^).  Figure  13  shows  a comparison  of  the  results  from  these 
three  studies.  When  the  statistical  nature  of  the  different  measure- 
ments and  analyses  used  in  these  independent  studies  is  considered,  the 
consistency  is  striking. 

20.  The  conclusion  of  this  investigation  is  that  the  use  of  winds 
around  the  lake  perimeter  appears  to  give  an  adequate  representation  of 
winds  over  a lake.  Certainly,  the  representation  is  quite  limited  in 
its  ability  to  reproduce  effects  of  mesoscale  and  smaller  variations  in 
the  wind  field.  What  is  obtained  in  this  analysis  is  a synthesized 
representation  of  the  synoptic-scale  winds  over  a lake.  Only  calibra- 
tion of  the  wave  models  using  lake  winds  estimated  from  land  winds  can 
answer  the  question  as  to  the  adequacy  of  this  representation  for  hind- 
casting purposes. 
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NORMALIZED 


vELOCi  i Y,  KNOTS 

Figure  13.  Velocity-dependent  behavior  of  the  ratio  of 
anemometer-level  wind  speeds  (V)  to  geostrophic  wind 
speeds  (G).  Values  for  V/G  have  been  normalized  to 
neutral  conditions  assuming  a value  of  V/G  for  land 

of  0.^5 
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PART  III:  APPLICATION  OF  WAVE  MODELS 


21.  Theoretical  treatments  of  wave  growth  due  to  the  action  of 

22 

the  wind  date  back  to  Helmholtz.  Since  then,  a proliferation  of 

theories  has  been  offered  to  explain  this  phenomenon.  It  was  not  until 

the  19^0's,  however,  that  practical  wave  hindcasting  and  forecasting 

11  12 

began.  Sverdrup  and  Munk  ’ presented  semi-empirical  methods  for 
calculating  wave  growth  due  to  winds.  These  methods  used  a single  mea- 
sure of  wave  height,  "significant  wave  height,"  and  a single  measure  of 
wave  period,  "significant  period,"  to  characterize  waves  on  the  ocean 
surface.  The  significant  wave  height  was  defined  as  the  average  of  the 
one-third  highest  waves  in  an  observation  period  and  was  intended  to 
correspond  to  that  wave  height  estimated  visually  by  an  observer.  The 
significant  period  was  meant  to  correspond  to  the  period  of  the  "sig- 
nificant"  waves.  Johnson  applied  the  Buckingham  Pi  theorem  to  obtain 
nondimens ional  relationships  between  wave  heights,  wave  periods,  wind 
speed,  fetch,  and  duration  of  winds.  The  form  of  these  relationships 
has  remained  essentially  unchanged.  The  only  real  modifications  have 

been  based  on  changing  the  values  of  some  of  the  empirical  coefficients 

8 2b 

when  warranted  by  additional  data  (Bretschneider  ’ ). 

22.  A different  type  of  representation  of  the  ocean  surface  also 

began  during  the  late  19^0's  when  investigators  of  random  noise  in 

circuits  (Rice,  Tukey  and  Hamming  ) began  formulating  spectral 

27 

relationships  for  stationary  time  series.  Pierson  took  these  concepts 

and  applied  them  to  ocean  waves.  In  this  characterization,  waves  are 

visualized  as  being  a sum  of  an  infinite  number  of  wave  trains  composed 

28 

of  varying  frequencies  and  directions  of  travel.  Longuet-Higgins^ 
showed  that,  under  certain  constraints,  the  distribution  of  wave  heights 
in  a spectrum  was  specified  by  a Rayleigh  distribution  and  formally 
equated  the  concept  of  "significant  wave  height"  with  the  energy  spec- 
trum by  the  relationship 


“i/3  - 2'83  * 


(3) 
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where  is  the  significant  wave  height  and  E is  twice  the  total 

variance  in  the  wave  spectrum.  Several  theories  have  since  been  ad- 
vanced on  the  growth  of  wave  spectra  under  the  action  of  the  wind 
(Appendix  B).  These  provide  the  framework  for  several  different  spec- 
tral hindcast  models  that  exist  today. 

23.  A definite  choice  had  to  be  made  in  this  study  whether  to  use 
a significant  wave  method  (sometimes  referred  to  as  an  SMB  method)  or  a 
spectral  model,  and  if  a spectral  model  were  to  be  used,  which  one  would 
be  selected. 

2k.  To  choose  an  optimal  model  for  the  Great  Lakes  requires  an 
examination  of  the  accuracy  required  by  users  of  the  wave  information 
and  the  accuracy  available  from  the  various  models.  The  high  cost  of 
dike  construction  for  dredged  material  retention  and  the  high  cost  of 
coastal  zone  construction  certainly  justify  selection  of  the  most  ac- 
curate model. 

25.  The  major  disadvantages  in  applying  the  SMB  technique  are 
linked  to  its  use  of  a single  parameter  to  represent  sometimes  complex 
phenomena.  If  the  wind  field  were  a known  constant  over  a uniform  fetch 
and  if  the  waves  were  truly  monochromatic,  the  SMB  technique  should  be 
an  excellent  method  for  estimating  wave  heights.  The  reduction  in  use- 
fulness and  accuracy  of  the  SMB  technique  for  the  Great  Lakes  is  related 
to  deviations  from  these  assumptions.  Variability  in  the  shoreline 
geometry  of  the  Great  Lakes  and  in  wind  velocities  and  directions  during 
the  passage  of  a storm  can  introduce  considerable  error  in  the  attempt 
to  simplify  the  winds  and  fetches  to  a form  commensurate  with  assump- 
tions inherent  in  SMB  hindcasts. 

26.  In  a wave  spectrum  specified  at  a given  time,  the  energy 
density  may  be  fetch-limited  for  some  frequencies  and  duration-limited 
for  others.  With  the  SMB  technique,  however,  waves  can  be  treated  only 
as  fetch-limited  or  duration-limited. 

27.  Recent  field  studies  into  the  growth  of  waves  tend  to  point 
out  this  problem  of  fetch  and  duration  equivalence.  According  to  the 
most  recent  deepwater  curves  published  for  the  significant  wave  method 

C 

(SPM')  for  a 30-knot  wind,  a 50-mile  fetch  is  equivalent  to  about 


2h 


a 7-hr  duration.  In  other  words,  waves  for  a 50-mile  fetch  will  he 


duration-limited  as  long  as  the  wind  has  blown  less  than  7 hr.  A field 

29 

study  by  Barnett  and  Wilkerson  demonstrated  that  waves  with  periods 
as  high  as  8.5  sec  were  fully  developed  within  50  miles  of  a coast  for 
winds  of  30  knots  blowing  straight  offshore  (Figure  lU).  The  group 
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DISTANCE  FROM  COAST  (NAUTICAL  MILES) 

Figure  lU.  Growth  of  frequency  components  of 
naves  for  an  offshore  wind  of  30  knots  (after 
Barnett 30 ) 

velocity  of  these  waves  is  around  15  mph,  so  the  total  possible  time 
under  the  action  of  the  wind  is  only  about  3.3  hr.  Thus,  the  waves 
would  become  fetch-limited*  after  only  3.3  hr,  not  7 hr  as  given  by  the 
SMB  method. 

28.  Because  of  these  apparent  inadequacies  in  the  SMB  method  for 
conditions  typical  of  the  Great  Lakes,  an  investigation  into  the  poten- 
tial of  numerical  models  was  initiated.  Such  models  have  been  in  use 

31 

since  the  early  1960's  (Baer  ) and,  at  various  stages  of  development. 


* This 
lower 

assumption  that  higher  frequencies  are  fully  developed  before 
frequencies  also  was  validated  in  the  same  field  study. 
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have  been  proven  capable  of  providing  good  estimates  of  wave  heights. 
These  models  differ  from  the  significant  wave  model  in  that  they  use  a 
system  of  grid  points  on  the  water  surface  for  calculations  of  energy 
transfer  into  and  out  of  the  wave  spectrum.  In  this  way,  they  can 
achieve  a much  more  detailed  representation  of  the  shoreline  and  wind- 
field  geometries.  Also,  since  the  model  calculates  energy  transfers  for 
small  increments  of  time,  the  wind  field  in  a numerical  model  can  be 
specified  more  exactly  in  time  than  the  wind  field  in  a significant  wave 
calculation. 

29.  Several  numerical  models  are  presently  available.  Basically, 
all  of  these  using  theoretical  representations  of  energy  transfer  mech- 
anisms compute  energy  spectra  at  each  grid  point.  The  ability  to  retain 
information  from  each  grid  point  is  certainly  an  advantage  of  these 
models.  Thus,  what  might  require  several  separate  hindcasts  with  an  SMB 
approach  can  often  be  done  in  one  step  with  a numerical  model. 

30.  Only  two  numerical  models  of  wind-generated  waves  have  been 

used  extensively  in  the  United  States — one  developed  at  New  York  Uni- 

13  32  19 

versity  by  Pierson  and  Moskowitz,  Inoue,  and  Cardone,  and  another 

30 

developed  at  Scnpps  Institute  of  Oceanography,  primarily  by  Barnett. 
These  two  models  are  discussed  in  more  detail  in  Appendix  B.  A third 
model,  the  French  spectro-angulaire  model,  also  has  been  used  in  the 
United  States,  but  will  not  be  examined  here  since  it  is  similar  in  many 
respects  to  the  model  developed  at  Scripps. 

31.  The  primary  difference  between  these  models  is  the  form  of 

the  energy  transfer  into  the  waves  attributed  to  different  theoretical 

mechanisms.  The  "NYU"  model  essentially  considers:  (a)  a linear  source 

due  to  the  interaction  of  the  water  surface  with  turbulent  pressure 

fluctuations  in  the  atmosphere  (Eckart,^  Phillips^);  (b)  a nonlinear 

mechanism  due  to  a resonant  interaction  between  waves  and  the  mean  wind 
35 

profile  (MilesJ^);  and  (c)  a wave-breaking  term  related  to  the  fully 

13 

developed  spectral  form  proposed  by  Pierson  and  Moskowitz.  The 

30 

model  developed  by  Barnett,  in  addition  to  the  Miles'  and  Phillips' 
mechanisms  and  wave  breaking,  considers  the  effect  of  wave-wave  inter- 

n/T  O 

actions  on  the  development  of  a spectrum  (Hasselmami  ~ ) . 
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32.  Even  though  the  basic  research  needed  to  clarify  the  energy 

transfer  mechanisms  is  far  from  complete,  both  models  have  been  tested 

39 

successfully  in  previous  studies.  Lazonoff  et  al.  documents  compari- 
sons between  hindcast  wave  heights  using  the  "NYU"  model  and  observed 
waves  on  the  Mediterranean  Sea.  For  this  set  of  comparisons,  80  percent 
of  the  observed  and  hindcast  wave  heights  were  within  3 ft  of  each  other. 

The  Barnett  model  has  not  been  as  widely  tested  in  deep  water.*  Fig- 

30 

ure  15  shows  a single  time  series  comparison  made  by  Barnett  for  an 
area  in  the  Atlantic. 
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Figure  15.  Hindcasted  significant  wave  heights  for 
l6-l8  Dec  1959.  using  Barnett  model  with  60-nautical- 
mile  grid  "0"  and  120-nauti cal -mile  grid  "0"  as  com- 
pared with  observed  wave  heights  from  the  OWS 
"Weather  Reporter" 


* Fleet  Numerical  Weather  Central  in  Monterey,  California,  has  now  com- 
pleted adapting  a version  of  the  wave-wave  interaction  model  to  their 
computer  systems.  Preliminary  results  from  hindcasts  in  the  Gulf  of 
Alaska  region  have  compared  extremely  well  with  wave  observations 
(personal  communication  from  Professor  Warren  Thompson  of  the  Naval 
Postgraduate  School,  Monterey,  California). 
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33.  In  previous  tests  on  Lakes  Erie  and  Ontario,  several  se- 
quences of  wave  spectra  were  selected  for  comparison  with  hindcasts  hy 
the  "NYU”  and  Barnett  models.  The  corresponding  significant  wave  heights 
for  the  periods  from  which  these  spectra  were  taken  started  at  some  low 
level  (about  1 to  2 ft),  progressed  to  some  maximum  or  maxima  (greater 
than  5 ft),  and  then  returned  to  near  the  initial  level.  With  these 
time  series,  the  two  important  aspects  in  the  model  were  (a)  the  compari- 
son with  the  complete  hydrograph  and  (b)  the  comparison  between  observed 
maxima  and  hindcast  maxima.  In  this  study,  the  behavior  of  both  the 
time  series  and  individual  maxima  hindcast  with  the  "NYU"  model  was 
found  to  be  more  erratic  than  the  hindcast  with  the  wave-wave  interaction 
model.  Also,  the  nondimensional  wave  heights  predicted  by  the  "NYU" 

model  were  not  consistent  with  recent  field  observations  such  as  those 

40  .. 

of  Hasselmann  et  al.  Figures  16  and  17  present  two  comparisons  be- 
tween sequences  of  wave  heights  hindcast  by  the  latter  model.  Figure  18 
gives  the  distribution  of  observed  maxima  relative  to  hindcast  maxima. 
Additional  results  from  model  calibrations  are  given  in  Appendix  B. 

34.  The  wind  input  to  these  test  runs  on  the  computer  consisted 
of  3-hr  winds  taken  from  wind  data  on  magnetic  tapes  for  Buffalo  and 
Rochester  on  Lake  Ontario  and  Toledo,  Cleveland,  Erie,  and  Buffalo  on 
Lake  Erie.  These  winds  were  input  into  the  program  directly,  and  all 
conversions  and  interpolations  of  the  wind  field  are  performed  within 
the  computer  program.  This  means  that  subjectivity  in  the  wind  input 
used  in  this  calibration  was  held  to  an  absolute  minimum. 

35.  At  the  time  of  the  calibration  tests  for  Lake  Huron,  there 
were  no  deepwater  wave  gage  data  available  to  this  study.  Consequently, 
no  additional  calibrations  are  included  in  this  report. 

36.  The  rms  error  for  the  hindcasts  shown  in  Figure  18  is  about 
1.5  ft.  The  tendency  of  the  model  to  overpredict  in  short-fetch  condi- 
tions, as  noted  in  Reports  1 and  2,  is  not  as  pronounced  as  indicated  by 
the  smaller  number  of  hindcasts  considered  in  those  earlier  reports. 

The  underprediction  during  the  summer  period  appears  to  remain,  however, 
and  is  most  likely  due  to  the  lack  of  resolution  of  the  small-scale 
summer  circulation  features. 
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Figure  l6.  Comparison  of  time  series  of  significant  wave  heights 
observed  at  Muskegon  for  storm  of  29  October  1965  and  wave 
heights  hindcast  by  numerical  model  for  same  storm 
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Figure  17.  Comparison  of  time  series  of  significant  wave  heights  observed  at 
Point  Pelee  for  storm  of  7 April  1973  and  wave  heights  handcast  by  numerical 

model  for  same  storm 


Figure  18.  Comparison  of  maximum  significant  wave 
heights  for  individual  storms  as  recorded  by  gages 
(Point  Pelee,  Lake  Erie;  and  Toronto,  Cobourg,  and 
Main  Duck,  Lake  Ontario)  and  significant  wave 
heights  hindcast  by  the  numerical  model  using 
Barnett's  parameterizations 

37.  Another  possible  contributor  to  the  error  in  short  fetch  con- 
ditions is  the  relative  importance  of  mesoscale  mhomogeneities  in  the 
wind  field  for  these  cases.  In  the  short  fetch  situation,  the  waves  are 
propagating  in  an  offshore  direction;  in  the  summer  storms,  actual  max- 
ima are  only  about  5 to  6 ft.  Neither  of  these  cases  are  critical  to 
the  practical  problem  of  forecasting  extreme  wave  heights  along  the 
coastlines  of  the  Great  Lakes.  The  agreement  between  hindcast  wave  max- 
ima and  observed  wave  maxima  for  cases  involving  fetches  over  20  miles 
and  for  waves  propagating  toward  shore  was  extremely  good.  All  of  the 
hindcast  maxima*  were  within  1.5  ft  of  observed  maxima.  The  rms  error 


* In  some  cases,  the  wave  gages  did  not  function  part  of  the  time  dur- 
ing the  high-wave  conditions.  For  these  instances,  comparisons  were 
made  between  wave  heights  at  the  time  of  the  highest  gage  record. 
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in  estimating  peak  wave  heights  for  this  set  of  conditions  is  about  1 ft. 

38.  The  agreement  between  the  wave-wave  interaction  model  and 
observed  waves  unequivocally  confirms  the  capability  of  state-of-the-art 
numerical  wave  models,  coupled  with  adequate  wind  data,  to  produce  re- 
sults of  sufficient  accuracy  for  most  engineering  applications.  One 
point  that  should  be  brought  out,  however,  is  the  lack  of  observed  wave 
heights  above  15  ft  with  which  to  extend  the  calibration  curve.  If 
strong  nonlinear  effects  are  present  in  waves  beyond  the  range  of  the 
calibration,  these  would  be  impossible  to  detect  with  the  present  data. 
However,  the  numerical  wave  models  have  been  calibrated  at  ocean  sites 
for  waves  considerably  larger  than  15  ft.  Thus,  unless  there  is  a large 
error  associated  with  the  friction  term,  there  is  little  reason  to  ex- 
pect a marked  increase  in  error  for  wave  heights  in  excess  of  15  ft  hind- 
cast  with  the  present  model. 
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PART  IV:  THE  ANALYSIS  OF  EXTREME  WAVE  HEIGHTS 


r~ 


39-  In  Part  III,  a capability  to  reproduce  wave  heights  from  in- 
dividual storms  was  demonstrated.  This  information  is  valuable,  but  it 
is  the  statistical  synthesis  of  this  information  into  the  estimation  of 
probabilities  of  extremes  that  is  most  important  in  design  and  planning 

applications.  The  analytical  treatment  of  extremes  is  well  developed 
/ l+l  1+2  1+3, 

(Jenkinson,  Gumbel,  and  Gringorten  ) and  can  be  used  very  effec- 
tively for  values  not  extrapolated  too  far  beyond  the  total  length  of 
record. 

1+0.  A major  facet  of  this  study  has  been  to  obtain  long  records 
of  reliable  wind  data.  For  this  purpose,  several  problems  had  to  be 
addressed.  First,  only  about  26  yr  of  the  wind  data  were  available  on 
magnetic  tape;  the  remainder  were  handcopied  from  autograph  forms  avail- 
able at  the  National  Climatic  Center,  Asheville,  North  Carolina.  Second, 
wind  data  are  not  available  at  a consistent  height  or  a consistent  loca- 
tion for  the  entire  period  of  record;  thus  a detailed  analysis  of  the 
wind  profile,  as  discussed  in  Appendix  A,  is  necessary.  Results  of  these 
analyses  are  transformations  of  the  wind  from  one  height  to  another. 

For  the  simple  case  in  which  the  wind  to  be  transformed  can  be  con- 
sidered totally  within  a boundary  layer  characteristic  of  land,  the  one- 
seventh  power  law  can  be  invoked.  These  winds  are  reduced  to  a common 
anemometer  height  (20  ft)  and  transformed,  using  the  empirical  results 
obtained  in  Part  II,  into  winds  over  the  lake.  For  stations  located 
within  a mile  or  so  of  the  lakeshore,  winds  blowing  off  the  lake  still 
retain  characteristics  of  the  marine  boundary  layer  (particularly  for 
anemometer  heights  typical  of  older  Great  Lakes  weather  stations,  200- 
300  ft).  These  winds  are  transformed  directly  to  winds  over  the  lake 
using  a theoretical  treatment  of  meteorological  boundary  layers.  In  this 
way,  a single,  quasi-homogeneous  population  of  winds  estimated  over  the 
lake  is  obtained.  Appendix  A details  the  techniques  used  to  formulate 
these  transformations. 

1+1.  In  order  to  obtain  design  information  for  all  directions  of 
wave  approach  for  coastal  structures,  the  wind  data  were  stratified  by 


33 


a 1+5-deg  direction  class  (i.e.,  N,  NW,  W,  etc.).  They  were  further  sub- 
divided by  separating  the  winds  into  seasonal  categories  (January-March, 
April-June,  July- Sept ember,  October-December ) . This  was  accomplished  to 
facilitate  treatment  of  seasonal  variations  in  water  level  and  icing  con- 
ditions on  the  Great  Lakes.  Seasonal  wind-speed  maxima  (in  each  direc- 
tion) for  each  wind  station  were  determined;  the  storms  during  which 
these  wind  speeds  occurred  were  keypunched.*  To  reduce  computer  time, 
only  the  larger  half  of  these  storms  (in  terms  of  the  wind-speed  crite- 
ria) was  input  into  the  numerical  model.  For  each  storm  (at  1-hr  inter- 
vals), two-dimensional  spectra  for  all  grid  points  adjacent  to  shore 
were  written  onto  magnetic  tapes  and  saved.  Maximum  wave  heights**  by 
direction  are  output  in  the  form  of  punched  cards.  It  is  the  latter 
information  that  is  used  in  the  calculation  of  wave  height  recurrence 
intervals . 

1+2.  Since  only  seasonal  maxima  were  computed  for  wave  heights, 
the  appropriate  method  to  use  in  estimating  recurrence  intervals  is 

1+2 

taken  from  the  asymptotic  theory  of  extremes.  In  their  texts,  Gumbel 

1+3 

gives  an  excellent  discussion  of  these  techniques,  and  Gringorten  pre- 
sents a good  demonstration  of  the  application  of  these  techniques. 
Basically,  the  asymptotic  forms  are  solutions  to  the  equation 

Fn(x)  = F(a  x + b ) (14) 

n n 


* A minimum  of  6 hr  before  the  maximum  and  3 hr  after  the  maximum  were 
included  in  the  storm  to  ensure  that  the  peak  wave  conditions  were 
calculated. 

**  Significant  wave  heights  for  nearshore  conditions  are  estimated  from 
the  two-dimensional  spectra  at  the  grid  points  by  the  relationship 
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where  fp  and  f^  are  the  lowest  and  highest  frequencies  considered 
in  the  numerical  model,  0p  and  0p  are  limiting  angles  for  waves 
propagating  toward  the  shore,  and  F(f,0)  is  the  two-dimensional 
spectrum. 


where  a and  b are  constants,  which  are  functions  of  sample  size 
n n r 

n , and  F is  the  cumulative  probability  function  of  x . This  theory 

is  based  on  the  premise  called  the  Stability  Postulate,  which  states 

that  the  distribution  of  the  largest  value  in  Nn  observations  will 

approach  the  same  asymptotic  expression  as  the  distribution  of  the 

largest  value  in  N samples  of  size  n . The  only  three  solutions 

possible  for  Equation  U are  termed  the  Fisher-Tippett  Types  I,  II,  and 

III  distributions.  Each  solution  is  derived  from  different  constraints 

111 

on  a and  b ; however,  Jenkinson  showed  that  all  three  of  these 
n n 

distributions  could  be  written  in  a common  form  as 


x = x 

o 


(5) 


where 

xq  = parameter  of  extremal  distribution 
a = parameter  of  extremal  distribution 
e = Napierian  base  (2.71828...) 

k = parameter  denoting  curvature  in  extremal  distribution 
y = reduced  variate  in  distribution  of  extremes 
Figure  19  shows  the  form  of  Fisher-Tippett  Types  I,  II,  and  III  dis- 
tributions when  k = 0 , k < 0 , and  k > 0 , respectively.  It  is  also 
important  to  note  in  this  figure  the  bounds  on  the  different  types: 

(a)  Type  I is  unbounded;  (b)  Type  II  has  a lower  limit;  and  (c)  Type  III 
has  an  upper  limit.  Unfortunately,  there  is  no  reliable,  objective 
method  for  determining  which  distribution  one  is  dealing  with  when  given 
only  a single  set  of  sample  characteristics.  Arguments  based  on  the 
stability  of  extrapolations  and  the  ease  of  application  usually  favor 
the  Fisher-Tippett  Type  I distribution  when  curvatures  are  relatively 
low. 

1*3.  In  the  analysis  of  wave  height  extremes,  the  first  step  was 
to  plot  the  wave  heights  computed  in  the  wave  model  against  the  log- 
arithm of  recurrence  intervals*  estimated  from  the  formula 

* The  terms  "recurrence  interval"  and  "return  period"  are  used  inter- 
changeably in  this  report. 
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THE  SCALE  OF  X IS  IN  UNITS 
OF  THE  PARAMETER  a , TAK- 
ING X0  AS  ORIGIN,  a IS  THE 
SLOPE  OF  THE  X,Y  CURVE 
AT  Y=0,  WHERE  X=X0 

GENERAL  EQUATION 
X=X0+a(t£±Z) 

WHICH  SIMPLIFIES,  FOR  h=0 
TO  x=x0+«y  (TYPE  r) 


Figure  19-  Schematic  representation  of  the  three 
Fisher-Tippett  asymptotic  distributions  on  x-y 
diagrams  (after  Jenkinson^) 


T(H) 


1 

1 - F(H) 


(6) 


where 

T = return  period 
H = wave  height 

F(h)  = cumulative  probability  function  of  H 
A typical  result  of  such  plots  is  shown  in  Figure  20.  The  curvature  at 
the  bottom  of  these  graphs  is  due  to  the  inclusion  of  samples  from 
storms  other  than  the  seasonal  maxima  for  each  direction  and  the  in- 
clusion of  only  the  larger  half  of  the  annual  maxima  in  the  numerical 
hindcasts. 
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RECURRENCE  INTERVALS  FOR  ANGLE  CLASS  2 
HARBOR  BEACH.  MICH.  GRID  LOCATION  16.  11 
LAT  = 43.88  LONG  = 82.52 


J I I I 

0 ’ 0 20  30  40  50  60 
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Figure  20.  Sample  plot  of  log  T versus 
significant  wave  height  as  calculated  from 
Equation  6 

UU.  The  upper  portion  of  these  graphs  represents  a region  iden- 
tical to  that  which  would  be  calculated  from  an  analysis  of  all  of  the 
annual  maxima.  As  discussed  in  Appendix  D,  this  section  of  the  curve 
can  be  used  to  reconstruct  estimates  of  parameters  for  the  entire 
curve.  These  results  are  based  on  calculations  designed  to  optimize  the 
estimation  of  large  return  periods  and  hence  should  not  be  used  to 
estimate  return  periods  of  less  than  a year.  Another  factor  to  consider 
in  applying  asymptotic  methods  to  small  recurrence  intervals  is  empha- 
sized  by  Dalrymple.  In  a series  of  annual  maxima,  the  recurrence 
interval  is  defined  as  the  average  interval  in  which  a wave  height 
greater  than  or  equal  to  a given  size  will  recur  as  an  annual  maximum. 

A partial-duration  series  must  be  used  to  define  the  average  interval 
between  wave  heights  of  a given  size  regardless  of  their  relationship  to 
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other  storms  of  the  same  year.  Relationships  between  recurrence  in- 
tervals defined  by  the  two  different  techniques  have  been  derived  by 
Langbein^  and  Chow*4"*"  (Table  5). 

45.  Figure  21  shows  the  locations  of  grid  points  used  to  define 
"deepwater"*  wave  heights  along  the  U.  S.  shoreline  of  Lake  Michigan. 
These  points  were  chosen  to  be  sufficiently  removed  from  the  effects  of 
shallow  water  and  local  coastline  geometry  to  peimit  them  to  be  represen- 
tative of  the  segments  of  coast  in  which  they  are  centered.  Also,  the 
effects  of  storm  surges,  seiches,  and  other  fluctuations  in  water  level 
become  increasingly  important  as  the  coast  is  approached.  Whereas  sta- 
tistics specified  off  the  coast  can  be  considered  independent  of  small 
changes  in  lake  level,  waves  and  water  levels  closer  to  shore  must  be 
superimposed  in  order  to  obtain  calculations  of  wave  heights.  For  storm 
waves  calculated  during  this  study,  the  two-dimensional  spectra  can  be 
retrieved  from  magnetic  tape  and  input  into  shallow-water  wave  routines 
along  with  detailed  specifications  of  water  levels  during  the  storm. 
Although  this  is  adequate  for  individual  storms,  it  might  be  more  ad- 
vantageous to  investigate  the  covariance  between  the  wave  spectra  and 
short-term  fluctuations  of  the  lake  level.  Thus,  a statistical  model 
could  be  derived  that  could  provide  nearshore  wave  information  at  mini- 
mal cost. 

46.  Ice  cover  significantly  influences  winter  wave  climates  on 
all  of  the  Great  Lakes  with  the  possible  exception  of  Lake  Ontario.  Two 
factors  must  be  considered  in  examining  the  effect  of  ice  cover  on  re- 
ducing wave  heights  incident  on  shores  or  structures:  (a)  reduction  of 
fetch  lengths  for  wave  generation,  and  (b)  protection  from  waves  af- 
forded by  ice  cover  near  the  shore.  Since  ice  typically  forms  around 
the  periphery  of  a lake  before  extending  toward  the  center,  the  second 
of  these  two  factors  is  more  important  for  most  nearshore  sites.  As  a 
reasonable  first  approximation,  meteorological  conditions  that  generate 


* These  depths  are  not  greater  than  one-half  the  wavelength  of  some  of 
the  waves  on  the  Great  Lakes,  but  they  are  in  locations  where  refrac- 
tion and  frictional  effects  would  be  minimally  affected  by  changes  in 
the  water  level  of  up  to  several  feet. 
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Figure  21.  Locations  of  grid  points  for  which  two-dimensional 
spectra  are  retained  on  magnetic  tape  and  for  which  design 
wave  information  is  listed  in  Appendix  h 
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large  waves  in  the  winter  (individual  extratropical  cyclones)  can  be 
considered  independent  of  meteorological  conditions  responsible  for  ice 
formation  (long-term,  large-scale  circulation  patterns).  A much  more 
detailed  study  would  be  required  to  establish  relationships  between 
storm  probabilities  and  large-scale  circulation  patterns,  and  the  in- 
crease in  the  accuracy  of  estimates  of  recurrence  intervals  could  very 
well  be  marginal.  Consequently,  the  assumption  of  independence  seems 
justifiable  in  this  study.  With  this  assumption,  the  probability  of  a 
particular  wave  height  (p(H))  occurring  at  a site  can  be  written  as  the 
product  of  the  probability  of  its  occurence  estimated  from  hindcasts  not 
considering  ice  cover  (p’(H))  and  the  probability  of  there  being  no  ice 
at  the  site  (p(e)) 

p(H)=  p'(H)  p(e)  (?) 

1*7.  Thus  the  recurrence  interval  for  conditions  with  possible  ice 
cover  can  be  expressed  as 

T'  = T 4>  (8) 

where  T is  the  recurrence  interval  for  waves  with  no  ice  cover  (as 
determined  by  the  computer  program)  and  $ is  defined  as 

(9> 

Equation  9 provides  a simple  method  of  converting  the  generalized  wave 
statistics  into  statistics  applicable  to  specific  locations  around  a 
lake. 

1*8.  If  the  effect  of  fetch  reduction  appears  necessary  for  rea- 
sonable wave  height  estimates  at  a particular  site,  a transformation 
can  be  applied  such  as 

Hl  = g(H,  F,  Fl)  (10) 

where 

= actual  wave  height  as  affected  by  fetch  reduction  due  to 
ice  cover 


1*0 


g = function  of  H , F , and  F 

L 

F = land-bounded  fetch 

F = fetch  during  partial  ice  cover 
Li 

Since  the  original  wave  height  is  calculated  with  the  complete  numerical 
program,  a simple  relationship  based  on  dimensional  considerations  can 
be  applied  to  reproduce  the  approximate  variation  of  wave  height  with 
fetch  without  introducing  too  much  error.  Hence,  the  curves  given  in 
SPM  for  wave  growth  with  fetch  (ignoring  the  duration  of  winds)  can 
be  used  to  scale  wave  heights  predicted  with  no  ice  cover  into  wave 
heights  with  reduced  fetches.  Figure  22  presents  these  curves  for 


Figure  22.  Growth  of  significant  wave  height  as  a function  of  wind 

speed  and  fetch  (after  SPM5) 

fetches  from  20  to  200  miles.  For  example,  if  the  original  fetch  were 
100  miles,  these  curves  indicate  that  a wind  velocity  of  approximately 
28  knots  would  be  required  to  generate  a significant  wave  height  of 
10  ft.  If  ice  cover  reduced  the  fetch  to  50  miles,  the  significant 
wave  height  for  this  velocity  is  expected  to  be  about  7 ft.  Hence,  the 
probability  of  exceeding  a 10-ft  wave  height  becomes  the  probability  of 


exceeding  a 7-ft  wave  height.  By  using  such  a technique,  the  recurrence 
interval  can  be  written  as 

T'(H)  = T(H')4’  (11) 

where  H is  the  wave  height  into  which  H'  is  transformed  by  fetch 
reduction. 

1*9*  In  Appendix  E,  results  of  numerical  computations  of  wave 
heights  with  recurrence  intervals  of  5,  10,  20,  50,  and  100  yr  are  given 
along  with  the  control  bands  for  each  point  in  Figure  21.  These  are 
tabulated  as  (a)  a function  of  season  and  approach  direction,  (b)  a 
total  of  all  seasonal  probabilities  for  each  approach  direction,  and 
(c)  a total  of  all  direction  and  seasonal  probabilities.  The  signifi- 
cant periods  associated  with  these  wave  heights  also  are  included  in 
Appendix  E. 


PART  V:  DISCUSSION  OF  RESULTS 


50.  This  study  represents  a synthesis  of  three  phases  of  wave 
hindeasting: 

&.  Estimation  of  winds  over  the  water, 
b.  Operationalization  of  a wave  model. 

c_.  Calculation  of  return  period  information  from  program 
outputs. 

Each  phase  of  the  study  was  treated  in  detail,  since  errors  in  any  phase 
could  alter  the  results  dramatically.  As  discussed  in  Appendix  D,  the 
total  error  in  wave  hindcasts  relative  to  the  variability  of  annual  max- 
ima is  sufficiently  small  that  the  distribution  of  waves  from  the  hind- 
casts  should  not  be  affected  by  it;  and  as  seen  in  the  calibration 
results  (Part  III),  the  waves  estimated  from  the  wave  model  do  not  ap- 
pear to  contain  any  determinable  bias.  Thus,  the  primary  source  of 
error  in  the  estimates  of  return  periods  is  due  to  expected  sampling 
variability  and  is  inversely  proportional  to  the  square  root  of  the 
number  of  years  in  the  sample. 

51.  It  is  important  for  any  user  of  information  reported  herein 
to  bear  in  mind  that  the  wave  heights  presented  in  Appendix  E refer  to 
deepwater  conditions.  The  accuracy  of  the  results  appears  to  be  quite 
good  for  these  conditions,  and  they  can  be  applied  as  a general  approxi- 
mation to  wave  heights  along  the  coast.  However,  when  detailed  informa- 
tion for  design  and  planning  purposes  is  required,  it  is  imperative  that 
the  two-dimensional  spectra  be  retrieved  from  magnetic  tapes  and  used 

in  conjunction  with  the  proper  shallow-water  transformations.  Although 
significant  wave  techniques  can  provide  some  idea  of  the  actual  trans- 
formations, shallow-water  wave  heights  estimated  from  these  methods  tend 
to  be  biased  toward  higher  waves  (Freeman  et  al.  °).  A better,  less 

biased  estimate  of  shallow-water  wave  heights  probably  could  be  obtained 
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from  numerical  shallow-water  routines  such  as  developed  by  Collins. 
Substantial  evidence  is  now  becoming  available  that  indicates  a marked 
energy  loss  in  moving  from  deep  to  shallow  water.  This  loss  of  energy, 


b3 


if  substantiated,  could  provide  a significant  reduction  in  the  wave 
heights  listed  in  Appendix  E. 

52.  As  a final  point,  it  should  be  emphasized  that  the  informa- 
tion presented  in  this  report  constitutes  only  a small  portion  of  that 
available  from  the  complete  set  of  hindcast  data  stored  on  magnetic 
tapes.*  For  many  purposes,  a natural  choice  for  design  considerations 
formed  by  navigation  constraints  does  not  always  coincide  with  seasons 
chosen  for  use  in  this  study.  However,  the  basic  data  set  can  be  re- 
processed within  any  arbitrary  seasonal  divisions.  Other  information, 
such  as  characteristic  durations  of  wave  heights,  estimates  of  nearshore 
currents,  and  estimates  of  individual  wave  height  probabilities,  can  all 
be  obtained  from  analyses  of  the  hindcast  data. 


* These  two-dimensional  spectra  are  stored  at  1-hr  intervals  during  each 
hindcast  storm  for  all  shoreline  points  and  selected  offshore  points. 
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Table  1 

Wind  Analyses* 


Comparisons  with  Observed  Winds 

Number  of 

Data  Pairs 

Correlation 

Coefficient 

Bret schneider  winds  versus  10-m  winds 

36 

0.63 

Jacobs  7.5-m  winds  versus  7-5-m  winds 

1+3 

0.55 

Jacobs  19-5-m  winds  versus  l6-m  winds 

1+9 

0.37 

Richards  winds  versus  l6-m  winds 

1+1+ 

0.36 

Richards  winds  versus  10-m  winds 

36 

0.21+ 

* From  Reference  7- 


Table  2 

Estimates  of  Wind  Ratios  "R"  for  the  Great  Lakes* 


1 7 

Richards 

15 

a 

Richards 

Dragert , ; 

Hunt 

Lemire 

Extension 

McIntyre- 

Jan 

1.96 



Feb 

1.9!+ 

— 

Mar 

1.88 

1.88 

— 

Apr 

1.81 

1.81 

— 

Spring 

1.35 

1.38 

May 

1.71 

1.71 

— 

Jun 

1.31 

1.31 

— 

Jul 

1.16 

1.16 

— 

Aug 

1.39 

1.39 

— 

Fall 

1.82 

1.87 

Sep 

1.78 

1.78 

— 

Oct 

1.99 

1.99 

— 

Nov 

2.09 

— 

Dec 

1.98 

— 

Average 

1.56* ** 

1.63** 

1 .66+ 

1.56+ 

* From  Reference  18. 

**  Navigation  season, 
t Annual . 


Table  3 


I 


Correlation  Between  Winds  Over  Land  and  Wind  Observations 
by  Ship  Stratified  by  Stability  Class  and  Number  of 
Observations  per  Three  Hours 


Stability 

Number  of  Observations  per  3 

hr 

Class 

1 

2 

3 

1+ 

1 

0.780(38) 

0.825(17) 

0.809(9) 

2 

0.528(119) 

0.695(50) 

0.717(30) 

0.707(13) 

3 

0.1+80(207) 

0.590(92) 

0.610(39) 

0.715(11) 

1+ 

0.589(227) 

0.698(83) 

0.898(31+) 

0.920(11) 

5 

0.370(97) 

0.631(36) 

0.672(15) 

0.825(3) 

6 

0.358(31) 

0.1+26(15) 

0.  59l+  ( 6 ) 

7 

0.730(10) 

0.097(1+) 

— 

— 

Note:  Number 
relation 

in  parentheses 
is  based. 

represents  number 

of  samples 

on  which  cor- 

Table  1+ 

RMS  Error  in  Wind  Estimates 
Over  the  Great  Lakes 


Land 

RMS 

Wind  Speed 

Error 

knots 

knots 

0-3 

5.2 

3-6 

5.9 

6-9 

5.3 

9-12 

1+.6 

12-15 

l+.o 

15-18 

3.9 

18-21 

3.9 

21-21+ 

3.6 

21+-27 

3.5 

Table  5 

Relationship  Between  Recurrence  Intervals 


as  Calculated  from  Partial-Duration 


Series  and  from  Annual  Maxima 


Partial-Duration 

Series 


Annual  Floods 


10.00 

20.00 

50.00 


10.50 

20.50 

50.50 


100.00 


100.50 


APPENDIX  A:  METHODOLOGY  OF  OVERLAKE  WIND  ESTIMATION 


1.  In  deep  water,  the  generation  of  waves  by  the  wind  can  be 
given  in  functional  forms  as  (Kitaigorodskii'’^) : 

F(oo)  = G(w,T,X)  (Al) 

where 

F(w)  = one-dimensional  spectral  energy  and  a function  of  angular 
frequency  w 

G = function 

T = nondimens ional  time  factor 
X = nondimensional  fetch  factor 

Both  T and  X are  in  turn  functions  of  the  friction  velocity 

X = *|  (A2) 

u* 

T = ^ (A3) 

u# 

where 

x = actual  fetch 
g = gravitational  acceleration 
us  = friction  velocity 
t = wind  duration 

The  dependence  of  wave  growth  on  the  friction  velocity  necessitates  a 
knowledge  of  the  wind  profile  rather  than  the  wind  at  only  one  level . 

2.  Prior  to  the  mid-1960's,  the  height  of  wind  instruments  above 
ground  varied  considerably  from  station  to  station  (20  to  > 300  ft). 
Instrument  heights  and  station  locations  at  most  cities  were  changed 
several  times. 

3.  The  methods  for  transforming  the  raw  wind-speed  data  into 
inputs  to  the  wave  hindcasting  program  will  be  considered  in  two  parts : 
(a)  data  taken  prior  to  19^8,  and  (b)  data  taken  after  19^8.  This  date 
is  chosen  as  a separator  because  the  National  Weather  Service  began 

t Raised  numbers  refer  to  similarly  numbered  items  in  "References"  on 
pp  k5-k9  at  end  of  main  text. 


consistently  recording  data  at  all  first-order  airport  weather  stations 
in  19^8. 

Post-19^8  Wind  Transformations 

If  the  assumption  is  made  that  the  wind  profile  at  an  instru- 
ment site  is  logarithmic  in  shape,  the  following  relationship  can  be 
used  to  relate  two  wind  speeds  (u^  , u^)  taken  at  different  elevations 
(Zx  , Z2): 


Over  land  where  a neutral  profile  can  usually  be  assumed  for  hourly  aver- 
age conditions.  Equation  AU  gives  a good  estimate  of  the  wind,  if  only 
the  instrument  height  varies  and  not  location  of  the  instrument  tower. 

5.  Changes  in  the  location  of  the  wind  instruments  generally  re- 
sult in  changes  in  roughness  for  the  local  boundary  layer  unless  the 
distances  are  small  and  the  environmental  characteristics  homogeneous 
over  a large  area.  Hence,  the  one-seventh  power  law  may  produce  errors. 
To  test  the  effect  of  movement  of  stations,  the  subsequent  approach  was 
used. 

6.  For  every  season  and  station,  the  wind-speed  distribution  func- 
tion was  estimated  for  each  period  where  instrument  location  and  height 
were  constant.  Recurrence  intervals  were  calculated  and  plotted  (Fig- 
ure Al)  and  the  linear  least-square  estimate  of  the  relationship  between 
the  logarithm  of  the  recurrence  interval  and  wind  speed  was  determined 
for  winds  in  the  15-  to  30-knot  velocity  range.* 

7.  To  calculate  the  transformation  to  the  standard  station  height 
since  19&5  (20  ft),  the  equations  between  recurrence  interval  and  wind 


* This  range  was  selected  on  the  basis  of  there  being  a sufficient  fre- 
quency of  occurrence  of  these  wind  speeds  to  provide  a large  sample 
for  the  regression  calculations.  The  nonlinear  effect  of  deviations 
due  to  sampling  in  higher  velocities  was  thus  avoided. 
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speed  for  each  period  are  expressed  as 


106  T*  = a*  + V 


(A5) 


where 


log  T = a + B u 
s s s s 


= recurrence  interval  at  elevation  above  20  ft 
T = recurrence  interval  at  elevation  of  20  ft 


(A6) 


= regression  coefficients  relating  wind  speeds  at  different 
elevations 
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If  population  characteristics  of  wind  speed  are  assumed  to  be  reasonably 
constant,  the  effect  of  changes  in  anemometer  location  can  be  estimated 
by  equating  recurrence  intervals.  Combining  Equations  A5  and  A6  gives 


and  hence 


a£  + V = l0g  T*  = log  Ts  = as  + BsUs 


u 

s 


aA  - a + B„u  aA  - a B 


s s s 


(at  ) 


(A8) 


with  Equation  A5  showing  a linear  relationship  between  wind  speeds  at 
20  ft  and  at  another  height . 

8.  Examination  of  the  distribution  graphs  (Figure  A1 ) indicates 
that,  in  the  range  of  velocities  of  interest,  the  shifts  have  in  general 
been  parallel.  If  major  population  changes  had  occurred,  the  lines 
should  show  greatly  differing  tilts.  It  is  also  evident  from  Figure  A1 
that  some  fluctuations  do  occur  for  very  high  wind  speeds  in  part  due 

to  sampling.  Estimates  of  the  relationship  properly  rely  on  inter- 
mediate values  of  wind  speed.  Each  seasonal  distribution  function  is 
based  on  approximately  5,000  to  25,000  observations.  Hence,  the  assump- 
tion that  population  characteristics  of  winds  over  land  are  statisti- 
cally stable  appears  to  be  reasonable. 

9.  For  changes  in  instrument  height  or  small  spatial  movements 


Ah 


of  site,  the  values  predicted  from  Equation  A 5 agree  well  with  the 
one-seventh  power  law.  If  the  spatial  change  is  several  miles,  however, 
the  results  differ  greatly. 

10.  In  summary,  transformation  following  the  method  of  Equa- 
tions Ah  and  A8  can  be  derived  to  adjust  wind  speeds.  The  primary  as- 
sumption that  must  be  met,  however,  is  that  the  wind  profile  is  neutral 
on  the  time  scale  of  average  hourly  winds.  If  the  site  is  a land  site 
away  from  immediate  lake  influence,  this  assumption  is  usually  assured 
and  is  also  probably  good  for  post-19^8  wind  stations,  as  shown  in  the 
next  section. 


Pre-19^8  Winds 

11.  Prior  to  19^8,  most  available  weather  data  were  taken  at 
sites  located  in  the  downtown  areas  of  the  cities  and  not  at  the  air- 
ports. Among  the  problems  that  had  to  be  resolved  before  transformation 
of  these  winds  were  (a)  the  proximity  of  these  stations  to  the  lake, 

(b)  the  shifting  of  station  heights  and  locations,  and  (c)  the  unavail- 
ability of  the  data  on  magnetic  tape  or  computer  cards. 

12.  The  pre-19^8  data  are  only  available  in  autograph  form. 

Hence,  the  cost  of  applying  the  distribution  function  arguments  of  the 
previous  section  was  prohibitive.  However,  the  shifts  in  station  loca- 
tion were  small  spatially  and  the  variations  in  height  were  large;  thus, 
the  one-seventh  power  law  could  be  reasonably  applied.  The  effect  of 
large  buildings,  typical  of  downtown  areas,  upon  the  wind  climate  must 
be  regarded  as  unknown  at  this  time.  However,  the  extreme  height  of  the 
site  locations  (100-300  ft)  atop  these  buildings  somewhat  alleviates  the 
problems  associated  with  funneling  between  buildings. 

13.  A major  problem  that  must  be  confronted  for  these  stations 
involves  their  proximity  to  the  lake.  For  winds  blowing  offshore,  the 
assumption  of  a neutral  land  profile  is  probably  justified;  the  winds 
can  be  transformed  to  20  ft  using  Equation  A8.  After  this  adjustment 
in  height,  the  land  and  lake  transformation  formulae  of  Part  II  are 
used  to  estimate  the  lake  winds. 
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lU.  For  winds  blowing  onshore,  the  boundary  layer  associated  with 
land  roughness  elements  grows  in  height  with  distance  inland.  Hence,  if 
the  instrument  site  is  sufficiently  high  and  near  the  lake,  it  is  in  the 
marine  boundary  layer;  therefore,  stability  effects  must  be  considered. 

15.  If  flow  from  the  lake-wind  regime  onto  land  is  regarded  as  an 
analogue  to  the  classical  problem  of  incompressible,  turbulent  flow  over 
the  leading  edge  of  a flat  plate,  the  depth  of  the  land  boundary  layer 
(6#)  can  be  expressed  as  (Shames'^) 

6*  = (v~)  /5*J/5(0.376)  (A9) 

where  ° 

v = dynamic  viscosity  of  air 

Vq  = wind  velocityt 

= distance  inland  from  the  lake-land  interface 
When  this  analogue  is  used  for  a velocity  of  50  mph  and  temperatures 
typical  of  the  Great  Lakes  environment,  the  thickness  of  the  land 
boundary  layer  is  calculated  to  be  25  ft  at  1 mile  from  shore,  h2  ft 
at  2 miles,  58  ft  at  3 miles,  and  88  ft  at  5 miles. 

16.  Many  stations  prior  to  19^8  were  located  within  2 miles  of 
the  lake.  Thus,  although  the  calculated  boundary-layer  thickness  is 
only  a crude  estimate,  it  appears  that  such  stations  were  in  the  marine 
boundary  layer.  Likewise,  all  of  the  post-19^8  land  sites  were  located 
more  than  5 miles  inland  and  have  anemometers  at  heights  well  within  the 
land  boundary  layer  as  predicted  by  Equation  A 9. 

17.  For  cases  in  which  stations  appeared  to  fall  within  the  ma- 
rine boundary  layer,  a theoretical  form  of  the  wind  profile  was  applied 
to  transform  winds  from  one  level  to  another.  In  this  transformtion, 
the  similarity  form  of  a wind  profile  in  the  presence  of  vertical  tem- 
perature gradients  is  given  by 

“(Z)  = p-  In  (A10) 

t In  the  classical  problem,  the  velocity  V is  constant  with  height  up 
to  the  edge  of  the  plate;  however,  in  this  case,  V must  be  replaced 
by  a velocity  averaged  over  height. 
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where 


K = Von  Karman  constant 
Zq  = roughness  height 

<f>  = modification  function  for  wind  profile  due  to  thermal  effects 
L'  = stability  length 

For  stable  conditions,  <p  has  been  found  to  be  representable  as  a 
linear  function  of  Z/L' 


<t> 


BZ 

L' 


(All) 


where  B , the  linear  constant  for  log-linear  wind  profile,  equals  -7  . 
In  unstable  conditions,  (J>  is  given  by 

*■/(■* + T*7  * 777  * 77?- *)"  (Al2) 

where  c is  nondimens ional  wind  shear  and  6 is  the  integration  vari- 
able. The  relationship  between  e and  Z/L'  is  specified  by  the  KEYPS 

52 

function  constant  k (Panofsky  ) 


k 

c 


kZ 

L' 


3 


e 


1 = 0 


(A13) 


with  k = 18  . 

18.  These  equations  can  be  solved  to  yield  u-  , Z , and  L' 

* o 

as  a function  of  wind  speed,  anemometer  height,  air-sea  temperature  dif- 
ference, and  height  of  air  thermometer  above  water  surface.  A numerical 

19 

solution  and  computer  routine  developed  by  Cardone  are  incorporated 

into  the  wave-hindcast  model.  Once  u#  , Zq  , and  L'  are  computed, 
the  wind  speed  at  another  level  can  be  obtained  by  the  use  of  Equa- 
tion A10;  hence,  wind  velocities  at  the  desired  anemometer  elevation 
(19.5  m)  can  be  calculated. 

Transformation  to  Winds  Over  the  Lake 

19-  After  winds  have  been  transformed  to  a standard  elevation 
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of  20  ft*  (by  the  one-seventh  power  law  for  pre-19^8  winds  and  by  the 
regression  relationships  for  post-19^7  winds),  the  transformation  to 
winds  over  the  lake  is  given  by 

= RU^  (All*) 

where  is  the  wind  velocity  over  water,  R is  a function  of  wind 

velocity  and  air-sea  temperature  difference  (T  - T ),  and  U0  is  the 

3.  S JC 

wind  velocity  over  land 

R = R(U.,T  - T ) (A15) 

X/  a,  S 

20.  Figure  A2  gives  the  relationship  between  R and  . The 
dependence  of  R on  air-sea  temperature  difference  is  given  in  Fig- 
ure 11  of  the  main  text. 


Figure  A2.  Velocity  ratios 
as  a function  of  wind  speed 
over  land 


* This  does  not  refer  to  winds  which  lie  in  the  marine  boundary  layer, 
as  given  by  Equation  A9. 
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APPENDIX  B:  TREATMENT  OF  GENERAL  ASPECTS  OF 
NUMERICAL  WAVE  MODEL 


1.  The  spectral  energy  density  in  a wave  spectrum  is  treated  as 
a five-dimensional  variable.  It  is  a function  of  surface  location 
x(x,y)  , wave  number  vector  kCk^kg)  , and  time  (t),  i.e.,  F(x,k,t)  . 
The  wave  number  vector  can  alternately  be  viewed  in  terms  of  frequency 
(f)  and  propagation  direction  (0).  The  total  energy1*  at  point  (xo,yo,tQ) 
can  be  obtained  by  integrating  it  over  wave  number  space 


2 

o 


11 

kl  k2 


F(xo’yo,to’kl’k2)  <*1  ^2 


(Bl) 


or  in  frequency  and  direction  space 


2 

a 


/ fF(x o’yo’Vf’0)  df  d6 
f e 


(B2) 


2.  Ocean  waves  are  a statistical  phenomenon.  The  complex  manner 
in  which  they  interact  with  other  media  at  boundaries  along  internal 
boundaries,  or  even  within  a spectrum,  can  be  described  only  in  terms  of 
statistical  averages  over  space  and  time.  When  k and  x are  given  as 
the  wave  number  vector  and  location  vector,  respectively,  the  transfer 

of  energy  into  and  out  of  a wave  component  can  be  written  in  the  form  of 

53** 

an  energy  balance  or  radiative  transfer  equation  (Gelci  and  Cazal£, 

5I+. 

Hasselmann  ) 


DF  _ 3F  . 3F_  . 3F 

Dt  3t  Xi  3x.  i 3k. 

l l 


S 


where 


3F/3t  = local  (in  space)  rate  of  change  of  F 
x^(3F/3x^)  = effect  of  advection  of  energy  into  F 


(B3) 


* In  actual  application,  the  term  F does  not  have  the  units  of 

energy  but  rather  of  length  squared.  It  also  represents  that  portion 
of  the  variance  of  the  surface  elevation  in  a particular  frequency 
band  and  direction  band  at  location  x . 

**  Raised  numbers  refer  to  similarly  numbered  items  in  "References"  on 
pp  1+ 5-1+9  at  end  of  main  text. 
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k^(9F/3k^)  = effect  of  changes  in  the  wave  number  on  F 

S = any  process  tlja^  transfers  energy  into  and  out  of  wave 
component  F(x,k,t) 

The  dot  over  and  k^  denotes  differentiation  with  respect  to  time, 

and  the  subscript  i represents  the  conventional  summation  over  all 
dimensions.  In  a Lagrangian  context,  DF/Dt  represents  the  rate  of 
change  of  F as  it  propagates  through  x to  k space  along  the  path 
determined  by  the  Hamiltonian  equations. 

xi  = a(x,k)  (BU) 

i 

ki  = cr(x,k)  (B5) 

i 


where  a is  the  circular  frequency  associated  with  F . In  deep  water 
or  in  water  of  constant  depth  d , the  refraction  term  (k^3F/3k^)  van- 
ishes. In  this  case.  Equation  B4  becomes 


x. 

1 


constant 


(B6) 


where  the  constant  is  the  group  velocity  of  waves  with  circular  fre- 
quency in  water  of  depth  d . 
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3.  Hasselmann  demonstrated  that  the  source  term  S can  be 
represented  in  a general  form  by 


S(k)  = S±  + s2  + s3  + + SJ  + S£  + + S£  + s5  + s6  + s7  (BT) 


where 


= a 

52  = BF(k)  in  which  B is  the  coefficient  of  coupling  be- 

tween surface  waves  and  mean  wind  profile 

53  = F(k)/  y(k,k')F(k')  dk' 

= -SF(k)  + /e(k,k')F(k')  dk' 

S|...S£  = same  functional  form  as  S^.-.S^ 


S5  = SSfl  T1F(k,)F(k")F(k  - k'  - k")  - T2F(k)F(k' )F(k" ) dk’  dk" 

in  which  and  T2  are  transfer  functions  and  are  depen- 

dent on  the  coupling  coefficients  between  spectral  components 
involved  in  the  interactions 

= -v. .k.k.F(k)  with  v..  as  the  viscosity  sensor 
6 ij  1 j lj 

= an  undetermined  dissipation  function  due  to  wave  breaking 
The  processes  embodied  in  these  terms  are  as  follows: 

a.  S (a)  - energy  transfer  to  the  wave  field  through  turbu- 
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lent  pressure  fluctuations  in  the  atmosphere  (Eckart, 

Phillips3*1). 
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b.  S2  - energy  transfer  due  to  Miles  instability  mechanism. 

c_.  - energy  transfer  in  a nonlinear  correction  to  Miles' 

theory. 

d.  S,  - energy  transfer  due  to  wave-turbulence  interactions 

CfZ. 

(Hasselmann 5 ). 

e_.  Sj  , , ^3  > - energy  transfers  due  to  analogous 

processes  to  , and  , except  the  inter- 

action is  between  waves  and  other  ocean  phenomena  rather 
than  waves  and  atmospheric  phenomena. 

_f . S,.  - energy  transfer  due  to  nonlinear  wave-wave  interac- 

tions (Hasselmann36,37,38,56). 

g_.  Ss  - loss  of  energy  from  waves  due  to  turbulent  bottom 
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friction  (Hasselmann  and  Collins^). 

U.  When  the  processes  are  presented  in  this  way,  the  multitude  of 
theoretical  treatments  of  the  wind  generation  of  waves  can  be  put  in 
some  context  and  thus  compared.  However,  the  theory  presented  by 
Hasselmann  covers  only  weak  interactions.  Any  processes  that  violate 
the  assumption  of  a slowly  varying  spectrum  and  some  nonexpansible  pro- 
cesses are  not  included  in  Equation  BJ.  It  is  still  a definite  im- 
provement over  the  unorganized  development  of  empirical  and  theoretical 
mechanisms  for  wave  growth  and  provides  a means  of  comparing  and  evalu- 
ating source  functions  proposed  by  different  investigations. 

5.  Source  functions  for  the  two  models  used  in  the  study  can  be 
written  as  follows: 


B3 
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a.  For  the  Inoue  model  using  modified  Miles-Phillips 
mechanism,! 


S,(l  - q ) + S*(l  - ) 


(B8) 


where 

sum  of  sources 

ratio  of  energy  density  in  one-dimensional 
spectrum  E(f)  to  that  in  fully  developed 
spectrum  E (f) 

CO 

As  given  by  Pierson  and  Moskowitz, 


where 

a = Phillips  equilibrium  constant  (0.008l) 
g = gravity 

8 = energy  transfer  due  to  interaction  between  mean 
wind  profile  and  waves  (as  a constant  = O.7I+) 

= wind  frequency  (ratio  of  g to  wind  speed  u ) 
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b_.  For  the  Barnett  model  using  parameterized  wave-wave  in- 
teractions. 


SB  ’ (S1  + S2)(l  " *)  + s5  (BIO) 

where  Sg  represents  the  sum  of  sources.  The  limiting 
function  <f>  for  wind  input  equals  d exp  [-d2(R*-  F)/F] 
with  dg  and  d^  constants  and  R , the  angular  spread- 
ing function,  defined  by  Phillips  equilibrium  spectrum  as 

2 

R(f,0)  = h( 0 ) (Bll) 


t The  actual  form  of  Sg  in  the  Inoue  model  was  a one-dimensional 
spectrum  in  place  of  a two-dimensional  spectrum.  S*  is  given  by 
BE(f)'  as  opposed  to  BF(f,9)  where  F(f,0)  is  the  two-dimensional 
spectrum. 


6.  For  numerical  processing,  Barnett ^ chose  a normalized  cos ^ 0 
distribution  of  energy  as  a function  of  direction  for  h(0)  , the  spread- 
ing factor;  that  is, 


h(6)  = I1  cos1*  0 (B12) 

7.  With  two  such  different  source  functions,  one  might  expect 
vastly  different  results  from  the  application  of  these  two  models;  how- 
ever, this  is  not  the  case.  The  inability  of  experimental  procedure  to 
adequately  separate  the  different  source  terms  has  led  to  a situation 
in  which  source  terms  in  hindcast  models  are  most  likely  an  empirical 

mixture  of  more  than  one  theoretical  source  term.  As  an  example  of  this, 

59 

the  same  data  from  the  limited-fetch  study  by  Snyder  and  Cox  have  been 
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used  by  both  Barnett  and  Inoue  to  justify  totally  different  func- 
tional forms  for  source  terms.  Yet,  in  most  cases  the  two  models  appear 
to  coincide  rather  well,  at  least  in  terms  of  gross-scale  parameters 
such  as  significant  wave  height  and  frequency  of  the  maximum  spectral 
energy  density.  It  is  easier  to  compare  the  functional  forms  of  the 
two  models  if  the  limiting  terms  are  temporarily  discarded.  The  full 
radiative  transfer  Equation  B3  is  solved  by  Barnett,  whereas  Inoue 
reduces  the  equation  to  its  linear  form 

at  ■ 5l  + s*  (B13) 

In  the  linear  form  of  this  model,  energy  is  "jumped"  from  one  grid  point 
to  another  in  a somewhat  erratic  manner. 

8.  During  the  preparation  of  some  preliminary  results  for  sites 
on  Lake  Erie  and  Lake  Ontario,  it  was  determined  that  the  one- 
dimensional nature  of  the  Inoue  source  function  could  not  adequately 
treat  variations  in  shoreline  geometry  typical  of  the  Great  Lakes.  The 
asymmetry  of  the  spectrum  with  respect  to  propagation  direction  could 
not  be  accounted  for  near  the  coast.  Winds  blowing  parallel  to  a coast 
thus  tend  to  generate  large  waves  in  proximity  to  the  coastline.  Other 
areas  of  weakness  were  the  difficulty  of  treating  a rapidly  changing 
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wind  direction  and  the  lack  of  consideration  given  to  the  variation  of 
fetch  in  cases  with  highly  irregular  shorelines.  For  most  wave  genera- 
tion sequences  on  the  Great  Lakes,  the  one-dimensional  treatment  still 
gave  very  reasonable  results;  however,  under  some  conditions  wave 
heights  appeared  to  grow  improbably  fast  and  at  other  times  improbably 
slow. 

9-  The  success  of  the  Inoue  model  in  deep  water  (where  the 
spatial  gradients  are  usually  somewhat  small)  makes  it  very  difficult 
to  dismiss  the  validity  of  the  source  function,  at  least  in  terms  of 
reproducing  significant  wave  heights  and  periods.  This  prompted  a mod- 
ification of  the  Inoue  source  function  to  the  form 

ST  = a + BF<f>(e,0  ) (BlU) 

1 w 

where  <}>  ( 0 , 0 ) represents  an  attempt  to  account  for  the  effect  of  the 
w 

angle  between  the  propagation  direction  0 and  wind  direction  0 . A 

w 

simple  form  of  $(0,0^)  equal  to  cos  |@  - 0^|  was  found  to  perform 

reasonably  well  in  generating  a spectrum  with  a spreading  factor  pro- 
2 

portional  to  cos  0 around  the  mean  propagation  direction.  Recent 

field  studies  (Hasselmann  et  al.^)  have  shown  that  such  a spreading 

factor  gives  the  best  fit  to  observations.  Deviations  from  this  form 

and  detailed  fine  structure  in  the  directional  distribution  of  wave 

energy  have  been  found  to  be  irreproducible  (Tyler  et  al.^). 

10.  By  using  Equation  B13,  the  Inoue  model  can  be  compared  with 
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the  Barnett  model.  Barnett  parameterized  the  net  transfer  of  energy 
due  to  wave-wave  interactions  and  expressed  the  general  form  of  this 
nonlinear  (n£)  energy  transfer  as 

(f)  nl  = r - <B15) 

where  T represents  the  gain  of  energy  of  wave  component  F from  all 
interactions  in  which  it  plays  a passive  role  and  -tF  represents  the 
loss  of  energy  due  to  all  the  interactions  in  which  it  plays  an  active 
role.  For  a slowly  varying  spectrum  (as  necessary  for  the  theory  of 


B6 


r 


weak  interactions  to  be  valid),  the  form  of  the  radiative  transfer 
equation  at  time  tQ  for  the  two  different  source  functions  is  (again 
neglecting  the  limiting  terms  in  each  model)  written  as 
Inoue  Model 

If  = Oj  t 6jP  - 00  • VF  (Blfi) 

and 

Barnett  Model 

||  = + T + (6B  - t)F  - CQ  • VF  (BIT) 

where  C,,  represents  group  velocity  and  the  subscripts  I and  B re- 

fer  to  the  specific,  empirically  determined  source  terms  used  by  Inoue 

and  Barnett,  respectively.  The  superficial  similarity  is  striking.  It 

would  seem  that,  if  a_  = a + T and  B_  = B_  - t , the  two  source  func- 

Id  Id 

tions  would  be  equivalent.  However,  whereas  aT  , a_  , BT  , and  3_ 

1 D 1 D 

are  functions  only  of  wind  speed,  gravity,  wave  direction,  and  wave  fre- 


quency,  the  functional  forms  of  r and  F , wave-wave  interaction 

pa- 

rameters,  are 

t = f (g,E,f ,fQ ,0 , 0q) 

(Bl8) 

and 

r = r(g,E,f,fo,e,eo) 

(B19) 

Thus  the  following 

three  new  parameters  E , f 

, and  0 are 

o 

introduced: 

E = j j F(f,e)  df  de 

(B20) 

fo\[t  JJr(f,e)r2  af  de] 

1/2  \ 

(B21) 

and 

J /F(f’0)  Sln  0 df 

0 = tan-1  — 

/ j F(f,e)  cos  e df 

\ # 

de  / 

(B22) 

de  ) 

* These  terms  are  given  as 

fQ  = | //  F(f,0)f  df  d0  and  0O  = E // F(f»0)0  df  de 
. . 30 

in  the  original  work  by  Barnett  but  cannot  be  correctly  computed 
from  these  forms. 


BT 


11.  These  new  parameters  couple  the  source  function  to  the  local 
spectrum,  something  not  done  in  the  Inoue  model.  This  change  is  very 
important  at  least  in  terms  of  evaluating  the  physical  interpretation 

of  the  energy  transfers.  In  the  Barnett  model,  the  source  function  pro- 
duces an  eventual  quasi-equilibrium,  whereas  the  source  function  must  be 
artificially  cut  off  once  the  spectral  energy  has  reached  a velocity- 
dependent  form.  Since  the  wave-wave  model  does  not  require  such  an  un- 
realistic cutoff,  it  should  be  interpreted  as  being  a better  represen- 
tation of  the  actual  physics  of  wave  growth. 

12.  The  dependence  of  the  wave  spectra  on  wave-wave  interactions 
has  received  additional  support  in  recent  years  by  the  documentation  of 
an  "overshoot”  effect  (Barnett,  ^ Hasselmann  et  al.^).  Figure  B1  shows 


a typical  time  sequence  of  wave  spectra  from  Toronto.  The  presence  of 
energy  in  the  spectral  peak  above  that  predicted  on  a strict  similarity 
basis  (~f  is  evidence  that  the  spectrum  is  also  affected  by  a self- 


similarity criterion  , f/f 

|_  V max , 

tween  wave  breaking  and  spectral  peak. 


suggesting  a relationship  be- 


13.  The  parameterizations  used  by  Barnett  to  represent  the  effect 


of  wave-wave  interactions  are 


7 2 7 
7.5  x lO'E  f1 


1 + 16  cos  (9  - 0 ) I (f  - 0.53  f )' 
o o 


(B23) 


for  f > 0.53  fQ  ; otherwise,  t = 0 


k.k  x l08E3f8 
o 


4 / 

cos  (9  - 


/ f \ 
( 1 - °-'*2  -) 

X exp  -h(l  - 


+ 0.1 


(B2k) 


for  f > 0 . i+2  f and  18  — 0 I < tt/2  ; in  all  other  conditions  F = 0 . 
o 1 o 

Comparisons  between  calculations  of  energy  transfer  rates  based  on  these 


parameterizations  have  shown  good  agreement  for  an  empirical  Neumann 
spectrum  (Barnett  ).  Barnett  contends  that  differences  in  spectral 
shape  do  not  alter  the  wave-wave  transfer  rates  excessively  since  they 
are  dependent  on  integral  expressions  of  spectral  energy.  However,  re- 
sults from  the  JONSWAP  (Joint  North  Sea  Wave  Project)  study  indicate 
that  spectral  shape  may  significantly  alter  rates  of  energy  transfer 
under  some  conditions.  Even  so,  the  time  necessary  to  evaluate  the 
theoretical  form  of  the  wave-wave  interactions  (30  min  on  a CDC  6600 
per  spectrum;  Sell  and  Hasselmann^ ) makes  their  practical  application 
to  hindcasting  infeasible  at  present. 


lk.  Up  to  this  point  only  source  terms  , 


S^  , and  have 


been  considered.  Terms  and  are  possibly  included  in  the  em- 
pirical forms  of  and  ; terms  Sj  , , and  are 

neglected  on  the  Great  Lakes.  The  interaction  between  the  spectrum 

and  the  bottom  (S^)  in  general  cannot  be  neglected.  Hasselmann  and 
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Collins  proposed  the  following  spectral  dissipation  function 


calculated  from  a quadratic  stress  law  using  a perturbation  technique  on 
the  equations  of  motion: 


friction 


4 - 


(B25) 


with 


v 


ij 


gc. 


2 2 
a cosh  kh 


6 . . <u>  + 
ij 


(B26) 


where 

= coefficient  of  bottom  friction 

6 . . = Kronecker  delta 
ij 

u = fluid  velocity  relative  to  the  bottom 

The  brackets  <>  indicate  averages  taken  over  the  entire  ensemble  of 

Uq 

al £.  possible  wave  fields  having  the  specified  spectrum.  Collins  uses 
a decoupled  form  of  Equation  B25  for  numerical  purposes.  In  this  ap- 
proximation, the  frictional  dissipation  function  is  defined  as 


— C gkC 

*(f,e)  = — 2^ — 2 — * F(f’0)  * <u>  (B2T) 

2t to  cosh  kd 


with 


<u> 


o cosh  kd 


(B28) 


2 2 2 2 

The  multiplication  of  E(f)  by  g k /(o  cosh  kd)  in  Equation  B28 

converts  the  surface  spectrum  to  a bottom  velocity  spectrum. 

15.  The  only  term  now  unaccounted  for  is  S7  , energy  loss  due 
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to  breaking.  Phillips  hypothesized  on  a similarity  basis  that  the 
equilibrium  form  of  the  spectrum  should  be  independent  of  wind  speed. 
Although  the  rate  of  energy  transfer  through  a fully  developed  spectral 
component  would  increase  with  an  increase  in  wind  speed,  the  net  gain 
would  still  be  zero.  The  Inoue  and  Barnett  models  both  have  limiting 

factors  in  their  original  presentation.  For  the  Inoue  model,  it  is 

(2  2 \ 1/2 

1 - Eq/EooJ  ; and  for  the  Barnett  model,  it  is  given  by 


r' 


1 - <j)  (at  least  for  the  wind  input).  In  light  of  recent  evidence  for 
the  tendency  of  an  "overshoot"  effect,  the  need  for  such  limiting  terms 
as  a spectral  component  approaches  the  spectral  energy  is  dubious. 

Thus , the  source  terms  are  given  instead  by  Equations  Bl6  and  B17  in 
which  there  are  no  limiting  terms.  This  facilitates  the  expansion  of 
Equations  Bl6  and  B17  and  the  numerical  evaluation  by  a Taylor  series 
resulting  in 


(S*  + Sn*F)  {Sl  + SnJtF)  2 

F + AF  = F + — - — At  + — ■ ol  - — At  S . 

1!  d\  nx. 


, (St*S»«F>  .,3.2  , , <St*S„tF>  . 

3!  nil  n!  nil 


(B29) 


for  F + AF  < E ; otherwise,  F + AF  = E . Terms  S„  and  S „ rep- 
resent  the  sum  of  linear  source  terms  and  the  sum  of  nonlinear  source 
terms,  respectively.  It  should  be  noted  that  both  one-dimensional  and 
two-dimensional  forms  are  involved  in  this  calculation.  Details  of 
these  numerical  methods  will  be  given  in  the  sixth  report  of  this  series. 

1 6.  Experimental  research  is  continuing  in  attempts  to  improve 

the  source  terms.  However,  the  forms  of  the  source  function  used  in 
this  report  have  performed  well  in  previous  comparisons  between  hind- 
casts  and  observed  waves  (Part  III).  Such  success  supports  the  conten- 
tion that,  if  other  processes  are  imported  besides  , s2  ’ an<*  S5  ’ 

they  are  effectively  included  in  the  empirical  evaluations  of  the  co- 
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efficients  of  , and  . Thus,  even  though  Hasselmann 

points  out  that  term  should  be  larger  than  S 2 , the  actual  form 

of  B used  in  hindcast  models  is  probably  the  sum  of  effects  from  both 
and  . Whether  numerical  models  treat  each  source  term  in  Equa- 

tion B7  separately  or,  in  fact,  compound  several  terms  is  beyond  the 
scope  of  this  study.  The  critical  question  is  the  overall  accuracy  of 
the  models  relative  to  the  needs  of  users  of  wave  information  on  the 
Great  Lakes.  In  the  end,  this  is  assessed  only  by  comparisons  with 
actual  data  and  the  needs  of  the  user. 

17.  In  Part  III,  some  results  from  the  wave  model  used  in  this 
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study  were  compared  with  wave  gage  records  at  Point  Pelee  on  Lake  Erie 
and  Toronto,  Cobourg,  and  Main  Duck  on  Lake  Ontario.  These  comparisons 
showed  the  relationship  between  maximum  observed  significant  wave 
heights,  a single  time  series,  and  a single  spectrum.  Additional  com- 
parisons are  given  in  Figures  B2-B10. 


DAYS.  NOVEMBER  1973 


DAYS,  SEPTEMBER  1973 

Figure  B2.  Comparisons  of  hindcast  and  ob- 
served sequences  of  significant  wave  height 
at  Point  Pelee,  with  the  best  agreement  shown 
at  the  top  and  the  worst  at  the  bottom 


2 13  V 

DAYS,  NOVEMBER  1972 
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Figure  B3.  Comparisons  of  hindcast  and  ob- 
served sequences  of  significant  wave  height 
at  Toronto,  with  the  best  agreement  shown  at 
the  top  and  the  worst  at  the  bottom 
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DAYS.  APRIL  1972 


DAYS,  NOVEMBER  1972 

Figure  Bh.  Comparisons  of  hindeast  and  observed 
sequences  of  significant  wave  height  at  Cobourg, 
with  the  best  agreement  shown  at  the  top  and  the 
worst  at  the  bottom 
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Figure  B5.  Comparisons  of  hindcast  and  ob- 
served sequences  of  significant  wave  height 
at  Main  Duck,  with  the  best  agreement  shown 
at  the  top  and  the  worst  at  the  bottom 
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rave  periods  at  Point  Breeze  with  the  number  of  wave  periods  at  Oswego  with  the  number  of  occur- 
iccurrences  of  significant  periods  at  Cobourg,  rences  of  significant  periods  at  Main  Duck,  both 
both  as  a function  of  significant  wave  height  as  a function  of  significant  wave  height 


18.  In  Figures  B2-B5,  hydrographs  of  significant  wave  heights 
for  selected  periods  of  time  are  shown  for  each  of  the  four  sites 
previously  specified.  For  each  site,  the  hydrograph  with  best  agreement 
and  that  with  worst  agreement  between  observed  and  hindcast  waves  are 
given.  In  general,  the  phasing  and  magnitude  of  the  two  series  agree 
rather  well,  with  larger  wave  heights  tending  to  show  a somewhat  higher 
correlation. 

19.  In  Figure  B6,  observed  frequencies  of  spectral  maxima  are 
compared  with  hindcast  frequencies  of  spectral  maxima  for  the  same  peak 
storm  conditions  included  in  Figure  18.  This  graph  indicates  a tendency 
to  underestimate  f at  high  frequencies  but  shows  good  agreement  at 
frequencies  less  than  0.15  Hz.  This  effect  is  due  in  part  to  some  of 
the  numerical  assumptions  included  in  the  wave  model  in  order  to  reduce 
computation  time  at  the  high-frequency  end  of  the  spectrum. 

20.  For  all  grid  points  shown  in  Figure  21,  the  relationship  be- 
tween significant  period  and  significant  wave  height  was  obtained  from 
model  outputs.  Since  both  significant  period  and  significant  height  are 
functions  of  nondimens ional  fetch,  they  must  also  be  highly  related  to 
each  other  in  situations  where  little  swell  is  present.  Consequently, 
sites  in  similar  geometric  situations  should  have  similar  relationships 
between  significant  height  and  significant  period.  To  test  this,  re- 
gression lines  were  constructed  for  significant  period  as  a function 

of  significant  height  for  the  three  angle  classes  of  wave  approach 
(Appendix  E).  These  lines  are  drawn  over  tables  of  peak  wave  period 
versus  significant  wave  height  taken  from  the  Canadian  gage  data  (Fig- 
ures B7-B9).  In  each  case,  the  sites  compared  are  across  the  lake  from 
each  other  and  should  have  similar  fetch  geometries.  The  numbers  in 
Figures  B7-B9  represent  the  number  of  occurrences  of  waves  in  each  1-ft 
by  1-sec  category.  These  figures  indicate  uhat  significant  periods 
hindcast  by  the  numerical  model  are  in  close  agreement  with  periods 
bserved  on  the  Great  Lakes. 

21.  In  addition  to  the  comparisons  between  gage  data  and  hindcast 

2 

• , the  rate  of  growth  of  nondimensional  wave  height  (H  = g/E/u-) 

2 

. • dimensional  fetch  (F  = gF/u#)  computed  with  the  wave  model  was 
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compared  with  field  evidence  from  Mitsuyasu  and  Hasselmann  et  al. 

As  seen  in  Figure  BIO,  the  agreement  between  the  wave  model  and  these 
data  is  excellent. 


APPENDIX  C:  WIND  INTERPOLATION  OVER  THE  LAKE 


1.  When  wind  speed  and  direction  are  given  at  a series  of  stations 
either  along  the  lake  periphery  or  on  the  lake  itself,  it  is  necessary 

to  interpolate  corresponding  wind  speed  and  direction  data  at  every  grid 
point  used  in  the  calculation  of  the  wave  spectra.  The  method  of  inter- 
polation  generally  follows  that  of  Platzman. 

2.  For  a grid  point,  the  distance  to  each  location  for  which  input 
data  are  available  is  calculated  and  is  denoted  by  d.  , i enumerating 
the  data  input  points  and  ranging  from  1 to  N.  The  following  N products 
are  formed: 


\ = Hj^dj  , (i  = 1,  N),  (j  = 1,  N) 
where  as  an  example, 

*1  = d2  * V--dN 


A normalizing  factor  is  calculated 

N 

NORM  = V JL 
i=l 

and  a series  of  N weights  is  formed 

l. 

W = — — 

i NORM 


where  NORM  is  the  sum  of  the  multiple  products  £.  and  W.  repre- 

* ^h1 

sents  the  weighting  of  the  interpolation  function  on  the  l wind 
station. 

3.  To  interpolate  wind  speed  at  the  grid  point  the  following 
weighted  sum  is  used: 


U = W.  U,  + W0U0  + . . . + W U 
11  2 2 n n 


where  the  IP  is  input  wind  data  at  the  i 


.th 


station . 


* Raised  numbers  refer  to  similarly  numbered  items  in  "References"  on 
pp  I5-I9  at  end  of  main  text. 
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U.  The  vectorial  wind  velocity  was  decomposed  into  U (east-wesc) 
and  V (north-south)  components.  Each  component  was  interpolated  over 
the  lake  and  recombined  to  give  a wind  speed  and  direction.  As  a further 
modification,  if  adjacent  input  data  stations  had  wind  directions  within 
90  deg  of  each  other,  directions  were  produced  using  the  Platzman  scheme, 
but  wind  speeds  were  assigned  the  average.  This  modification  is  to  allow 
for  the  rotation  of  the  winds  without  divergence  and  a lowering  of  the 
wind  speed. 


APPENDIX  D:  ANALYSIS  OF  EXTREMES 


1.  Since  curvatures  in  the  upper  portions  of  the  graphs  of  wave 
height  H versus  return  period  T*  were  small,  only  the  problem  of 
estimation  of  parameters  for  the  Fisher-Tippett  Type  I distribution  need 
be  considered  here.  For  this  distribution,  the  asymptotic  distribution 
function  is  given  by 


F(H)  = exp  -exp  £-a(H  - u)J 


(Dl) 


where  a and  u are  parameters**  of  the  distribution  and  must  be  es- 
timated from  the  sample  characteristics.  From  Equation  Dl,  the  reduced 
variate  y can  be  defined  as 


y = o(H  - u) 


(D2) 


U2+ 


and,  as  shown  by  Gumbel,  the  reduced  variate  and  the  return  period 
are  related  by 


y = -In  Jin 


T - 1 


(D3) 


For  increasing  values  of  T , Equation  D3  converges  to 


y = Jin  T - -g  T 


(DU) 


which  holds  with  an  error  of  less  than  0.7  percent  from  T = 7 onward. 
Solving  for  T , one  obtains 


T = ey  + | (D5) 


where,  for  large  T , combining  Equations  D2  and  D5  results  in 


* Calculated  by  the  relationship,  T = l/l  - F(H)  . 

**  The  units  of  a and  u are  l/length  and  length,  respectively, 
t Raised  numbers  refer  to  similarly  numbered  items  in  "References"  on 
pp  U5-U9  at  end  of  main  text. 
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T = ea(H-u) 


where  e represents  the  Napierian  base  (2.71828...)  or 

H = u + —5)  (D7) 

Equation  D7  clearly  demonstrates  the  expected  increase  in  extreme  wave 
heights  as  the  logarithm  of  the  return  periods.  Empirical  results  sup- 
porting this  conclusion  can  be  found  in  many  early  studies  of  flooding 
on  rivers.  It  is  this  relationship  that  allows  plots  of  large  values  of 
a variate  (wave  heights,  water  levels,  river  discharges,  etc.)  against 
return  periods  to  be  fit  by  a straight  line  on  semilogarithmic  paper. 
Hence,  the  parameters  u and  a can  be  calculated  by  fitting  a 
straight  line  to  the  upper  region  of  the  semilogarithmic  plots  of  wave 
height  versus  return  period.  The  slope  of  these  lines  is  an  estimate  of 
l/a  , and  the  extrapolated  value  of  H , for  T equal  to  1 , can  be 
used  as  an  estimate  of  u . This  method  is  somewhat  equivalent  to  the 
method  of  moments  for  estimating  a and  u . The  method  of  moments  is 
simple,  straightforward,  and  reasonably  efficient  (Gringorten^’^'’ ) ; and 
although  several  other  methods  for  calculating  a and  u are  available, 
the  method  of  moments  was  used  in  this  report. 

2.  The  problem  of  constructing  control  bands  for  estimates  of 
return  periods  is  relatively  simple  in  the  method  of  moments.  For  a 
Fisher-Tippett  Type  I distribution,  the  expected  rms  error  of  an  esti- 
mated return  period  is  given  by 


S = S’ 
x 


/l.lOOOy^  + 1.1396y  + 1 
N 


where  S is  the  standard  deviation  of  the  annual  wave  height  maxima  and 
N is  the  number  of  years  in  the  sample.  It  is  the  term  Sx  that  is 
presented  along  with  all  return  intervals  given  in  Appendix  E.  In  the 
following  discussion,  it  is  important  to  bear  in  mind  that  predictions 
of  extremes  are  based  on  characteristics  of  the  distribution  of  a vari- 
ate, not  on  individual  values  of  the  variate.  Hence,  errors  that  result 
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in  hindcasts  sometimes  too  low  and  sometimes  too  high  may,  in  fact,  be 
tolerable  if  the  range  of  errors  is  sufficiently  small.  This  range  is 
determined  by  the  sum  of  all  the  components  of  error  inherent  in  hind- 
casting: (a)  misrepresentation  of  wind  fields,  (b)  inadequate  specifi- 

cation of  boundary  conditions,  (c)  simplification  of  assumptions  in 
numerical  calculations,  and  (d)  lack  of  properly  formulated  source  terms 
in  the  radiative  transfer  equation.  The  total  error  must  be  separated 
into  two  parts,  bias  and  random  error.  The  bias  is  assumed  to  be  negli- 
gible after  proper  calibration  of  the  wave  model.  What  remains,  pos- 
sibly nonnegligible,  is  the  random  error;  and  it  is  random  error  and  the 
relationship  between  this  error  and  errors  in  estimates  of  return  pe- 
riods that  are  addressed  in  this  appendix. 

3-  In  previous  investigations  of  properties  of  extremes,  the  vari- 
ate was  assumed  to  be  a known,  measured  quantity  and,  hence,  contained 
zero  random  error.  In  wave  hindcasts,  the  computed  probability  of  wave 
heights  depends  on  the  combined  effects  of  the  actual  wave  height  prob- 
ability in  nature  and  errors  present  in  wave  hindcasts.  This  can  be 
considered  as  the  probability  statement 

p(H ' ) = J p(H' |H)p(H)  dH  (D9) 

where  H'  represents  hindcast  wave  heights  and  H represents  actual 
wave  heights.  If  the  actual  wave  heights  follow  a double  exponential 
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distribution,  the  result  is  written  as  (Gumbel  ) 

-a(H-u) 

p(H)  = — (DIO) 

n 

where  n is  the  number  of  independent  observations  in  each  annual  or 
seasonal  set  of  observations. 

k.  If  it  is  assumed  that  there  is  negligible  bias  and  that  many 
independent  factors  create  random  errors  in  wave  heights  hindcast  by 
numerical  techniques,  a loose  application  of  the  central  limit  theorem 
may  be  invoked  to  support  a Gaussian  distribution  of  these  errors. 

Thus,  we  obtain 


D3 


(Dll) 


where  a is  the  standard  deviation  of  the  errors  in  hindcast  wave 
heights.  Substituting  Equations  D9  and  DIO  into  D8  yields 


p(H'  ) 


(D12) 


Rearranging  Equation  D12  shows  that  the  probability  of  H'  can  be 
written  as  a product  of  a double  exponential  probability  density  and 
an  integral  function  of  the  parameters  a and  a 


p( H ' ) 


'a(H~U^  f _1_  e~H*" /2-ZH 

n J 6 


(D13) 


where  the  reduced  wave  height  variate  H is  expressed  as 


H* 


(DlU) 


and  parameter  Z of  the  combined  exponential  and  Gaussian  distribution 
is  defined  as 


Z = ac 


(D15) 


Integrating  Equation  D12  yields 


-cx(H— u) 

p(H')  = — e 


(di6) 


where  4>(x)  is  the  probability  integral 

x 

? r 

4>(x)  = — 

/iT  J 

o 

Thus,  the  distribution  function  can  be  calculated  from  the  relationship 


2 

-x 

e dx 


(D17) 


F(H' ) = 1 - - e-0^H-U^(Z) 


dU 


(Dl8) 


where  H’(Z)  is  defined  as 


H'(Z)  = e 


Z 2 1J2 


I - 4>  ( — 

\/2; 


(D19) 


Figure  D1  shows  the  relationship  between  fCz)  and  Z . 

5.  At  this  point,  it  is  necessary  to  recall  the  definition  of 
the  term,  Z = <xa  , in  order  to  return  to  the  original  point  of  this 
discussion,  the  effects  of  random  error  in  hindcasts  on  estimates  of 
return  periods.  In  this  definition,  the  parameter  c represents  the 
standard  deviation  of  the  errors  in  hindcasts,  and  a represents  an 
inverse  function  of  the  variation  of  the  annual  maxima.  If  we  choose  a 
new  parameter  8 to  be  proportional  to  the  variation  in  annual  maxima. 
Equation  D15  becomes 


Z = 


(D20) 


In  this  form,  the  meaning  of  Z becomes  more  intuitive.  Its  value 
represents  the  ratio  of  the  coefficients  of  variation  of  the  Gaussian 
and  double  exponential  distributions  contributing  to  p(H)  . If  the 
asymptotic  variation  is  much  larger  than  the  random  error  present  in 
hindcasts,  the  effect  of  the  random  error  is  very  small.  On  the  other 
hand,  if  the  random  error  is  larger  than  the  variation  of  annual  maxima, 
the  distribution  p(H' ) is  significantly  different  from  p(H)  . 

6.  In  Part  III,  the  rms  error  in  the  hindcasts  was  given  as 
0.5  ft.  When  the  distributions  of  extremes  hindcast  by  the  wave  model 
were  examined,  values  of  8 were  found  to  lie  between  1 and  10.  If  the 
rms  error  is  not  dependent  on  fetch,  the  values  predict  up  to  a 50  per- 
cent factor  of  conservatism  for  small  waves.  For  large  variations  in 
annual  maxima,  the  effect  of  the  random  error  in  hindcasts  can  clearly 
be  neglected. 
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APPENDIX  E:  TABLES  OP  DESIGN  WAVE  INFORMATION 


1.  Tables  in  this  appendix  contain  calculated  estimates  of  return 
periods  from  techniques  discussed  in  Part  IV  and  Appendix  D.  Raw  data 
input  into  the  calculations  are  wave  height  maxima  from  individual  storms 
These  storms  were  selected  with  objective  criteria  as  discussed  in 

Part  IV;  and  a numerical  model,  based  on  the  wave-wave  interaction  parame 
terizations  by  Barnett,  was  used  to  calculate  waves  from  each  storm. 

2.  Figure  21  in  Part  IV  shows  the  locations  of  all  the  deepwater 
sites  for  which  tables  are  listed.  The  wave  information  is  separated 
into  two  sets:  (a)  Tables  El  and  E2  with  wave  height  information,  and 
(b)  Table  E3  with  wave  period  information. 

3.  The  format  for  Table  El  is  as  follows: 

a.  A separate  page  is  given  for  each  point  shown  in  Figure  21. 

b.  With  each  site,  5-,  10-,  20-,  50-,  and  100-yr  significant 
wave  heights  are  listed  for  each  season. 

c_.  Control  band  (one  standard  deviation)  estimates  for  each 
wave  height  are  included  inside  parentheses. 

ji.  The  three  angle  classes  are  defined  as  viewed  by  an  ob- 
server standing  on  shore  (Figure  El): 

(l)  Angle  Class  1 - Mean  wave  approach  angle  greater  than 
30  deg  to  right  of  normal  to  shore. 
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Figure  El.  Definition  sketch  of  the 
angle  classes 


El 


r 1 


(2)  Angle  Class  2 - Mean  wave  approach  angle  within  30  deg 
to  either  side  of  normal  to  shore. 

(3)  Angle  Class  3 - Mean  wave  approach  angle  greater  than 
30  deg  to  left  of  normal  to  shore. 

Originally,  a finer  stratification  of  wave  angles  was  at- 
tempted; however,  it  was  determined  that  even  though  the 
wind  might  be  blowing  parallel  to  shore,  the  mean  angle  of 
approach  shoreward  almost  always  fell  within  60  deg  of 
normal  to  shore.  Thus,  additional  angle  classes  were  not 
necessary. 

i*.  Table  E2  gives  the  significant  wave  heights  for  all  approach 
angles  and  seasons  combined  for  return  periods  of  5,  10,  20,  50,  and 
100  yr. 

5.  Table  E3  wave  period  information  presents  the  data  as  a func- 
tion of  location  (one  per  page),  wave  height  (rows),  and  angle  of  wave 
approach  ( columns ) . 

6.  Table  EU  presents  the  azimuths  of  the  normal  to  shore  vectors 
shown  in  Figure  21. 
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Table  El 


TABLE  OF  EXTHEMES  ESTIMATES 
GRID  LOCATION  22,12  LATM3.Q2  lON=82.33 

SHORELINE  GRID  POINT  l 
WINTER 

angle  classes 


SPRING 

angle  classes 


ANGLE  CLASSES 
2 


(Continued) 


E3 


port  Huron 

Ml 


1 

2 

3 

all 

5 

8 . 9 ( 

1.6) 

19- 7( 

1.6) 

9 . 5 ( 

1.6) 

19,7' 

0.4) 

10 

10. 8( 

2.1) 

20. 3( 

2.1) 

10  • 2 < 

2.1) 

20. 3< 

0.5) 

20 

12. 1( 

2.7) 

21. 0< 

2.7) 

10- 8< 

2,7) 

21, 0( 

0,6) 

50 

13. 4< 

3.3) 

21. 6( 

3.3) 

11. 8< 

3,3) 

2 1 , 6 * 

0,7) 

100 

14. 1< 

3,8) 

22. 0< 

3.8) 

12-11 

3,8) 

22.0' 

0,9) 

1 

2 

3 

all 

5 

2 . 0 ( 

1.5) 

9 . 8 ( 1.5) 

5 • 1 ( 

1.5) 

9,9< 

0,6) 

10 

2 . 8 ( 

2,0) 

11  - 4 < 2.0) 

6 • 3 ( 

2.0) 

11.5* 

0,7) 

20 

3 . 1 ( 

2.5) 

12. 2(  2.5) 

7 . 5 < 

2,5) 

12. 3* 

0,9) 

50 

3 . 5 ( 

3.0) 

13- 4(  3,0) 

8.71 

3.0) 

13,5'. 

1.1) 

iOo 

3 . 9 ( 

3.5) 

14. 2(  3,5) 

9 . 1 < 

3.5) 

14,3' 

1.3) 

SUMMER 

angle  classes 

1 

2 

3 

all 

5 

1 . 3 < 

1,6) 

9.8'  1.6) 

6 * 9 < 

1.6) 

9.8' 

0,4) 

10 

2 . 6 ( 

2.1) 

10. 8(  2.1) 

7 . 9 < 

2,1) 

1 0 . 8 < 

0.6) 

20 

3 . 6 ( 

2.7) 

11. 5<  2.7) 

8 • 9 ( 

2.7) 

11.5* 

0.7) 

i>0 

4 . 6 ( 

3.3) 

12- 1(  3.3) 

9 . 8 < 

3,3) 

12.2': 

0.9) 

100 

5 . 2( 

3.8) 

12. 5(  3.8) 

10. 2< 

3,8) 

12. 6 < 

1.0) 

fall 

all 


5 

6 . 9 < 

1.5) 

17. 7( 

1.5) 

8 . 9 ( 

1.5) 

17.8' 

0.4) 

10 

9 . 8 ( 

2.0) 

19. 0( 

2.0) 

9.8C 

2,0) 

19.1' 

0.6) 

20 

11.8c 

2.5) 

19. 7( 

2.5) 

10-8C 

2.5) 

19,7' 

0.7) 

50 

14. If 

3.0) 

20. 3( 

3.0) 

1 1 . 5 < 

3,0) 

20. 3< 

0,9) 

100 

15. 4( 

3.5) 

2 1 . 0 ( 

3,5) 

11 . 8 < 

3,5) 

2 1 , 0 ' 

1.1) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 
GRID  LOCATION  21,12  LAT=43,16  LQN=82.33 

shoreline  GRID  POINT  2 
WINTER 


lake^ort 

Ml 


angle  classes 


1 

2 

3 

ALl- 

5 

6 . 2 ( 

1.6) 

14. 8( 

1.6) 

18. 4< 

1.6) 

18.5^ 

0.4) 

10 

6 . 9 ( 

2.1) 

16- 4( 

2.1) 

19 . 0 C 

2.1) 

19, 1< 

0.5? 

20 

7 . 5 ( 

2.7) 

17. 7( 

2.7) 

19. 7( 

2.7) 

19, 7< 

0.6) 

50 

8 . 2 ( 

3.3) 

18. 7( 

3.3) 

20  • 3 C 

3.3) 

20.3* 

0.7) 

100 

8 . 5 ( 

o*8) 

19. 4( 

3.8) 

2o  • 7 ( 

3.8) 

20. 7' 

0.9) 

SPRING 

angle  classes 


1 

2 

3 

AL»- 

5 

1 . 6 C 

1.5) 

4 . 7 ( 

1.5) 

8 . 7 ( 

1.5) 

8,8' 

0,6) 

10 

1 .6< 

2.0) 

6 . 7 ( 

2.0) 

10  • 6 < 

2.0) 

10. 7'. 

0.8) 

20 

2 • 0 ( 

2.5) 

8 . 3 ( 

2.5) 

11. 8< 

2.5) 

1 1 . 9 ' 

1.0) 

50 

2 . 0 ( 

3,0) 

9 . 8 ( 

3,0) 

13. O' 

3,0) 

13.lt 

1.2) 

100 

2 . 0 < 

3.5) 

11. 0< 

3.5) 

13. 8( 

3,5) 

13. 9< 

1.4) 

1 


5 

1.0( 

1,6) 

10 

1.0< 

2.1 ) 

20 

1. 0( 

2.7) 

50 

1 . 0 < 

3,3) 

lOQ 

1.0{ 

3,8) 

1 


5 

4 . 9 { 

1.5) 

10 

6 . 2 ( 

2,0) 

20 

6 . 9 ( 

2.5) 

50 

7 . 9 ( 

3.0) 

1°0 

8 . 5 ( 

3.5) 

SUMMER 

angle  classes 
2 3 


3 . 3 ( 

1.6) 

9 • 2 ' 

6 . 6 ( 

2.1) 

9.8< 

8 . 9 ( 

2.7) 

10. 2< 

11 . 5( 

3.3) 

10. 8< 

13. 1( 

3.8) 

11 . 2 ( 

FALL 

angle  classes 

2 

3 

13- 4( 

1.5) 

16. 7( 

16. 4( 

2.0) 

17. 7( 

18. 7( 

2.5) 

18- 4< 

2 1 • 0 ( 

3.0) 

19. O' 

22. 3( 

3.5) 

19. 4( 

1.6) 

alL 

9.3' 

0.4) 

2.D 

10.2' 

0.5) 

2.7) 

10,5' 

0.7) 

3,3) 

11. 6 < 

0,8) 

3,8) 

13.2' 

1.0) 

1.5) 

all 

17. 1< 

0,6) 

2.0) 

18, 4( 

0,8) 

2.5) 

19,4' 

1.0) 

3,0) 

21. 1< 

1.2) 

3,5) 

22,4' 

1.4) 

(Continued ) 
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Table  El  (Continued) 


TABLE  Or  EXTREMES  ESTIMATES 

GRID  LOCATION  20.12  LAT=43,3l  L0N=82.33  LEXINGTON 

Ml 


SHORELINE  GRID  POINT  3 
WINTER 

angle  classes 


1 

2 

3 

alL 

5 

6 . 9 ( 

1.6) 

13. 8( 

1.6) 

13. 1* 

1,6) 

14.  l* 

0.5) 

10 

7.9( 

2.1) 

15. 1( 

2.1) 

13. 8< 

2,1) 

15.1* 

0.6) 

20 

8 . 5 ( 

2.7) 

16. 1( 

2.7) 

14- 1< 

2,7) 

16. 2( 

0.8) 

5>0 

8 . 9 ( 

3,3) 

17. H 

3.3) 

14. 4( 

3,3) 

17 . 2 1 

1.0) 

lOQ 

9 . 5 < 

3.8  ' 

17. 7( 

3.8) 

14. 8( 

3.8) 

17  »8  * 

1.1) 

1 


5 

2 . 4 ( 

1.5) 

10 

2 . 8 ( 

2.0) 

20 

3 . 1 ( 

2.5) 

50 

3 . 5 ( 

3.0) 

100 

3 . 5 ( 

3,5) 

1 


5 

2 . 3 ( 

1,6) 

10 

3 . 3 ( 

2.1) 

20 

3 . 9 ( 

2.7) 

50 

4 . 6 ( 

3.3) 

1°0 

5 . 2 ( 

3.8) 

1 


5 

7.2C 

1.5) 

10 

8 . 5 ( 

2.0) 

20 

9 . 8 ( 

2.5) 

50 

10-8C 

3.0) 

l00 

11.5( 

3.5) 

SPRING 

angle  classes 
2 3 


5.91 

1.5) 

6 - 3 C 

7 . 5 < 

2,0) 

7 • 5 < 

8.71 

2.5) 

8 . 7 < 

9 . 8 ( 

3.0) 

9.8* 

10.2* 

3.5) 

10. 2< 

SUMMER 

angle  classes 

2 

3 

1 • 6 ( 

1.6) 

6 • 2 C 

3 . 0 ( 

2.1) 

7 . 9 C 

3 . 9 ( 

2.7) 

a . 9 < 

4 . 6 ( 

3.3) 

9 . 8 ( 

5 . 2 C 

3,8) 

10  -5( 

FALL 

angle  classes 

2 

3 

13. 1( 

1.5) 

12.  B< 

14  ■ 4( 

2.0) 

14. 8* 

15. 4( 

2.5) 

16. 1( 

16. 4( 

3,0) 

17. 4< 

16.71 

3.5) 

18. 4( 

1.5) 

ALl- 

6,4* 

0,6) 

2,0) 

7,6* 

0,8) 

2,5) 

8,8* 

1.0) 

3.0) 

9,9* 

1,2) 

3,5) 

10,3* 

1.4) 

1.6) 

al(- 

6.3* 

0.7) 

2,1) 

7.9* 

0.9) 

2.7) 

9.0* 

1.1) 

3.3) 

9,8* 

1.4) 

3,8) 

10.5* 

1.6) 

1,5) 

ALC 

13.8* 

0.8) 

2,0) 

15.4* 

1.0) 

2.5) 

16.7* 

1.3) 

3,0) 

18.0* 

1.6) 

3.5) 

19.0* 

1.8) 

(Continued) 
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Table  El  (Continued) 


TABLE  or  EXTREMES  ESTIMATES 
GRID  LOCATION  19.11  LAT=43,45  l0N=82.52 

SHORELINE  GRID  POINT  4 
WINTER 

angle  CLASSES 

12  3 


5 

10. 5( 

1.6) 

14. 4< 

1.6) 

13. 8( 

1,6) 

14.81 

0,3) 

10 

11. 5( 

2.1) 

15. 4( 

2.1) 

14.4* 

2.1) 

15.4* 

0,4) 

20 

12.1* 

2,7) 

15. 7( 

2.7) 

14. 8( 

2.7) 

15, 7* 

0,5) 

50 

12. 8( 

3,3) 

16. 4( 

3,3) 

15.1* 

3,3) 

1 6 * 4 1. 

0.7) 

l00 

13. 1( 

3,8) 

16. 7( 

3.8) 

15. 1* 

3,8) 

16, 7 < 

0.8) 

SPRING 

ANGLE  CLASSES 

1 

2 

3 

alU 

5 

3 . 5 ( 

1,5) 

7. 1( 

1.5) 

5 • 9 ( 

1.5) 

7.2* 

0,5) 

10 

5 . 1 ( 

2.0) 

8 . 7 ( 

2.0) 

7.5( 

2,0) 

8.8* 

0,7) 

20 

5.9< 

2.5) 

9 . 4 ( 

2.5) 

8 • 3 ( 

2,5) 

9.5* 

0.9) 

50 

7.1< 

3.0) 

10 . 6( 

3,0) 

9 . 4 < 

3.0) 

10.7* 

1.1) 

100 

7 . 9 ( 

3,5) 

1 1 . 4 ( 

3.5) 

9 . 8 ( 

3.5) 

11.5* 

1,2) 

SUMMER 

ANGLE  CLASSES 

1 

2 

3 

AL*. 

5 

2 . 6 ( 

1.6) 

2 • 3 ( 

1.6) 

6 • 9 ( 

1.6) 

6.9* 

0,6) 

10 

4 . 3 ( 

2.1) 

4 . 3 < 

2.1) 

8.2* 

2,1) 

8.2* 

0,9) 

20 

5 . 2 ( 

2.7) 

5 • 6 ( 

2.7) 

9-2* 

2,7) 

9.3* 

1.1) 

50 

6 . 2 ( 

3.3) 

6 . 9 ( 

3.3) 

10. 2* 

3,3) 

10.2*. 

1.3) 

100 

7 . 2 ( 

3,8) 

7 • 5 ( 

3.8) 

10  • 8 ( 

3,8) 

10.8* 

1.5) 

fall 

angle  CLASSES 

1 2 3 ALL 


POBT  SANILAC 
Ml 

ALl- 


5 

9 . 8 ( 

1.5) 

13. 8( 

1.5) 

13.4* 

1.5) 

14,8* 

0,9) 

10 

11.2* 

2,0) 

15. 4( 

2,0) 

15.7* 

2.0) 

16.7* 

1,2) 

20 

12. 5( 

2,5) 

16. 7( 

2.5) 

17.4* 

2.5) 

1 8 » 4 * 

1,5) 

50 

13. 4( 

3,0) 

18-0* 

3,0) 

19.4* 

3.0) 

20.0* 

1.8) 

100 

14.1* 

3,5) 

18. 7( 

3.5) 

20.3* 

3,5) 

2l,0< 

2.1) 

( Continued) 
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Table  El  (Continued) 


table  or  extremes  estimates 

GRID  LOCATION  18.11  LATM3.59  l0N=82,52  FOReST^LLE 

Mi 


SHORELINE  GRID  POINT  5 
WINTER 

angle  classes 


1 

2 

3 

ALL 

5 

10. 8( 

1.6) 

14. 8(  1.6) 

14. 4< 

1.6) 

15, 4* 

0,3) 

10 

12. 1( 

2.1) 

15.7<  2.1) 

15. 4( 

2.1) 

16. 1* 

0.4) 

20 

12.81 

2.7) 

16.41  2.7) 

16  • 1 1 

2.7) 

16.4* 

0,5) 

50 

13. 8( 

3.3) 

17.11  3,3) 

17. 1( 

3.3) 

17.lt 

0,7) 

100 

14. 1( 

3.8) 

17.41  3.8) 

17.41 

3,8) 

17.4* 

0.8) 

SPRING 

angle  classes 

1 

2 

3 

alI- 

5 

5 . 1 ( 

1.5) 

7 . 1 ( 1.5) 

6.3< 

1.5) 

7.2* 

0,6) 

10 

5 . 9 ( 

2.0) 

8 . 7 ( 2,0) 

7 • 9 < 

2,0) 

8.8* 

0,8) 

20 

6 . 7 ( 

2.5) 

9 . 8 ( 2,5) 

9* 1< 

2.5) 

9.9* 

1.0) 

50 

7 . 5( 

3.0) 

11. 0<  3,0) 

10. 2( 

3.0) 

H.lt 

1.2) 

100 

7 . 9( 

3.5) 

11- 8(  3.5) 

11- 0< 

3.5) 

11.9* 

1.4) 

SUMMER 

ANGLE  CLASSES 

1 

2 

3 

al<- 

5 

2 . 6 < 

1.6) 

3 • 9 ( 1,6) 

6 • 6 < 

1.6) 

6*6* 

0,7) 

10 

3.9( 

2.1) 

5 . 2 C 2,1) 

8 . 2< 

2.D 

8 , 2 1 

1,0) 

20 

4 , 9 ( 

2,7) 

6 . 2 ( 2,7) 

9 • 5 * 

2,7) 

9f  6 1 

1,2) 

*0 

5.9( 

3,3) 

7 . 2 < 3,3) 

10- 5( 

3,3) 

10- 5t 

1.5) 

lOo 

6 . 6 ( 

3.8) 

?. 9 ( 3.8) 

11. 2< 

3,8) 

11,2* 

1.7) 

FALL 

ANGLE  CLASSES 

1 

2 

3 

ALl- 

5 

10. 2( 

1.5) 

1 4 '•  8 ( 1.5) 

14. 1< 

1.5) 

15,7* 

1.0) 

10 

11. 8< 

2.0) 

16. 7(  2,0) 

16. 7( 

2.0) 

18. Ot 

1.3) 

20 

13. 1( 

2,5) 

18. 0(  2,5) 

18. 7( 

2,5) 

19,7* 

1.7) 

50 

14. 1( 

3.0) 

19. 4<  3.0) 

20  • 7< 

3,0) 

21,3* 

2.1) 

100 

15.1< 

3.5) 

20« 3(  3,5) 

21, 6< 

3,5) 

22,3*. 

2,4) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 
GRID  LOCATION  17,11  LAT*43,74  LQN=82.52  HELENA 

M1 


shoreline  GRID  POINT  6 
WINTER 

angle  classes 


1 

2 

3 

aLU 

5 

11. 8( 

1.6) 

15- 4( 

1.6) 

14. 8( 

1,6) 

15.4* 

0,4) 

10 

12. 1( 

2.1) 

16. 4( 

2,1) 

15. 7( 

2.1) 

16.4* 

0.6) 

20 

12. 8( 

2.7) 

17. 1( 

2.7) 

16. 4( 

2,7) 

17.1* 

0,7) 

so 

13. 1< 

3,3) 

18. 0( 

3.3) 

17. 1( 

3,3) 

18.1* 

0,9) 

100 

13. 4( 

3.8) 

18. 4( 

3.8) 

17- 4( 

3,8) 

18.5* 

1,0) 

SPRING 

angle  classes 


1 

2 

3 

all 

5 

5. 1C 

1.5) 

7.1* 

1.5) 

5.9* 

1,5) 

7.2* 

0.5) 

10 

6 . 3 < 

2.0) 

8 . 7 ( 

2.0) 

7.5* 

2,0) 

8.8* 

0.7) 

20 

7 . 1 ( 

2,5) 

9 • 8 ( 

2.5) 

8.7* 

2.5) 

9.9* 

0,9) 

so 

7 . 9 ( 

3.0) 

11.  oc 

3.0) 

9.8* 

3.0) 

11.1* 

1.1) 

100 

8 . 7 ( 

3,5) 

11.8* 

3.5) 

10-6* 

3.5) 

11,9* 

1.2) 

SUMMER 

ANGLE  CLASSES 

1 

2 

3 

alL 

5 

3 . 6 ( 

1.6) 

3-9* 

1,6) 

6.6* 

1.6) 

6.6* 

0.8) 

10 

4 . 3 < 

2.1) 

5-2* 

2.1) 

8-2* 

2.1) 

8.2* 

1.1) 

20 

5 . 2 ( 

2.7) 

5.9* 

2.7) 

9.5* 

2.7) 

9,6* 

1.3) 

so 

5.9C 

3.3) 

6.6* 

3.3) 

10-8* 

3,3) 

10.8* 

1.7) 

100 

6.2C 

3,8) 

7.2* 

3.6) 

11.5* 

3,8) 

11.5* 

1.9) 

fall 

ANGLE  CLASSES 

1 

2 

3 

AlL 

5 

1 0 . 2 ( 

1.5) 

14.8* 

1,5) 

13.8* 

1,5) 

16.1* 

0.9) 

10 

12. 1( 

2.0) 

16. 7( 

2,0) 

16-7* 

2,0) 

18.0* 

1.2) 

20 

13.8C 

2.5) 

18. OC 

2.5) 

18. 7C 

2,5) 

19.7* 

1,5) 

50 

15. 1( 

3.0) 

19.7  C 

3,0) 

20. 7( 

3,0) 

21,3* 

1,9) 

100 

16. 1C 

3.5) 

20  - 3 C 

3.5) 

22.0* 

3,5) 

22,3* 

2.2) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATE* 
GRID  LOCATION  16.11  LATM3.88  LQN=82.52 

SHORELINE  GRID  POINT  7 
WINTER 


HARBOR  BEACH 
Ml 


angle  classes 


1 

2 

3 

al>- 

5 

10. 2( 

1,6) 

16. 1( 

1.6) 

16. 1< 

1,6) 

16, 7*. 

0.4) 

10 

11.8C 

2.1) 

17. 11 

2.1) 

16. 7( 

2,1) 

17.71 

0.5) 

20 

13. 1( 

2.7) 

18  - 0 < 

2.7) 

17  • 4 < 

2,7) 

18.4* 

0,7) 

50 

14. 4< 

3.3) 

19. 0( 

3.3) 

18. 0( 

3,3) 

19.1* 

0,8) 

100 

15. 1( 

3.8) 

19. 4( 

3.8) 

18. 4( 

3,8) 

19.4*: 

1.0) 

spring 

ANGLE  CLASSES 

1 2 3 AlL 


5 

5 . 9 ( 

1.5) 

6 • 3 ( 

1.5) 

7. 1C 

1,5) 

7.2* 

0,6) 

10 

7 . 9 ( 

2.0) 

8 . 7 ( 

2,0) 

9.1* 

2,0) 

9.2* 

0.8> 

20 

9 . 1 ( 

2.5) 

10  • 2 ( 

2.5) 

10 . 6 < 

2,5) 

10.7* 

1.0) 

50 

10. 6( 

3.0) 

11. 8< 

3.0) 

11. 8( 

3,0) 

11.9* 

1.3) 

100 

11 ,4( 

3.5) 

13  - 0 < 

3.5) 

12  • 6 < 

3,5) 

13.1* 

1.5) 

SUMMER 

ANGLE  CLASSES 

1 

2 

3 

all 

5 

1.6( 

1.6) 

3 • 6 < 

1.6) 

8 • 2 ( 

1.6) 

8,9*. 

0,6) 

10 

2 . 3 ( 

2.1) 

6.2* 

2.1) 

9.5* 

2,D 

10.2* 

0.9) 

20 

3 . 0 ( 

2.7) 

8.2( 

2.7) 

10. 5( 

2,7) 

H.2* 

l.D 

50 

3 . 6 ( 

3.3) 

10  • 2 ( 

3.3) 

11 . 5C 

3,3) 

12,1* 

1.3) 

l00 

3 . 9 ( 

3.8) 

11*5( 

3.8) 

12>1< 

3,8) 

12.8* 

1.5) 

fall 

angle  classes 

1 

2 

3 

al^ 

5 

8 . 2 ( 

1.6) 

16- 7( 

1.5) 

14. 8< 

1,5) 

17.1* 

0.9) 

10 

ii.ec 

2.2) 

19. 0< 

2.0) 

16.1* 

2.0) 

19,1* 

1.2) 

20 

14. 8( 

2.7) 

20. 7< 

2.5) 

17. 1< 

2.5) 

20.8* 

1.5) 

50 

17. 4( 

3.3) 

22. 3 ( 

3.0) 

18.0* 

3,0) 

22,4* 

1.8) 

100 

19. 0( 

3.9) 

23. 3( 

3.5) 

18. 7( 

3.5) 

23,4* 

2,1) 

( Continued) 
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Table  El  (Continued) 


TABLE  of  EXTREMES  ESTIMATES 
GRID  LOCATION  15, iQ  LAT=44,03  LQN=82,71 

shoreline  GRID  POINT  8 
winter 

angle  CLASSES 


heron  CITY 
Ml 


1 

2 

3 

alL 

5 

15. 7( 

1,6) 

18. 0( 

1.6) 

16. K 

1.6) 

18.1* 

0.5) 

10 

17. 1( 

2,1) 

19. 0< 

2.1) 

17*  7( 

2,1) 

19. 1* 

0,6) 

20 

17. 7( 

2.7) 

20- 0< 

2.7) 

19. 4< 

2.7) 

20. 0* 

0.8) 

50 

18. 4( 

3,3) 

20. 7( 

3.3) 

20  • 7 < 

3,3) 

20. 7* 

1,0) 

100 

19. 0( 

3.8) 

21. 3< 

3.8) 

2l . 6 ( 

3,8) 

2 1 . 7 < 

l.D 

SPRING 

angle  CLASSES 

1 

2 

3 

alI- 

5 

5 . 5 ( 

1.5) 

8 • 3 ( 

1.5) 

7 . 5 ( 

1,5) 

8 • 4 « 

0.5) 

10 

7 . 5 ( 

2.0) 

9 . 4 ( 

2.0) 

9 • 1 C 

2.0) 

9 . 5 C 

0,6) 

20 

9 . 1 ( 

2.5) 

10  • 6 ( 

2.5) 

1 0 • 2 C 

2,5) 

1 0 • 7 < 

0.8) 

50 

10. 2( 

3,0) 

11. 4( 

3.0) 

11 . 4 < 

3,0) 

U.5* 

1.0) 

l00 

11. 0( 

3.5) 

12. 2( 

3.5) 

12. 2( 

3,5) 

12. 3*. 

1.1) 

SUMMER 

ANGLE  CLASSES 

1 

2 

3 

ALU 

5 

2 . 3 ( 

1.6) 

9 . 5 ( 

1.6) 

9 • 8 C 

1,6) 

10.5* 

0,8) 

10 

4 . 3 ( 

2.1) 

11. 2( 

2,1) 

11. 8( 

2,1) 

12*1*. 

1.0) 

20 

5 . 9 ( 

2,7) 

12- 5( 

2.7) 

13. 1C 

2,7) 

13. 4< 

1,3) 

50 

7 . 5 ( 

3.3) 

13- 4( 

3.3) 

14 . 4 C 

3,3) 

14.4* 

1.6) 

100 

8 . 5 ( 

3,8) 

14. 1( 

3,8) 

15. 1< 

3,8) 

15.4* 

1.8) 

F ALL 

ANGLE  CLASSES 

1 

2 

3 

AL*> 

5 

14. 4( 

1.5) 

16. 7( 

1.5) 

15. 4( 

1.5) 

17.4* 

0,7) 

10 

16. 7( 

2.0) 

18. 7( 

2.0) 

16. 7< 

2,0) 

19.4* 

1.0) 

20 

18. 4( 

2.5) 

20  - 3 C 

2.5) 

17. 7C 

2.5) 

20.3* 

1.2) 

50 

20. 0< 

3.0) 

22- 0( 

3.0) 

19. 0C 

3.0) 

22.1* 

1.5) 

100 

21. 0< 

3.5) 

23- 0< 

3.5) 

19. 7< 

3.5) 

23.1* 

1.8) 

(Continued) 
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Table  F.1  (Continued) 


table  or  extremes  estimates 

GRID  LOCATION  15.  8 LATs44.03  L0N  = 8-J.11 

SHORELINE  GRID  point  9 
WINTER 

angle  classed 

12  3 


5 

16. 4( 

1.6) 

18. 0( 

1.6) 

11. 5C 

1.6) 

10 

17. 4( 

2.1) 

19. 4C 

2,1) 

12. 8< 

2.1) 

20 

18. 0( 

2.7) 

20. 31 

2.7) 

14. 1C 

2.7) 

50 

18. 7( 

3,3) 

21. 3( 

3.3) 

15. 1< 

3,3) 

lOo 

19. 0( 

3.8) 

22. 0( 

3.8) 

15. 7C 

3,8) 

SPRING 

angle  classes 

1 

2 

3 

5 

6 . 3 ( 

1.5) 

9. 1C 

1.5) 

7 • 5 C 

1.5) 

IQ 

8.3C 

2.0) 

10- 2( 

2.0) 

8 . 3 < 

2.0) 

20 

9 . 4 ( 

2.5) 

11  - 0 ( 

2.5) 

9. 1C 

2,5) 

50 

11. 0( 

3,0) 

11. 8C 

3,0) 

9 . 8 C 

3,0) 

100 

1 1 - 8 ( 

3.5) 

12  • 6 ( 

3.5) 

10. 2C 

3.5) 

SUMMER 

ANGLE  classes 

1 

2 

3 

5 

5 . 2 ( 

1,6) 

9 • 8 ( 

1.6) 

8 • 9 C 

1.6) 

10 

7.9( 

2.1) 

11  - 5 < 

2.1) 

9 . 5 < 

2.D 

20 

9.8C 

2.7) 

12- 5( 

2.7) 

10  • 2 C 

2.7) 

50 

12. 1C 

3.3) 

13. 8( 

3.3) 

1 0 * 8 ( 

3,3) 

lO0 

13. 1C 

3.8) 

14. 4< 

3.6) 

H • 2C 

3,8) 

PALL 

ANGLE  classes 

1 

2 

3 

5 

15.4C 

1.5) 

17. 1C 

1.5) 

10- 2C 

1.5) 

10 

17. 7C 

2.0) 

17. 7C 

2.0) 

ll . 8 < 

2,0) 

20 

19. 4( 

2.5) 

18. 4< 

2.5) 

13  • 1 < 

2.5) 

50 

21. OC 

3.0) 

19. 0( 

3.0) 

14. 4C 

3,0) 

lOQ 

22.3c 

3.5) 

19. 4( 

3.5) 

15. 1C 

3,5) 

(Continued) 


PORT  CRESCENT 

Ml 


all 

18. 1<  0.6) 
19.41  0,7) 
2 0 . 3 < 0.9) 
21.31  1.2) 
22. 1<  1.3) 


ALL 

9 .2  1.  0.4) 
10. 3(  0,5) 
il.lt;  0.7) 
11.91  0.8) 

12.71  1.0) 


all 

1 0 » 8 I 0,4) 
ll , 8 1 0.6) 
12.61  0.7) 
13. 9«  0.9) 
14.51  1.0) 


ALL 

17. ll  0,7) 

18.71  0.9) 

19. 7 1 1.1) 
21. Ol  1.4) 
22.4«  1.6) 
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Table  El  (Continued) 


TABLF  of  EXTREMES  ESTIMATES 

GRID  LOCATION  15,  7 LAT*44.Q3  L0N=83.30  ENTRANCE  SAGINAW  BAY 

M* 


shoreline  GRID  POINT  10 
WINTER 

angle  classes 


1 

2 

3 

all 

5 

18. 0( 

1.6) 

18. 0( 

1.6) 

10. 2C 

1,6) 

18.7* 

0.6) 

10 

19. 0( 

2.1) 

19. 7( 

2.1) 

11. 2( 

2,1) 

20,0( 

0.8) 

20 

19. 7( 

2.7) 

21. 0< 

2.7) 

12. 1( 

2,7) 

21, 0‘ 

1.0) 

50 

20. 3( 

3.3) 

22- 6( 

3.3) 

12. 8< 

3,3) 

22, 7* 

1.2) 

lOQ 

20. 7( 

3.8) 

23. 3( 

3.6) 

13. 1< 

3,8) 

23 . 4 « 

1,4) 

1 

5 

6 . 7 ( 

1.5) 

10 

8 . 7 ( 

2.0) 

20 

9 . 8 ( 

2.5) 

5o 

11. 0( 

3.0) 

100 

11  .Sc 

3.5) 

1 


5 

4 . 6 ( 

1.6) 

10 

6 . 9 ( 

2.1 ) 

20 

8 . 5 ( 

2.7) 

50 

10. 2( 

3.3) 

lOQ 

11. 2( 

3.8) 

1 


5 

17. It 

1.5) 

10 

18. 4< 

2.0) 

20 

19. 7( 

2.5) 

50 

20. 7( 

3.0) 

lOQ 

21, 3( 

3.5) 

SPRING 

angle  classes 

2 3 


7 . 1 ( 1.5) 

5.1( 

8 . 7 < 2.0) 

6 « 7 ( 

9 . 8 < 2.5) 

7 • 9 ( 

11. 4(  3.0) 

9 . 1C 

12. 2(  3.5) 

9 . 8 ( 

SUMMER 

angle  CLASSES 

2 

3 

8 . 2 ( 1.6) 

6 > 6 ( 

10- 2(  2,1) 

8-2( 

11. 2(  2.7) 

9.2< 

12 • 5 ( 3.3) 

10. 2< 

13. 4(  3.8) 

10- 8* 

fall 

angle  classes 

2 

3 

16. 4(  1.5) 

11 . 8 C 

18. 7(  2.0) 

13  • 4 ( 

20- 3(  2.5) 

14. 8< 

22. 3 ( 3.0) 

15. 7( 

23. 3(  3.5) 

16. 7< 

AlL 

1,5) 

7,2* 

0,5) 

2,0) 

8,8  < 

0.6> 

2.5) 

9,9* 

0,8) 

3.0) 

11.5* 

1.0) 

3.5) 

12.3* 

l.D 

alu 

1.6) 

8,5* 

0,8) 

2.D 

10.2* 

1.0) 

2.7) 

11.5* 

1.3) 

3,3) 

12,6* 

1.6) 

3.8) 

13,5* 

1.8) 

alL 

1.5) 

18,0* 

0.7) 

2.0) 

19, 7< 

0.9) 

2.5) 

20,7* 

1.2) 

3,0) 

22, 4^ 

1.4) 

3,5) 

23,4* 

1,7) 

(Continued) 


(Sheet  10  of  28) 


El  2 


Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 

GRID  LOCATION  14.  7 LAT=44.17  L0N=83.3Q  TARAS  CITY 

Ml 


SHORELINE  GRID  POINT  U 
WINTER 

angle  classes 


1 

2 

3 

alL 

5 

9 . 2 ( 

1.6) 

15. 7( 

1.6) 

19. 4( 

1,6) 

19.4* 

0,4) 

10 

10. 2( 

2.1) 

16- 4< 

2.1) 

20- 3* 

2,1) 

20, 3* 

0,6) 

20 

10. 8( 

2,7) 

17. 1( 

2.7) 

21. 0* 

2,7) 

21. 0* 

0,7) 

50 

11.51 

3,3) 

17. 7( 

3.3) 

2l . 6 C 

3,3) 

2 1 * 6 *' 

0.9) 

lO0 

11 .81 

3.8) 

18. 0( 

3.8) 

22*  0* 

3,8) 

22. 0* 

1.0) 

SPRING 

angle  classes 

1 

2 

3 

AlL 

5 

3 . 9 ( 

1.5) 

7 . 5 ( 

1.5) 

6. 7( 

1.5) 

7,6* 

0,7) 

10 

5.1< 

2.0) 

9 . 4 ( 

2.0) 

8 • 3 < 

2.0) 

9,5* 

1.0) 

20 

5 . 9 ( 

2.5) 

11- 4 ( 

2.5) 

9 • 4 ( 

2.5) 

11, 5* 

1.2) 

50 

6 . 7 ( 

3.0) 

13. 0< 

3,0) 

10-61 

3,0) 

13.1* 

1.5) 

100 

7 . 1 ( 

3.5) 

14- 2( 

3.5) 

U-0* 

3.5) 

1 4 . 3 * 

1.8) 

SUMMER 

angle  classes 


1 

2 

3 

ALl 

5 

3 . 9 ( 

1.6) 

4 . 3 ( 

1.6) 

5 • 6( 

1.6) 

6.2< 

0.7) 

10 

5 . 2 ( 

2.1) 

5 . 9 ( 

2.1) 

7 • 2 * 

2.D 

7.5* 

0.9) 

20 

6 . 2 ( 

2.7) 

6 • 9 ( 

2.7) 

8 • 5 ( 

2.7) 

8.5* 

1.1) 

50 

7 . 2 < 

3.3) 

7 • 9 ( 

3.3) 

9.8< 

3,3) 

9.8* 

1.4) 

lOo 

7 . 9 ( 

3.8) 

8 . 5 ( 

3.8) 

10- 5< 

3.8) 

10.5* 

1,6) 

FAlL 

ANGLE  CLASSES 

1 

2 

3 

AL*- 

5 

9 . 2 ( 

1.5) 

14. 1( 

1.6) 

17. 1( 

1.5) 

17.1* 

0.9) 

10 

11 .2( 

2.0) 

17. 7< 

2.1) 

19. 4( 

2,0) 

19.4* 

1.2) 

20 

12. 5( 

2.5) 

2 0 • 3 ( 

2.6) 

2 0 • 7 ( 

2,5) 

21.0* 

1.5) 

5 0 

14. 1( 

3.0) 

23. 3( 

3.3) 

22  • 3 ( 

3,0) 

23.4* 

1.9) 

lOO 

14 ,8( 

3.5) 

24. 9( 

3.8) 

23  • 3 < 

3,5) 

25.0* 

2.2) 

(Continued) 
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Table  El  (Continued) 


table  or  extremes  estimates 

GRID  LOCATION  13*  8 LAT=44.3l  LON=03.11  OSCGDA  Au  SAbLfc 

Ml 


SHORELINE  grid  point  12 
WINTER 

ANGLE  CLASSES 


1 

2 

3 

AlL 

5 

11. 5( 

1.6) 

16. 4 ( 

1.6) 

18. 0( 

1.6) 

18,  lC 

0.6) 

10 

12. 5( 

2.1) 

17. 1* 

2.1) 

19  - 4 < 

2.D 

1 9 . 4 * 

0.7) 

20 

13. 1( 

2.7) 

17. 7( 

2.7) 

20 . 3 ( 

2,7) 

2 0 » 3 * 

0,9) 

50 

13. 8( 

3.3) 

1 8 . 0 ( 

3,3) 

2l.3< 

3,3) 

2l,4* 

1.2) 

lO0 

14. 1( 

3.8) 

18. 4( 

3.8) 

2i.6< 

3,8) 

21.6* 

1,3) 

spring 

ANGLE  CLASSES 

1 

2 

3 

al*- 

5 

5 . 1 ( 

1.5) 

7 . 9 ( 

1.5) 

5. 1C 

1.5) 

a.o* 

0,6) 

10 

6 . 7 ( 

2,0) 

9 . 8 ( 

2.0) 

7.5* 

2.0) 

9.9* 

0.8) 

20 

7 . 5 ( 

2.5) 

11. 0< 

2.5) 

9.4* 

2,5) 

H.  1* 

1.0) 

50 

8 . 7 ( 

3.0) 

12  • 6 ( 

3.0) 

11. 4< 

3,0) 

12. 7< 

1.2) 

lOQ 

9 . 1 < 

3.5) 

13. 8( 

3.5) 

12- 2( 

3,5) 

13.9* 

1,4) 

SUMMER 

angle  classes 


1 

2 

3 

all 

5 

5 . 6 < 

1.6) 

4 • 9 ( 

1.6) 

4.3* 

1,6) 

5.9* 

0,6) 

10 

6 . 6 ( 

2.1) 

5 • 9 ( 

2.1) 

6-2* 

2.D 

7.2* 

0,9) 

20 

7 . 2 ( 

2.7) 

6.6* 

2.7) 

7.5* 

2,7) 

8.2* 

1.1) 

50 

7 . 9 ( 

3,3) 

7 ■ 5 C 

3.3) 

S . 9 C 

3,3) 

9.5* 

1,3) 

l00 

8 . 2 ( 

3.8) 

7 . 9 < 

3.8) 

9 . 8 ( 

3,8) 

10-2* 

1.5) 

FALL 

ANGLE  CLASSES 

' 1 

2 

3 

AL*- 

5 

11. 8< 

1.5) 

15  ■ 1 ( 

1.5) 

16. 4( 

1,5) 

16.7* 

0,9) 

10 

13. 8( 

2.0) 

18. 0( 

2.0) 

18-0* 

2.0) 

19,0* 

1.2) 

20 

15. 1( 

2.5) 

2 0 • 3 C 

2.5) 

19.4* 

2,5) 

20,3* 

1.5) 

50 

16. 4( 

3,0) 

22. 6( 

3.0) 

20  • 7< 

3,0) 

22.7* 

1.8) 

l00 

17. 4( 

3.5) 

24. 3( 

3.5) 

2l . 3 ( 

3,5) 

24.4* 

2.1) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 
GRID  LOCATION  12.  8 LAT=44,46  LON'S.*,!! 

SHORELINE  GRID  POINT  it 


6REENBUSH 

Ml 


WINTER 

angle  classes 


1 

2 

3 

all 

5 

12. 8< 

1.6) 

16. 4( 

1.6) 

18. 0( 

1.6) 

18  *4  ( 

0.4) 

10 

13. 4( 

2.1) 

17. 41 

2.1) 

19. 4( 

2.D 

19, 4< 

0,5) 

20 

14. 1( 

2.7) 

18. 0( 

2.7) 

20  • 0 < 

2.7) 

20»0< 

0,7) 

50 

14. 8( 

3.3) 

18. 7( 

3.3) 

2l  • 0 < 

3.3) 

21, 1( 

0,8) 

100 

15. 4( 

3,8) 

19. 0< 

3.8) 

2l.3< 

3,8) 

21, 41. 

1,0) 

SPRING 

angle  classes 

1 

2 

3 

all 

5 

5 . 9 ( 

1.5) 

7 . 5 ( 

1.5) 

3 • 9 ( 

1.5) 

7,6< 

0,6) 

10 

7 . 1 ( 

2.0) 

9 . 4 ( 

2.0) 

6 • 3 < 

2.0) 

9.5< 

0.8) 

20 

8 . 3 ( 

2.5) 

11  - 0 < 

2.5) 

7 • 9 ( 

2,5) 

lltl* 

1.0) 

50 

9 . 4 ( 

3.0) 

12  - 6 ( 

3.0) 

9.8< 

3,0) 

12»7< 

1.2) 

100 

9 . 8 ( 

3.5) 

13. 8( 

3.5) 

10. 6( 

3,5) 

13, 9( 

1.4) 

SUMMER 

angle  classes 

1 

2 

3 

ALl- 

5 

5 . 6 ( 

1.6) 

5 . 2 ( 

1.6) 

4 • 3 ( 

1.6) 

6.2* 

0,6) 

10 

6 . 2 ( 

2.1) 

6 • 9 ( 

2.1) 

6 • 2 t 

2,1) 

7r9‘ 

1,0) 

20 

6 . 9 ( 

2.7) 

8 • 2 ( 

2.7) 

7 . 9 ( 

2.7) 

9.2< 

1,3) 

50 

7 . 5 ( 

3.3) 

9-2( 

3,3) 

9 • 2 < 

3,3) 

10, 5< 

1.6) 

100 

7 . 9 ( 

3.8) 

10  • 2 ( 

3.6) 

10«2( 

3,8) 

U,2< 

1.8) 

fall 

ANGLE  CLASSES 

1 

2 

3 

All- 

5 

13. 1( 

1.5) 

14. 4( 

1.5) 

16. K 

1.5) 

16 . 7 < 

0,9) 

10 

14. 8( 

2.0) 

18- 0< 

2,0) 

17. 7( 

2.0) 

18, 7< 

1.2) 

20 

15. 7( 

2.5) 

20  * 3 ( 

2.5) 

19. 0< 

2,5) 

20, 3‘ 

1.5) 

50 

17. 1( 

3.0) 

23. 0< 

3.1) 

2 0 • 0 < 

3,0) 

23. 1< 

1,9) 

lOO 

17. 7( 

3.5) 

24. 6< 

3,6) 

2l.  0< 

3,5) 

24, 7< 

2.2) 

(Continued) 
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Table  El  (Continued) 


TAPLE  OF  EXTREMES  ESTIMATES 
GRID  LOCATION  11,  8 LAT=44,oO  L0N=83,il 

SHORELINE  GRID  POINT  14 


HARRlSviLLfc 

Ml 


WINTER 

angle  classes 


1 

2 

3 

all 

5 

13. 11 

1.6) 

19.  DC 

1.6) 

14. 8< 

1.6) 

19,1* 

0,6) 

10 

14. 8< 

2.1) 

20- OC 

2.1) 

16. 4< 

2.D 

20.0* 

0.7) 

20 

15. 7( 

2.7) 

21 . 0( 

2.7) 

17. 4( 

2.7) 

21,0* 

0.9) 

50 

16. 7( 

3.3) 

22.0< 

3,3) 

18  • 7 ( 

3,3) 

22.0* 

1,2) 

100 

17. 7( 

3.8) 

22- 3( 

3.8) 

19. 4( 

3,8) 

22.3* 

1,3) 

SPRING 

ANGLE  CLASSES 

1 

2 

3 

AL*- 

5 

6 . 7 ( 

1.5) 

a . 3 ( 

1.5) 

5 • 1 ( 

1.5) 

8.4* 

0.7) 

10 

7 . 5 ( 

2.0) 

10- 2( 

2.0) 

7.1( 

2,0) 

10,3* 

0.9) 

20 

7 . 9 ( 

2.5) 

11  - 4 < 

2.5) 

8 . 7 ( 

2.5) 

11.5* 

1.1) 

50 

8 . 7 ( 

3.0) 

13 . 0 < 

3.0) 

1 0 • 2 < 

3,0) 

13.1* 

1.4) 

mo 

8 . 7 ( 

3.5) 

13. 8( 

3.5) 

11. 0< 

3.5) 

13.9* 

1.6) 

SUMMER 

ANGlE  CLASSES 

1 

2 

3 

ALL 

5 

3 . 9 ( 

1.6) 

5 . 9 ( 

1,6) 

6 . 2 C 

1,6) 

6»6* 

0.8) 

10 

5 . 2 ( 

2.1) 

7.2< 

2.1) 

8 • 2 ( 

2,D 

8,2* 

1.1) 

20 

6 . 6 ( 

2.7) 

8 . 2 ( 

2.7) 

9 . 8 ( 

2,7) 

9.9* 

1.3) 

5 0 

7 . 5 ( 

3.3) 

9 . 2 < 

3,3) 

11. 2< 

3,3) 

1 1 . 3 < 

1.7) 

100 

8 . 2 ( 

3.8) 

9 . 8 ( 

3.8) 

12  • 1 < 

3,8) 

12.2* 

1.9) 

fall 

ANGlE  classes 

1 

2 

3 

all 

5 

11. 2( 

1.5) 

17. 1( 

1.5) 

14. 1< 

1.5) 

17.7* 

0.8) 

10 

13. 1( 

2.0) 

19. 4< 

2.0) 

17. 1< 

2,0) 

19.7* 

1.1) 

20 

14. 4( 

2.5) 

21- 0( 

2.5) 

19. 0( 

2,5) 

2l.O* 

1,3) 

50 

16.  It 

3.0) 

22- 6( 

3.0) 

21. 0* 

3,0) 

22,7* 

1,7) 

l00 

17. 1( 

3.5) 

23  - 6 ( 

3.5) 

22 . 3 < 

3,5) 

23,7* 

1.9) 

(Continued) 
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Table  El  (Continued) 


table  OF  EXTREMES  estimates 
GRID  LOCATION  10,  8 lAT=44,74  L0N=83,ll 

SHORELINE  GRID  POINT  15 
U INTER 

angle  classes 


black  River 

Ml 


1 

2 

3 

A Li- 

5 

13. 1( 

1.6) 

18. 7< 

1.6) 

14-8< 

1.6) 

18  . 8 i 

0.5) 

10 

14.  It 

2.1) 

19. 7( 

2.1) 

16.11 

2.1 ) 

19,71 

0.6) 

20 

14. 8( 

2.7) 

20.71 

2.7) 

17.11 

2.7) 

20.7* 

0,8) 

50 

15. 4( 

3.3) 

21. 3( 

3.3) 

18.01 

3,3) 

2 1 . 3 1 

1.0) 

100 

16. 1( 

3,8) 

22.01 

3,6) 

18.41 

3,8) 

22.0* 

1,1) 

SPRING 

ANGLE  CLASSES 


1 

2 

3 

ALl 

5 

7 . 1 ( 

1.5) 

8 . 3 ( 

1.5) 

4.71 

1.5) 

8 . 4 < 

0,7) 

10 

7 . 9 ( 

2.0) 

10.21 

2.0) 

6.71 

2,0) 

10.31 

0.9) 

20 

8 . 3 < 

2.5) 

11 . 8 ( 

2.5) 

7.91 

2,5) 

11.91 

1.2) 

5 0 

9 . 1 ( 

3.0) 

13. 4( 

3.0) 

9.41 

3,0) 

13.51 

1.4) 

100 

9 . 4 ( 

3.5) 

14.21 

3.5) 

10.21 

3,5) 

14.31 

1.7) 

SUMMER 

angle  classes 


1 

2 

3 

all 

5 

3.9C 

1.6) 

5.91 

1.6  ) 

5-91 

1.6) 

6 ,6 1 

0.8) 

10 

5.21 

2.1) 

6-91 

2.1) 

7.91 

2,1) 

8 . 2 < 

1.0) 

20 

5 . 9 ( 

2.7) 

7.91 

2.7) 

9.5< 

2.7) 

9 , 6 1 

1.3) 

^0 

6 . 9 ( 

3.3) 

8.51 

3.3) 

11-21 

3.3) 

11. 3« 

1,6) 

100 

7 . 5 ( 

3.8) 

9-21 

3.8) 

12-1< 

3,8) 

12. 2( 

1,8) 

fall 

ANGLE  CLASSES 

1 

2 

3 

all 

5 

12.11 

1,5) 

17.11 

1.5) 

12-81 

1,6) 

17, 7< 

0.8) 

10 

13.81 

2.0) 

19.41 

2.0) 

16-41 

2,2) 

19.71 

1.1) 

- u 

15.11 

A C « 

C.  • J ! 

n , ■*  • 

Cl  • 0 V 

2,5) 

iy . 4 i 

2.7) 

2l,4l 

1,4) 

50 

16.41 

3,0) 

23.31 

3.0) 

22-01 

3.3) 

23,41. 

1,7) 

1O0 

17.41 

3,5) 

24.31 

3.5) 

23-61 

3,9) 

24 » 4 * 

2.0) 

(Continued) 
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Table  El  (Continued) 


table  OF  EXTREMES*  estimates 

grid  LOCATION  9,  8 LATs44,89  LQN=83,11  oSSlNEKE 

M* 


SHORELINE  GRID  POINT  16 
winter 

angle  CLASSES 


1 

2 

3 

al*- 

5 

13.  «< 

1.6) 

18. 7( 

1.6) 

14. 4< 

1.6) 

18.8* 

0.6) 

10 

14. 4( 

2.1) 

20- 0* 

2.1) 

15. 4< 

2,1) 

2 0 . 0 (. 

0.9) 

20 

15. 1( 

2.7) 

21. 0< 

2.7) 

16. 4< 

2,7) 

2l,0< 

1.1) 

*0 

15. 4< 

3.3) 

22. 0( 

3.3) 

17. 1( 

3,3) 

22.0* 

1.31 

lOO 

16. 1( 

3.8) 

22- 6( 

3.8) 

17. 7( 

3,8) 

2 2 * 6 * 

1.5) 

SPRING 

ANGLE  CLASSES 


1 

2 

3 

all 

5 

7 . 5 ( 

1.5) 

7 . 9 < 

1.5) 

4 ■ 7 ( 

1,5) 

8,0* 

0.7) 

10 

9 . 1 ( 

2.0) 

10. 2< 

2,0) 

6 • 3 ( 

2,0 

10,3* 

1,0) 

20 

10. 2( 

2.5) 

11- 8< 

2.5) 

7.9  ( 

2.5) 

11,9* 

1,2) 

50 

11. 4< 

3.0) 

13. 0* 

3.0) 

9.K 

3,0) 

13.1* 

1,5) 

lO0 

12. 2( 

3.5) 

14.21 

3.5) 

9 . 8 ( 

3,5) 

14,3* 

1,8) 

1 


5 

4 . 9 ( 

1.6) 

10 

5 . 6 ( 

2.1) 

20 

6 . 2 ( 

2.7) 

50 

6 . 9 ( 

3.3) 

100 

7 . 2 ( 

3.8) 

1 


5 

11. 5( 

1.5) 

10 

13. 8( 

2,0) 

20 

i5 . 4 ( 

2.5) 

50 

17. 1( 

3.0) 

lOQ 

18. 0( 

3.5) 

SUMMER 

angle  classes 

2 3 


5 . 2 ( 

1.6) 

5 • 6 C 

6 . 6 ( 

2.1) 

7.5( 

7 . 2 ( 

2,7) 

9.2* 

8 . 2 ( 

3.3) 

10  -8< 

8 . 5 ( 

3 . b ) 

H.8< 

fall 

angle  CLASSES 

2 

16. 4( 

1.5) 

3 

12-1* 

19. 4( 

2.0) 

15- 7( 

21- 3( 

2.5) 

18. 4( 

23. 3( 

3.0) 

21. 3( 

24. 6( 

3.5) 

23-0* 

A'L*- 

1.6) 

6.2* 

0.8) 

2,1) 

7.9* 

1.1) 

2,7) 

9.3* 

1,3) 

3,3) 

10,9* 

1,7) 

3,8) 

11,9* 

1,9) 

ALL 

1.6) 

17.1* 

1,0) 

2.2) 

19.7* 

1.4) 

2,7) 

21,3* 

i,7i 

3,3) 

23,4* 

2.1) 

3.9) 

24.7* 

2.5) 

(Continued) 
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Table  El  (Continued) 


table  or  extremes  estimates 

GRID  LOCATION  8.  8 LAT=45.Q4  L0N  = 8«S.ll  NORTH  >*0INT 

Ml 


SHORELINE  GRID  POINT  17 
WINTER 

ANGLE  CLASSES 


1 

2 

3 

Ml 

5 

13. 1( 

1.6) 

18. 4< 

1.6) 

12. K 

1,6) 

18,5* 

0.6) 

10 

14. 4( 

2.1) 

19. 7( 

2.1) 

14. 8< 

2,D 

19,7* 

0,9) 

20 

15. 1( 

2,7) 

20. 7< 

2.7) 

17. 1< 

2,7) 

20.7* 

1.1) 

50 

16. 1( 

3.3) 

21-61 

3.3) 

19. 0< 

3,3) 

21, 6* 

1.3) 

100 

16. 7( 

3,8) 

22- 3( 

3.8) 

2q. 3( 

3,8) 

22,3* 

1.5) 

SPRING 

angle  classes 

1 

2 

3 

ALL 

5 

7 . 5 ( 

1.5) 

7 . 9 ( 

1.5) 

4.31 

1.5) 

8 . 0 I 

0.7) 

10 

9 . 4 ( 

2,0) 

9 • 8 < 

2.0) 

5 • 9 < 

2.0) 

9 , 9 1 

1.0) 

20 

1 0 . 6 ( 

2.5) 

11. 4< 

2.5) 

7 . 5 ( 

2.5) 

11. 5< 

1.2) 

5*0 

11. 8( 

3.0) 

13. 0( 

3,0) 

9.1< 

3.0) 

13,1* 

1,5) 

1°0 

12. 6( 

3.5) 

13- 8( 

3.5) 

9.8< 

3,5) 

13, 9< 

1.8) 

summer 

ANGLE  CLASSES 

1 

2 

3 

all 

5 

4 . 6 < 

1.6) 

4 . 9 ( 

1.6) 

4 • 6 < 

1,6) 

5,6  * 

0.8) 

10 

4 . 9 ( 

2.1 ) 

6 • 2 < 

2.1) 

6 • 9 ( 

2.D 

7.2< 

1.0) 

20 

5 . 2 ( 

2.7) 

7 . 5 < 

2.7) 

8 . 9 C 

2,7) 

9,0* 

1,3) 

50 

5 . 6 < 

3,3) 

8 . 5 ( 

3.3) 

10. 5( 

3,3) 

10.6* 

1.6) 

lOo 

5 . 9 ( 

3,8) 

9 . 2 < 

3.8) 

11 . 8 C 

3,8) 

11.9* 

1.8) 

FALL 

angle  classes 

1 

2 

3 

alL 

5 

11. 5( 

1.5) 

16. 1( 

1.5) 

11. 5< 

1.6> 

16.7* 

1.1) 

10 

13. 8( 

2.0) 

19. 4( 

2.0) 

15. 1( 

2.2) 

19, 4< 

1.5) 

2n 

l<5.4f 

2,5) 

21  .M 

2.5) 

17. 7< 

2.7) 

2 1 .71 

1,9) 

50 

17. 4( 

3.0) 

23  ( 

3.0) 

20. 7( 

3,3) 

24io« 

2.4) 

100 

18. 4( 

3.5) 

25. 3( 

3.5) 

22. 3( 

3,9) 

25, 4< 

2.7) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 
GRID  LOCATION  7.  7 LAT*45.18  L°N=83.30 

SHORELINE  GRID  POINT  18 
WINTER 

angle  classes 

12  3 


5 

14.  i( 

1,6) 

18. 0( 

1.6) 

13,1* 

1,6) 

10 

15. 1( 

2.1) 

19. 7( 

2.1) 

14.8* 

2,1) 

20 

15. 4( 

2.7) 

21.0* 

2,7) 

16.1* 

2,7) 

50 

16. 1( 

3.3) 

22- 3( 

3.3) 

17. 4( 

3,3) 

l00 

16. 4( 

3.8) 

23.0( 

3.8) 

18. 4< 

3,8) 

SPRING 

angle  classes 

1 

2 

3 

5 

6 . 7 ( 

1.5) 

5 • 5 ( 

1.5) 

5.1* 

1.5) 

10 

9 . 1 ( 

2.0) 

6 . 7 ( 

2.0) 

7.1* 

2,0) 

20 

10. 6( 

2.5) 

7 . 5 ( 

2.5) 

8 . 7 ( 

2,5) 

50 

12. 2( 

3.0) 

8 • 3 ( 

3.0) 

9.8* 

3,0) 

lOQ 

13. 0( 

3.5) 

9.1* 

3.5) 

11.0* 

3.5) 

summer 

ANGLE  CLASSES 

1 

2 

3 

5 

3 . 6 ( 

1.6) 

5-2* 

1.6) 

6.6< 

1,6) 

10 

4 . 6 < 

2.1) 

7 . 5 < 

2.1) 

7.9* 

2.1) 

20 

5 . 2 < 

2.7) 

9.2( 

2.7) 

8 . 9 ( 

2,7) 

^0 

5 . 9 ( 

3.3) 

1 0 - 8 * 

3.3) 

9.8* 

3,3) 

100 

6 . 2 ( 

3.8) 

11.8* 

3.8) 

10. 5( 

3,8) 

FALL 

ANGlE  CLASSES 

1 

2 

3 

5 

12. 5( 

1.8) 

16. 7( 

1.5) 

11-2* 

1,5) 

10 

16. 4( 

2.4) 

18. 4( 

2.0) 

11.8* 

2,0) 

20 

19. 71 

2.9) 

20  - 0 C 

2.5) 

12.5* 

2.5) 

50 

22. 6( 

3.7) 

21- 3( 

3.0) 

13.1* 

3,0) 

lOQ 

24. 3( 

4.2) 

22  - 3 * 

3.5) 

13.4* 

3,5) 

(Continued) 


Rockport 

Ml 


all 

18. 1*  0.8) 
19, 7*  1,0) 
21, 1,3) 
22,4*  1,6) 
23,0*  1.8) 


all 

6,8*.  0.7) 
9,2*  0.9) 
10.7*  1,1) 
12.3*  1,4) 
13.1*  1.6) 


ALL 

7.2*  0.8) 
8.9*  1.0) 
10,2*  1.3) 
11.5*  1.6) 
12.1*  1.8) 


all 

16.7*  0,8) 
18,7*  1,1) 
20,0*  1,4) 
22.7*  1.7) 
24.4*  2.0) 
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Table  FI  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 

GRID  LOCATION  6.  7 LATM5.32  LON=83.30  STONE  PORT 


shoreline  GRID  point  19 

RINTER 

angle  classes 


1 

2 

3 

ALl- 

5 

15. 1( 

1.6) 

13. 4( 

1.6) 

15.  K 

1,6) 

16.4* 

0.3) 

10 

17. 1( 

2.1) 

15. 1( 

2.1) 

16. 1* 

2,1) 

17.1* 

0,4) 

20 

18. 4< 

2.7) 

16.4( 

2.7) 

17- 1< 

2,7) 

18.5* 

0,5) 

50 

19. 7( 

3.3) 

17. 7( 

3.3) 

17- 7< 

3,3) 

19.8* 

0.7) 

100 

20 .3( 

3.8) 

18. 4( 

3.8) 

18.  Q< 

3,8) 

2 0 , 4 < 

0,8) 

SPRING 

ANGLE  CLASSES 

1 

2 

3 

alL 

5 

5 . 1 ( 

1.5) 

3 • 5 ( 

1.5) 

7*1* 

1.5) 

7.2* 

0,7) 

10 

6.7( 

2.0) 

6 . 3 ( 

2.0) 

9 . 4 ( 

2,0) 

9,5* 

0.9) 

20 

7 . 9 ( 

2.5) 

7 . 9 ( 

2.5) 

11.  o* 

2.5) 

11.1* 

1,2) 

50 

9 . 1 ( 

3.0) 

9. 4< 

3.0) 

12  • 6 < 

3.0) 

12.7* 

1.4) 

100 

9 . 8 ( 

3.5) 

10 .6( 

3.5) 

13. 4( 

3.5) 

13,5* 

1.7) 

SUMMER 

angle  classes 


1 

2 

3 

AL*. 

5 

2 . 3 ( 

1.6) 

3 • 3 < 

1.6) 

9.2* 

1,6) 

9.5* 

0,6) 

10 

4 . 3 ( 

2.1) 

6 . 2 < 

2.1) 

10- 5( 

2.1 ) 

10.5* 

0.7) 

20 

5 . 6 ( 

2.7) 

8 . 2 ( 

2.7) 

11. 5< 

2.7) 

11.5* 

0.9) 

50 

6 . 9 ( 

3.3) 

10. 2( 

3.3) 

12. 5< 

3.3) 

12.6* 

1.2) 

100 

7 . 9 ( 

3.8) 

11 .5( 

3 , b ) 

13. 1( 

3.8) 

13.2* 

1.3) 

FALL 

ANGLE  CLASSES 

1 

2 

3 

al*- 

5 

13. 4( 

1.5) 

12.51 

1.5) 

13  • 8 ( 

1,5) 

15.1* 

0.5) 

10 

15. 1( 

2.0) 

15.  K 

2.0) 

14. 1( 

2.0) 

16.1* 

0.6) 

2 n 

16 . 1 ( 

2 . 5 ) 

1 7 . i ( 

2,5) 

14- 4( 

O C * 
fc  * ' ' 

AMI* 

A r>  % 

U . O t 

50 

17. it 

3.0) 

19. 4< 

3.0) 

15-1* 

3.0) 

19.5* 

1.0) 

1 0 0 

17. 7( 

3.5) 

2 0 • 7 ( 

3.5) 

15-1* 

3,5) 

2 0 • 8 < 

1.1) 

(Continued) 
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Table  El  (Continued) 


TABLE  of  EXTREMES  ESTIMATES 
GRID  LOCATION  5,  6 LAT=45.46  L0N=83.50 

shoreline  GRID  POINT  20 

WINTER 

ANGLE  CLASSES 


adams  Point 

Ml 


1 

2 

3 

all 

5 

13 .1< 

1.6) 

12- 1( 

1.6) 

13. 8* 

1.6) 

14.4* 

0.4) 

10 

15. 1( 

2.1) 

13. 8( 

2.1) 

15. K 

2.1) 

15. 4* 

0,5) 

2 0 

16. 4( 

2.7) 

15. 1( 

2.7) 

16.  K 

2.7) 

16. 5* 

0.7) 

so 

17. 7( 

3.3) 

16- 4( 

3.3) 

1 6 7 * 

3.3) 

17.8* 

0,8) 

lOO 

18. 4( 

3.8) 

17.  It 

3,8) 

17. 4( 

3,8) 

18-5* 

1.0) 

SPRING 

angle  classes 


1 

2 

3 

all 

5 

4 . 7 ( 

1.5) 

3 • 9 ( 

1.5) 

6.3* 

1.5) 

6.4* 

0.6) 

10 

5 . 9 ( 

2.0) 

5 • 9 ( 

2.0) 

8.7* 

2.0) 

8.8* 

0.8) 

20 

6 . 7 < 

2.5) 

7.5* 

2.5) 

10-2* 

2.5) 

10.3* 

1.0> 

50 

7 . 9 ( 

3.0) 

8 • 7 ( 

3,0) 

11.8* 

3.0) 

11,9* 

1,2) 

1*>0 

8 . 3 ( 

3.5) 

9 . 4 C 

3.5) 

12.6* 

3,5) 

12.7* 

1.4) 

SUMMER 

angle  classes 


1 

2 

3 

all 

5 

2 . 3 ( 

1.6) 

2.0* 

1.6) 

8.5* 

1,6) 

8.5* 

0,5) 

10 

3 . 9 < 

2.1) 

4.3* 

2.1) 

9.5* 

2.1) 

9.6* 

0,7) 

20 

5 . 2 ( 

2,7) 

6.2* 

2.7) 

10-5* 

2.7) 

10.5* 

0.9) 

5o 

6 . 6 ( 

3.3) 

7.9* 

3.3) 

11.2* 

3.3) 

11.2* 

l.D 

1*>0 

7 . 2 ( 

3.8) 

8 . 9 ( 

3.8) 

11.8* 

3,8) 

11,9* 

1.2) 

fall 

angle  CLASSES 


1 

2 

3 

ALL 

5 

12.1* 

1,5) 

11.5* 

1.5) 

12.8* 

1,5) 

13.8* 

0.3) 

10 

13. 4( 

2.0) 

13.8* 

2.0) 

13.4* 

2.0) 

14.8* 

0.4) 

20 

14.4* 

2.5) 

15 . 4 ( 

2.5) 

14.1* 

2.5) 

15,4* 

0.6) 

50 

15. 4( 

3.0) 

17.4* 

3.0) 

14.4* 

3,0) 

17.5* 

0,7) 

100 

16. 1( 

3,5) 

18.4* 

3.5) 

14.6* 

3.5) 

18,5* 

0.8) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTKEMEb  ESTIMATES 
GRID  LOCATION  5.  5 LAT=45.46  LON=83,70 

shoreline  GRID  POINT  21 
WINTER 

angle  classes 


ROSERS  city 
Ml 


1 

2 

3 

AL*- 

5 

14. 1< 

1.6) 

13. 4( 

1.6) 

15. 7( 

1.6) 

16.  ll 

0,4) 

10 

15. 4( 

2.1) 

15. 1( 

2.1) 

16. 7( 

2.D 

16, 7< 

0,5) 

20 

16. 4( 

2.7) 

16. 71 

2.7) 

17- 4( 

2.7) 

17,5l 

0.6) 

50 

17. 4( 

3.3) 

18.01 

3.3) 

18- 4( 

3,3) 

1 8 . 5 < 

0.7) 

lO0 

18. 0( 

3.8) 

19.01 

3.8) 

18. 7( 

3.8) 

19, 1< 

0.9) 

SPRING 

ANGLE  classes 

1 

2 

3 

ale 

5 

5 . 1 ( 

1.5) 

4 . 3 < 

1.5) 

6 . 7 ( 

1.5) 

6.8‘ 

0.6) 

10 

6 . 7 ( 

2.0) 

6 . 3 ( 

2.0) 

9.1< 

2.0) 

9 ,2  l- 

0,7) 

20 

8 . 3 ( 

2,5) 

7 . 9 < 

2.5) 

10. 2( 

2,5) 

10. 3< 

0,9) 

50 

9 . 4 ( 

3.0) 

9 . 4 ( 

3.0) 

11. 8< 

3,0) 

11 .9  C 

1.1) 

100 

10 .2< 

3.5) 

10 .2( 

3.5) 

13.0< 

3,5) 

1 3 * 1 ^ 

1.3) 

SUMMER 

angle  classes 

1 

2 

3 

all 

5 

3 . 0 ( 

1.6) 

3 . 3 ( 

1.6) 

9 • 8 ( 

1,6) 

9.8< 

0.4) 

10 

4 . 9 { 

2.1 ) 

5 . 9 ( 

2,1) 

1 0 • 8 ( 

2,1) 

10. 8< 

0,5) 

20 

6.2( 

2.7) 

7 -9  ( 

2.7) 

11. 5( 

2.7) 

11.51 

0,7) 

50 

7 . 5 ( 

3.3) 

9 . 8 ( 

3.3) 

12  • 1 1 

3,3) 

12, 2l 

0.8) 

100 

8 . 2 ( 

3.8) 

11. 2C 

3,8) 

12  • 5 ( 

3,8) 

12. 6* 

1.0> 

FALL 

ANGLE  classes 

1 

2 

3 

all 

5 

12. 8( 

1.5) 

12. 8( 

1.5) 

14. 4( 

1.5) 

15, 4( 

0.2) 

10 

14. 8( 

2.0) 

14. 8( 

2.0) 

15. 1< 

2,0) 

15, 7 l 

0.3) 

20 

16. 4( 

2.5) 

16  • 4 i 

2.5) 

15  • 4 ( 

2,5) 

At  4 f 

40  • ' 

0,4) 

50 

18. 0( 

3.0) 

18. 0( 

3.0) 

15. 7( 

3,0) 

18, 1* 

0,5) 

100 

19. 0( 

3.5) 

19. 0( 

3.5) 

16- 1( 

3,5) 

19. 1< 

0.5) 

(Continued) 
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Table  El  (Continued) 


table  of  extremes  estimates 

GRID  LOCATION  4.  4 LAT=45.61  LON  = 8«J.90 

SHORELINE  GRID  POINT  22 
WINTER 

angle  classes 


hammoNd  bay 
m! 


1 

2 

3 

al<- 

5 

12. 5( 

1.6) 

13- 1( 

1.6) 

14. 1( 

1.6) 

1 4 . 8 I. 

0.4) 

10 

13. 8( 

2.1) 

14. 4( 

2.1) 

15. 4( 

2.1) 

15.41 

0.5) 

20 

14  ,B( 

2.7) 

15- 4( 

2.7) 

16- 4< 

2.7) 

16. 5* 

0,6) 

50 

15. 7( 

3.3) 

16. 4( 

3.3) 

17. 4( 

3.3) 

17,5* 

0.7) 

1 o o 

16. 1( 

3.8) 

17. 1( 

3.8) 

18. 0* 

3,8) 

18,1* 

0.9) 

spring 

angle  classes 

1 

2 

3 

alI- 

5 

4 . 7 ( 

1.5) 

3 • 5 ( 

1.5) 

6 • 3 ( 

1.5) 

6.4* 

0 . 5 > 

10 

5 . 9 ( 

2.0) 

5 . 5 ( 

2.0) 

8 . 7 ( 

2.0) 

8.8( 

0.6) 

20 

7 . 1 ( 

2.5) 

6 - 7 < 

2.5) 

10*2  < 

2.5) 

10.3< 

0,8) 

5 0 

7 . 9 ( 

3.0) 

8 • 3 ( 

3,0) 

11. 8< 

3,0) 

H.9l 

1.0) 

100 

a . 7 ( 

3.5) 

9.1( 

3.5) 

12. 6( 

3.5) 

1 2 . 7 * 

1.1> 

SUMMER 

ANGLE  classes 

i 

2 

3 

al*- 

5 

2 . 6 ( 

1.6) 

3 . 3 ( 

1.6) 

9 . 2 C 

1.6) 

9,3* 

0.3) 

10 

3 . 9 < 

2.1) 

5 . 6 ( 

2.1) 

9 • 8 C 

2.1) 

9.8* 

0.4) 

20 

4 . 9 ( 

2.7) 

7 . 2 ( 

2.7) 

10. 5( 

2.7) 

10.5* 

0.5) 

50 

5 . 9 ( 

3.3) 

8 • 9 ( 

3.3) 

10. 8( 

3,3) 

10.8* 

0.7) 

100 

6 . 6 ( 

3.8) 

10. 2( 

3.8) 

11 . 2 ( 

3,8) 

11,2* 

0,8) 

fall 

ANGLE  CLASSES 

1 

2 

3 

AL*- 

5 

11. 2< 

1.5) 

12.  It 

1.5) 

13  • 8 ( 

1.5) 

14,4*. 

0,2) 

10 

13. 4( 

2.0) 

13  - 8 ( 

2.0) 

14  • 4 < 

2.0) 

14,8* 

0.2) 

20 

14. 8( 

2.5) 

14. 8( 

2.5) 

15. 1< 

2.5) 

15.1* 

0.3) 

5 0 

16. 4( 

3.0) 

15. 7( 

3.0) 

15. 4( 

3.0) 

16.5* 

0.4) 

100 

17. 4( 

3.5) 

16. 4( 

3.5) 

15- 7( 

3,5) 

17.5* 

0.4) 

(Continued) 
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Table  El  (Continued) 


TABLE  or  EXTREMES  ESTIMATES 

GRID  LOCATION  3.  3 LAT*45.75  lQN=84,09  CORCWOOW  POINT 

Ml 


SHORELINE  GRID  POINT  23 
WINTER 

angle  classes 


1 

2 

3 

AlG 

5 

11. 2< 

1.6) 

12 .5  ( 

1.6) 

11. 8( 

1,6) 

13, 4* 

0.2) 

10 

11. 8( 

2.1) 

13. 8( 

2.1) 

13. 1* 

2.D 

13, 8* 

0,3) 

20 

12. 5( 

2.7) 

14. 8( 

2.7) 

14. 1* 

2,7) 

14,9* 

0,4) 

so 

13. 1< 

3.3) 

15. 4( 

3.3) 

15. 1* 

3,3) 

15.5*. 

0,5) 

100 

13. 4( 

3.8) 

16. 1* 

3.8) 

15- 7< 

3,8) 

16. 2* 

0,6) 

SPRING 

angle  classes 

1 

2 

3 

all 

5 

4 . 3 ( 

1.5) 

3 • 1 ( 

1.5) 

5.5( 

1,5) 

5,6* 

0.5) 

10 

5.1* 

2.0) 

4 . 7 ( 

2.0) 

7 . 5 < 

2,0) 

7-6* 

0,7) 

20 

5.5* 

2.5) 

5 • 9 ( 

2.5) 

8. 7< 

2,5) 

8.8* 

0.9) 

50 

6 . 3 ( 

3.0) 

7 . 1 ( 

3,0) 

1 0 • 2 < 

3.0) 

10.3* 

1.1) 

l00 

6 . 3 ( 

3.5) 

7.5< 

3.5) 

11  - 0 « 

3.5) 

ll-l* 

1.2) 

SUMMER 

angle  classes 

1 

2 

3 

all 

5 

2 . 3 ( 

1.6) 

4.61 

1.6) 

7 . 9 ( 

1.6) 

8.2* 

0.3) 

10 

3.9( 

2.1) 

6 • 9 ( 

2.1) 

8.5C 

2.1) 

8 . 9 I 

0,4) 

20 

4 . 9 ( 

2.7) 

8 • 5 ( 

2.7) 

9 • 2 < 

2.7) 

9.3* 

0.5) 

50 

5 . 9 ( 

3.3) 

10. 5( 

3.3) 

9 . 5 < 

3,3) 

10,6* 

0,6) 

l00 

6 . 6 ( 

3.8) 

11 . 5 < 

3 , b ) 

9 . 8 ( 

3,8) 

11.6* 

0.7) 

FALL 

angle  classes 

1 

2 

3 

AL*- 

5 

10. 2( 

1.5) 

ll-5t 

1.5) 

12- 1( 

1.5) 

12,8*. 

0.3) 

10 

11. 5( 

2.0) 

12- 5( 

2.  U) 

13. 1< 

2,0) 

13.4* 

0,3) 

20 

12. 5( 

2.5) 

13- 1( 

2.5) 

13. 8( 

2.5) 

13.8* 

0.4) 

=>0 

13.  M 

3.0) 

13  • 8 i 

O , V ; 

AM  ft  / 

X H > O 1 

3,0) 

14.9* 

A r t 

U . 2 / 

too 

14.1* 

3,5) 

14. 1< 

3.5) 

15. 1< 

3,5) 

15, 2< 

0,6) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXT«E*Ei>  ESTIMATES 
GRID  LOCATION  2,  3 LAT=45.89  lON=84.09 

SHORELINE  GRID  POINT  24 
WINTER 

angle  classes 


POINT  DOLOMITE 
Ml 


1 

2 

3 

ALl 

5 

9 . 5 ( 

1.6) 

15- 7( 

1.6) 

19  • 5( 

1,6) 

15, 7* 

0.4) 

10 

1? . 5 ( 

2.1) 

16-  4( 

2.1) 

11- 2( 

2.1> 

16. 4* 

0.5) 

20 

14 . 4( 

2.7) 

17. 1( 

2.7) 

ll-8( 

2.7) 

17.1* 

0.6) 

so 

16. 7( 

3.3) 

17. 7( 

3.3) 

12  • 5 ( 

3,3) 

17. 8* 

0,7) 

ion 

18. 0( 

3.8) 

18  - 0 ( 

3.8) 

12- 8< 

3,8) 

18.  1* 

0,9) 

SPRING 

ANGLE  CLASSES 

1 

2 

3 

all 

5 

4 . 7 ( 

1.5) 

9 • 4 ( 

1.5) 

1 - 2 < 

1.5) 

9.5* 

0.7) 

10 

5 . 5 ( 

2.0) 

11.41 

2.0) 

1 . 6 < 

2,0) 

U.5* 

0.9) 

20 

6 . 3 ( 

2.5) 

13-  0( 

2.5) 

2-  0( 

2,5) 

13, 1* 

1.2) 

so 

6 . 7 ( 

3.0) 

14.21 

3.0) 

2 • 4 ( 

3,0) 

14. 3* 

1.4) 

100 

7 . 1 ( 

3.5) 

15-  0( 

3.5) 

2 • 8 < 

3.5) 

15.1* 

1.7) 

SUMMER 

angle  classes 

1 

2 

3 

AlU 

5 

4 . 6 ( 

1.6) 

6 > 2 ( 

1,6) 

1 • 0 c 

1,6) 

6.9* 

0,4) 

10 

6 . 6 ( 

2.D 

7 . 5 ( 

2.1) 

1 . 3 * 

2,1) 

7*5* 

0.6) 

20 

7 . 9 < 

2.7) 

8 • 5( 

2.7) 

1 . 3 ( 

2,7) 

8.5* 

0.7) 

so 

9 , 2 ( 

3.3) 

9 . 2 ( 

3.3) 

1 • 6 ( 

3,3) 

9.3*' 

0.9) 

lOQ 

10 ,2( 

3.8) 

9 . 8 < 

3.8) 

1 . 6 < 

3,8) 

10.3* 

1.0) 

FALL 

ANGLE  CLASSgS 


1 

2 

3 

al>- 

5 

7 . 9 ( 

1.5) 

15- 7( 

1.5) 

7 • 9 ( 

1,5) 

15.7* 

1.0) 

10 

10 .5( 

2.0) 

18 . 0 ( 

2.0) 

8 • 9 ( 

2,0) 

18.1* 

i.3) 

20 

12. 5( 

2.5) 

19 . 7 ( 

2.5) 

9.8< 

2.5) 

19.7* 

1.7) 

so 

14. 4( 

3.0) 

21 . 3 C 

3.0) 

IQ . 5 ( 

3,0) 

21. 3^ 

2.1) 

lOQ 

15. 7( 

3,5) 

22  - 3 < 

3.5) 

11-2* 

3.5) 

22.3* 

2.4) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTIMATES 

GRTD  LOCATION  2.  4 LAT=45,89  LQN=83,90  DETOUR  REEF 

Ml 


SHORELINE  GRID  POINT  25 
WINTER 

angle  CLASSES 


1 

2 

3 

ALt- 

5 

7 . 9 ( 

1.8) 

15. 7(  1.6) 

13. 8( 

1.6) 

15,7* 

0.2) 

10 

11. 5< 

2.3) 

16. 4(  2.1) 

14. 8< 

2,1> 

16, 4* 

0.3) 

20 

14. 4< 

2.9) 

16. 7(  2.7) 

15- 7( 

2.7) 

16, 7 * 

0.4) 

50 

17. 1< 

3,6) 

17. 1(  3.3) 

16. 4( 

3,3) 

17, 1* 

0,5) 

100 

18. 7( 

4,2) 

17. 4(  3.8) 

17. 1( 

3,8) 

18, 8* 

0,6) 

SPRING 

angle  classes 

1 

2 

3 

all 

5 

5.9C 

1.5) 

8 • 3 ( 1.5) 

7 • 9 ( 

1.5) 

8.4* 

0,8) 

10 

7 . 1 ( 

2.0) 

1 0 - 2 C 2.0) 

9 . 8 ( 

2,0) 

10.3* 

1.0) 

20 

8 . 3 ( 

2.5) 

12. 2(  2.5) 

11. 0< 

2,5) 

12.3* 

1.3) 

50 

9 . 1 ( 

3.0) 

13. 8(  3.0) 

12. 6( 

3,0) 

13.9* 

1.6) 

1°0 

9 . 4 ( 

3,5) 

15. 0(  3.5) 

13. 4( 

3.5) 

15.1* 

1,8) 

SUMMER 

angle  classes 

1 

2 

3 

AL*. 

5 

6 . 9 ( 

1.6) 

8 . 2 < 1.6) 

2 > 6 ( 

1,6) 

8,2* 

0.4) 

10 

8 . 5 ( 

2,1) 

8 • 9 ( 2.1) 

4.61 

2.1) 

9.0* 

0.6) 

20 

9 . 8 ( 

2.7) 

9 . 5 ( 2.7) 

6 • 2 < 

2,7) 

9,9*: 

0,7) 

50 

11. 2< 

3.3) 

10- 2(  3.3) 

7 • 5 ( 

3,3) 

11,3* 

0.9) 

lO0 

1 1 . 8 C 

3.8) 

10 • 8 ( 3.8) 

8 • 5 ( 

3,8) 

11.9* 

1.0) 

f all 

angle  CLASSES 

1 

2 

3 

ALL 

5 

10. 8< 

1.5) 

15. 7(  1.5) 

14. 1( 

1.5) 

16.1* 

0,6) 

10 

12. 'U 

2.0) 

16. 7{  2,0) 

16. 7( 

2.0) 

17.4* 

0.8) 

20 

14. 4( 

2.5) 

17. 4(  2.5) 

18  • 4 C 

2.5) 

18,5* 

1,0) 

50 

16. 1( 

3.0) 

18. 0(  3,0) 

2 0 • 3 C 

3,0) 

20, a* 

1.3) 

lOQ 

17. 1( 

3.5) 

13 . 7 ( 3,5) 

2 1 • 3 < 

3.5) 

2 1 » 4 * 

1.5) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTATES 

GRID  LOCATION  3.  5 LATe45.75  lON=83,70  WEST  ENO  DRUMMOND  ISLAND 

M l 


SHORELINE  GRID  POINT  26 
WINTER 

angle  classes 


1 

2 

3 

all 

5 

10 .8< 

1.6) 

15. 1( 

1.6) 

13. 4C 

1,6) 

15, 4* 

0,2) 

10 

13. 8( 

2.1) 

15. 7( 

2.1) 

14. 8< 

2.1) 

16.  l* 

0,3) 

20 

15. 7( 

2.7) 

16. 41 

2.7) 

15. 7( 

2,7) 

16. 4* 

0,4) 

50 

18. 0( 

3.3) 

16. 7( 

3.3) 

16. 4( 

3,3) 

18.1* 

0.5) 

100 

19. 4< 

3.8) 

17.  l( 

3.8) 

17. 1( 

3,8) 

19.5* 

0,6) 

SPRING 

angle  classes 


1 

2 

3 

A L 

5 

6 . 3 ( 

1.5) 

7 • 9 { 

1.5) 

7 • 5 < 

1,5) 

8*0  < 

0,7) 

10 

7 . 1 ( 

2,0) 

10- 2( 

2.0) 

9 . 4 ( 

2,0) 

10,3* 

1.0) 

20 

7 . 9 ( 

2.5) 

ii  - at 

2.5) 

11. 0< 

2,5) 

11*  9* 

1.2) 

50 

8 . 7 ( 

3.0) 

13. 8( 

3.0) 

12. 6< 

3.0) 

13.9* 

1.5) 

100 

9 . 1 ( 

3,5) 

15. 0( 

3.5) 

13- 4< 

3.5) 

15.1* 

1,8) 

SUMMER 

angle  classes 

1 

2 

3 

ALL 

5 

6.9C 

1.6) 

7 . 5 ( 

1.6) 

3 ■ 0 < 

1,6) 

7,5* 

0,5) 

10 

8 . 5 ( 

2.1) 

8 • f ( 

2.1) 

4 • 9 ( 

2,1) 

8,5*. 

0.7) 

20 

9 . 8 ( 

2,7) 

9. 5( 

2.7) 

6 • 2 < 

2,7) 

9 » 9 * 

0,9) 

50 

11. 2( 

3.3) 

10  2( 

3.3) 

7 . 5 < 

3,3) 

11*3* 

1.1) 

l00 

11. 8( 

3.8) 

10  • 8 ( 

3.6) 

8 • 2 ( 

3,8) 

11*  9* 

1.2) 

fall 

angle  classes 

1 2 3 ALU 


5 

11. 2( 

1.5) 

14. 8( 

1.5) 

13. 4( 

1.5) 

15,4* 

0.6) 

10 

13. 1( 

2.0) 

16. 4< 

2.0) 

16. 1( 

2,0) 

17*1* 

0.8) 

20 

14. 8( 

2,5) 

17. 4( 

2.5) 

18. 0( 

2.5) 

18,1* 

1,0) 

50 

16. 4( 

3.0) 

18. 4( 

3.0) 

19. 7< 

3,0) 

19,8* 

1.3) 

lOo 

17. 4( 

3.5) 

18. 7( 

3.5) 

2l  • 0 * 

3,5) 

21.1* 

1.5) 

(Continued) 
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Table  El  (Continued) 


TABLE  OF  EXTREMES  ESTATES 

GRID  LOCATION  3.  6 LAT=45.75  LON=83.50  FALSE  DETOUR  CHANNEL 

Ml 


SHORELINE  GRID  POINT  27 
WINTER 

angle  classes 


1 

2 

3 

alI- 

5 

12. 5( 

1.6) 

16. 1( 

1.6) 

13. 4( 

1,6) 

16,4': 

0.4) 

10 

14. 8( 

2.1> 

17. 4( 

2.1) 

14. 8< 

2,D 

17,4' 

0.5) 

20 

16. 4( 

2.7) 

18. 0< 

2.7) 

15- 7( 

2.7) 

18,1' 

0.6) 

so 

18. 0< 

3.3) 

18. 7( 

3.3) 

16. 7( 

3,3) 

18.8' 

0.7) 

100 

19. 0( 

3.8) 

19-  0( 

3.6) 

17. 4( 

3,8) 

19.1' 

0.9) 

SPRING 

angle  classes 


1 

2 

3 

all 

5 

6 . 7 ( 

1.5) 

9 • 1 ( 

1.5) 

7 . 9 < 

1.5) 

9.2' 

0.7) 

10 

7 . 5 ( 

2.0) 

11. 4( 

2.0) 

9 • 8 ( 

2.0) 

11.5' 

0.9) 

20 

7 . 9 < 

2.5) 

13.01 

2.5) 

11,0' 

2,5) 

13.1' 

1,2) 

50 

8 . 7 < 

3.0) 

15-  0( 

3.0) 

12  • 6 ( 

3,0) 

15.1' 

1.4) 

lOQ 

8 . 7 ( 

3.5) 

15. 7( 

3.5) 

13  • 4 C 

3,5) 

15.8' 

1.7) 

1 


5 

6 . 9 ( 

1.6) 

10 

8 . 9 ( 

2,1) 

20 

10. 5( 

2.7) 

50 

1 1 . 8 < 

3.3) 

1 □ 0 

12. 8( 

3.8) 

1 


5 

14. 1( 

1.5) 

10 

15. 7( 

2.0) 

20 

17. 1( 

2.5) 

50 

18. 4( 

3.0) 

lOQ 

19. 0( 

3.5) 

SUMMER 

angle  classes 
2 3 


7 . 9 ( 

1.6) 

3- O' 

8 • 9 ( 

2.1) 

4.9( 

9 . 5 ( 

2.7) 

6 • 6 < 

10- 5( 

3.3) 

7 • 9 C 

1 0 • 8 ( 

3.6) 

8 • 9 < 

fall 

angle  classes 

2 

3 

15. 4( 

1.5) 

13- 1( 

17. 4( 

2.0) 

16.4' 

18 . 7( 

2,5) 

18. 7( 

20- 0< 

3,0) 

2l . 0 ' 

20  • 7 ( 

3,5) 

22. 3( 

(Continued) 


1.6) 

ALi- 

7.9' 

0.5) 

2,1) 

9.2' 

0,7) 

2,7) 

10.6' 

0.9) 

3,3) 

11.9' 

1.1) 

3,8) 

12,9' 

1.2) 

1.5) 

alC 

16.1' 

0.7) 

2.0) 

17,7' 

1.0) 

2,5) 

19,0' 

1.2) 

3,0) 

21, 1< 

1,5) 

3.5) 

22.4' 

1.6) 
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Table  El  (Concluded) 


table  or  extremes  estimates 

GRTD  LOCATION  3.  7 LAT=45.75  L*)N  = 83.30  COCKbURN  ISLAND 

Mi 


SHORELINE  GRID  POINT  28 


1 


5 

13. 8< 

1,6) 

10 

15. 4< 

2,1) 

20 

16. 7( 

2.7) 

50 

17. 7< 

3.3) 

lOQ 

18. 7( 

3.8) 

1 


5 

7 . 1 ( 

1.5) 

10 

7 , 9 ( 

2.0) 

20 

8 . 7 ( 

2.5) 

50 

9 . 1 ( 

3.0) 

lOQ 

9 . 4 ( 

3.5) 

winter 

ANGLE  CLASSES 
2 3 


17. 1C 

1.6) 

12  • 8 ( 

18. 0< 

2.1) 

14. 1( 

18. 7( 

2.7) 

15. 1( 

19. 4( 

3.3) 

16. 1< 

19. 7< 

3 . B ) 

16- 4( 

SPRING 

ANGLE  CLASSES 

2 

3 

9 . 1 ( 

1.5) 

7 . 5 ( 

11.  o< 

2.0) 

9 . 1 ( 

12. 6( 

2.5) 

10. 2* 

14. 2( 

3,0) 

11. 4( 

15. 4( 

3.5) 

12  • 2 * 

ALL 

1,6) 

17. 1* 

0.4) 

2,1) 

1 8 . 1 * 

0.6) 

2,7) 

18.8* 

0,7) 

3,3) 

19.41 

0.9) 

3,8) 

19, 7 < 

1,0) 

1,5) 

ALL 

9.2* 

0.7) 

2.0) 

H.l( 

0,9) 

2.5) 

12. 7* 

1.2) 

3.0) 

14,3* 

1.4) 

3.5) 

15,5* 

1.7) 

1 


5 

7 . 9 ( 

1.6) 

10 

9 . 5 ( 

2.1) 

20 

10. 8( 

2.7) 

50 

11. 8( 

3.3) 

lOQ 

12. 5( 

3.8) 

1 


5 

i4 . 8 ( 

1.5) 

10 

16. 1C 

2.0) 

20 

17. 1( 

2.5) 

50 

18. 0< 

3,0) 

100 

18. 7( 

3.5) 

SUMMER 

*ngle  classes 

2 3 


8 • 5 ( 

1.6) 

2-0* 

10.2* 

2.1) 

3 • 3 1 

11 . 2 ( 

2.7) 

4 . 3 ( 

12- 5( 

3,3) 

5 - 2 < 

13. 1( 

3 . b ) 

5 • 9 * 

FAlL 

ANGLE  CLASSES 

2 

3 

16- 4( 

1.5) 

1?  • 5 < 

17. 7( 

2.0) 

16- K 

18. 7< 

2.5) 

18- 7( 

19. 7( 

3.0) 

2l-3( 

20- 3( 

3.5) 

23. 0* 

1.6) 

ALL 

8.9* 

0.7) 

2.1) 

10.2* 

0,9) 

2,7) 

11.2* 

1,1) 

3,3) 

12.6* 

1.4) 

3,8) 

13.2* 

1.6) 

all 

1.5) 

16,7* 

0.5) 

2.D 

18.0* 

0,7) 

2,6) 

18.8* 

0,9) 

3,2) 

2l.4< 

1,1) 

3,7) 

23.1* 

1.2) 
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Table  E2 


wave  heights  tor  approach  directions  and  seasons  combined 


site  return  per i ops 


5 

10 

20 

50 

100 

1 

20.2 

20,7 

21.1 

21.7 

22,1 

2 

19.1 

19.7 

20.3 

21.5 

22.5 

3 

15.5 

16,4 

17.3 

18.5 

19,4 

4 

16.3 

17.3 

16.4 

19.9 

21.0 

5 

12.0 

12,0 

12,0 

21.2 

22,3 

6 

17,6 

18.7 

19.8 

21.2 

22,3 

7 

10.8 

19.7 

20,8 

22.3 

23,4 

6 

19.1 

20.1 

21.1 

22.4 

23,5 

9 

19.2 

20.1 

21.0 

22.2 

23,2 

10 

19.7 

20.8 

21.9 

23.4 

24,5 

11 

20.4 

21.1 

21.8 

23.5 

25,1 

12 

19.5 

20.4 

21.4 

23.0 

24,5 

13 

19.2 

20.1 

21.3 

23.1 

24,8 

l4 

20.2 

21.0 

21,8 

23.0 

24,0 

15 

20.1 

21. 0 

21.9 

23.3 

24.5 

16 

20.1 

21.1 

22.1 

23.7 

25,0 

1? 

20.0 

21.1 

22,3 

24,0 

25,5 

1® 

19.6 

20.7 

21.8 

23.5 

25.0 

1« 

17.1 

18.1 

19.1 

20.6 

21.8 

20 

15.2 

16.1 

17.1 

18.5 

19.5 

21 

16.3 

17.2 

18.0 

19.2 

20,0 

22 

15.3 

16.1 

16.9 

17.9 

18,8 

2® 

13.9 

14.5 

15.1 

15.9 

16,5 

24 

17,6 

18.6 

19.7 

21,2 

22,3 

25 

16.5 

17.7 

18.9 

20.5 

21,7 

26 

16.2 

17.5 

18.8 

20.5 

21,8 

27 

17.7 

18.5 

19.5 

21,1 

22,5 

28 

18.3 

19,0 

19.7 

21,4 

23,1 

E31 


1 


Table  E3 


GRID  LOCATION  22,12  lAT=43.02  lCN»E?.33 

GRID  POINT  NUMBER  1 


PORT  HURON 
MI 


significant  period  by  angle  class  ami  wavf  height 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

1 

1*6 

1.6 

2.0 

2 

3.2 

2.8 

3.2 

3 

4.4 

3.9 

4.3 

4 

5.3 

4.7 

5.2 

5 

6.0 

5.4 

5.7 

6 

6.2 

5.7 

5.9 

7 

6,4 

6.1 

6.2 

8 

6.7 

6.4 

6.4 

9 

6.9 

6.7 

6.7 

ID 

7.1 

7.1 

6.9 

11 

7.3 

7.4 

7.1 

12 

7.5 

7.7 

7.4 

13 

7.8 

8.0 

7.6 

14 

8.0 

8.4 

7.9 

15 

8.2 

8.7 

8.1 

16 

8.4 

9.0 

8.3 

17 

8.6 

9.4 

8.6 

18 

8.9 

9.7 

8.8 

19 

9.1 

10.0 

9.1 

20 

9.3 

10,4 

9.3 

21 

9,5 

10.7 

9.5 

22 

9.7 

11.0 

9.8 

23 

10.0 

11.3 

10.0 

24 

10.2 

11.7 

10.3 

25 

10.4 

(Continued) 

12.0 

(Sheet  1 

10.5 

of  28) 
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Table  E3  (Continued) 


GRID  LOCATION  21,12  LAT»43.16  L0bia8?.33 

GRID  POINT  NUMBER  2 


lakeport 

MI 


significant  period  by  angle  class  ak n wave  height 

WAVE  HEIGHT  (FT)  ANGLE  CLASS 


1 

2 

3 

1 

3,0 

2.6 

2.3 

2 

4.4 

4.2 

3.6 

3 

5.2 

5.0 

4.5 

4 

5.6 

5.6 

5.4 

5 

5.8 

6.0 

6.0 

6 

5.9 

6.3 

6.3 

7 

6.0 

6.5 

6.7 

8 

6.1 

6.8 

7.0 

9 

6.2 

7.0 

7.3 

10 

6,3 

7.3 

7.6 

11 

6.4 

7.5 

8.0 

12 

6.5 

7.6 

8.3 

13 

6.6 

8.0 

8.6 

14 

6.7 

8.3 

9.0 

15 

6.8 

8.5 

9.3 

16 

6.9 

6.8 

9.6 

17 

7.0 

9.0 

10.0 

18 

7.1 

9.3 

10.3 

19 

7.2 

9.5 

10.6 

20 

7,3 

9.8 

10.9 

21 

7.4 

10.0 

11.3 

2? 

7.5 

10.3 

11,6 

23 

7.6 

10.5 

11.9 

24 

7.7 

10.8 

12.3 

25 

7.8 

11.0 

12.6 

(Continued) 
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Table  E3  (Continued) 


GRID  LOCATION  20.12  LAT=43,3l  L0N«8?.33  LEXINGTON 

Ml 

GRID  POINT  NUMBER  3 


significant  period  by  angle  class  ak>  wave  height 


WAVE  HEIGHT  (FT) 

ANGLE  CLASS 

1 

2 

3 

1 

3.4 

3.3 

1.4 

2 

4.8 

4 . 4 

3.2 

3 

5.6 

5.0 

4.6 

4 

6.1 

5.5 

5.8 

5 

6.5 

6 . 0 

6.6 

6 

6.7 

6.3 

6.9 

7 

6.9 

6.5 

7.2 

8 

7.2 

6.8 

7.6 

9 

7.4 

7.0 

7.9 

m 

7.6 

7.3 

8.2 

n 

7.8 

7.6 

8.5 

12 

8.0 

7.8 

8.8 

13 

8.3 

8.1 

9.2 

14 

8.5 

8.3 

9.5 

15 

8.7 

8.6 

9.8 

16 

8.9 

8.9 

10. 1 

17 

9.1 

9.1 

10.4 

18 

9.4 

9.4 

10.8 

19 

9.6 

9.6 

11.1 

20 

9.8 

9.9 

11.4 

21 

10.0 

10.2 

11.7 

22 

10.2 

10.4 

12.0 

23 

10.5 

10 . 7 

12.4 

24 

10.7 

10.9 

12.7 

25 

10.9 

11.2 

13.0 

(Continued) 
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Table  E3  (Continued) 


GRID  LOCATION  19,11  l ATa43 . 45  L0h)*e?.52 
GRID  POINT  NUMBER  4 


PORT  SANILAC 
MI 


SIGNIFICANT  PFRIOD  BY  ANGLE  CLASS  Ain  WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

1 

3.6 

2.8 

1.8 

2 

5.0 

4.2 

3.4 

3 

5.6 

5.0 

4.5 

4 

6.0 

5.5 

5.5 

5 

6.2 

5.9 

6.2 

6 

6.3 

6.1 

6.5 

7 

6.5 

6.4 

6.8 

ft 

6.6 

6.6 

7.0 

9 

6.8 

6.9 

7.3 

in 

6.9 

7.1 

7.6 

11 

7.0 

7.3 

7.9 

12 

7.2 

7.6 

8.2 

13 

7.3 

7.8 

8.4 

14 

7.5 

8.1 

8.7 

15 

7.6 

8.3 

9.0 

16 

7.7 

8.5 

9.3 

17 

7.9 

8.8 

9.6 

18 

8.0 

9.0 

9.8 

19 

8.2 

9.3 

10. 1 

2H 

8.3 

9.5 

10.4 

21 

8.4 

9.7 

10.7 

22 

8.6 

10.0 

11.0 

23 

8.7 

10.2 

11.2 

24 

8.9 

10.5 

11.5 

25 

9.0 

10.7 

11.  f 

(Continued) 
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Table  E3  (Continued) 


r 


GRID  LOCATION  18,11  LAT=43.59  L0&«82.52 

GRID  POINT  NUMBER  5 


forestville 

MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AK  D WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

l 

3,2 

3.0 

1.6 

2 

4.7 

4.4 

3.4 

3 

5.6 

5.3 

4.6 

4 

6.0 

5.7 

5.6 

5 

6.4 

6.1 

6.4 

6 

6.6 

6.3 

6.7 

7 

6,8 

6.5 

6.9 

8 

7.0 

6.8 

7.2 

9 

7.2 

7.0 

7.4 

10 

7.4 

7.2 

7,7 

11 

7.6 

7.4 

8.0 

12 

7.8 

7.6 

8.2 

13 

8.0 

7.9 

8.5 

14 

8.2 

6.1 

8.7 

15 

8.4 

8.3 

9.0 

16 

8.6 

8.5 

9.3 

17 

8.8 

8.7 

9.5 

18 

♦ .0 

9.0 

9.8 

19 

9.2 

9.2 

10.0 

20 

9.4 

9.4 

10.3 

21 

9.6 

9.6 

10.6 

22 

9.8 

9.8 

10.8 

23 

10,0 

10.1 

11.1 

24 

10.2 

10.3 

11.3 

25 

10.4 

10.5 

11.6 

(Continued) 
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Table  E3  (Continued) 


GRID  LOCATION  17,11  lAT=43.74  LCNs6?;52 

GRID  POINT  NUMBER  6 


HELENA 

MI 


significant  period  by  angle  class  ami  wave  height 


WAVE  HEIGHT  (FT) 


angle  class 


1 

2 

3 

1 

3,2 

2.4 

2.6 

2 

4,7 

4.0 

4.2 

3 

5,6 

4.8 

5.0 
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Table  E3  (Continued) 


GRID  LOCATION  16,11  LAT»43.08  LCS«fc?.52 

GRID  POTNT  NUMBER  7 

SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  Ah n WAVE 


HARBOR  BEACH 
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Table  E3  (Continued) 


GRID  LOCATION  15,10  LATa44,03  L0hjse?.71 

GRID  POINT  NUMBER  8 


HURON  CITY 
MI 


SIGNIFICANT  PpRlOD  BY  ANGLE  CLASS  AND  WAVE  HEIGHT 
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Table  E3  (Continued) 


GRID  LOCATION  15.  8 lAT*44.03  L0Klsfr3.il 

GRID  POINT  NUMBER  9 


PORT  CRESCENT 
MI 


significant  period  BY  ANGLE  CLASS  AND  nave  height 
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Table  E3  (Continued) 


GRID  LOCATION  15,  7 
GRID  POINT  NUMBER  10 


LATs44,03  LOblsE3,30  ENTRANCE  SAGINAW  BAY 

MI 


SIGNIFICANT  PFRIOD  BY  ANGL c CLASS  AM)  WAVE  HEIGHT 
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Table  E3  (Continued) 


GRID  LOCATION  14,  7 LATs44.17  Leb|s63.30 

GRID  POINT  NUMBER  11 


TAWAS  CITY 
HI 


significant  period  by  angle  class  ami  wave  height 
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Table  E3  (Continued) 


GRID  LOCATION  13,  8 LAT=44,3l  L0RbE3.11 

grid  point  NUMBER  12 


OSCODA  AU  SABLE 
MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AMI  WAVE  HEIGHT 


WAVE  HEIGHT  (FT)  ANGLE  CLASS 
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Table  E3  (Continued) 


GRID  LOCATION  12,  8 LAT*44.46  L0Saf3.11 

GRID  POtNT  NUMBER  13 


GREENBUSH 

MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AMI  WAVE  HEIGHT 
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Table  E3  (Continued) 


GRID  LOCATION  11,  8 LATs44.60  L0Hsti3.ll 

GRID  POTNT  NUMBER  14 


HARRJSVTLLE 

HI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AID  RAVE  HEIGHT 


WAVE  HEIGHT  (FT) 
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Table  E3  (Continued) 


GRID  LOCATION  10,  8 |_AT  = 44.74  L0Ns63.ll 
GRID  POTNT  NUMBER  15 


black  river 

MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  Af*n  WAVE  HEIGHT 
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Table  E3  (Continued) 


GRID  LOCATION  9.  8 LATa44.89  L0hls83.ll 

GRID  POINT  NUMBER  16 


OSS  1 Nf-KE 
MI 


SIGNIFICANT  PpRIOD  BY  ANGLE  CLASS  Ah D WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 
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Table  E3  (Continued) 


GRID  LOCATION  8,  8 lATs45.04  LObiafJ.il 

GRID  POINT  NUMBER  17 


NORTH  POINT 
MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AKn  WAVE  HEIGHT 


WAVE  HEIGHT  (FT)  ANGLE  CLASS 
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Table  E3  (Continued) 


GRID  LOCATION  7,  7 LATs45.l8  LOSsES. 30 
GRID  POINT  NUMBER  18 


ROCKPORT 

MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  A^ n WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 
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GRID  LOCATION  6,  7 LATa45.32  LOS|s£3;30 

GRID  POINT  NUMBER  19 


STONE  PORT 
MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AKD  WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 
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Table  E3  (Continued) 


GRID  LOCATION  5,  6 L*T*45.46  LOfi«E3.50 

GRID  POINT  NUMBER  20 


ADAMS  POINT 
MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AMI  WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

1 

3.6 

3.4 

2.8 

9 

4.8 

4.4 

4.2 

3 

5.3 

5.0 

5.0 

4 
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5.5 

5.4 

5 

6.1 
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5.8 

6 

6.4 
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6.1 

7 
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6.4 

6.4 
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6.6 

9 
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8.6 
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16 

8.9 

8.9 

8.9 
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9 . 4 
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Table  E3  (Continued) 


GRID  LOCATION  5,  5 LAT=45.46  LOK|af3)70 

GRID  POINT  NUMBER  21 


ROGERS  CITY 
MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AM}  WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 
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6.9 

10 

7.3 

7.1 

7.2 

11 
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Table  E3  (Continued) 


GRID  LOCATION  4.  4 LAT345.6l  L0b*83i90 

GRID  POINT  NUMBER  22 


HAMMOND  BAY 
MI 


significant  period  by  angle  class  ami  wave  height 


WAVE  HEIGHT  (FT) 

angle  class 
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4.7 
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16 

9.3 
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8.9 
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8.9 
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9.1 
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9.6 
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21 
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22 
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10.6 

10.9 

25 

11.8 

10.8 

11.2 

(Continued) 


(Sheet  22  of  28) 


Table  E3  (Continued) 


GRID  LOCATION  3,  3 

LAT  = 45.75  L0&|*P4 

709  C0RDW00D 

POINT 

MI 

GRID  POINT 

NUMBER  23 

SIGNIFICANT 

period  by 

ANGLE  CLASS  AM) 

WAVE  HEIGHT 

WAVE  HEIGHT 

(FT) 

angle  class 

1 

2 

3 

1 

3.8 

2.6 

3.6 

2 

4.7 

4.0 

4.4 

3 

5,2 
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5.0 

4 
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5 
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6.0 

5.9 

6 

6.5 
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6.2 

7 

6.8 
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6.4 

8 

7.0 
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6.7 

9 

7.3 

7.1 

7.0 

10 

7.6 

7.4 

7.3 

11 

7.9 

7.6 
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12 
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7.9 

7.8 

13 

8.4 

8.2 

8.1 

14 

8.7 

8.4 

8.3 

15 

9.0 

8.7 

8.6 

16 

9.3 

9.0 

8.9 

17 

9.6 

9.2 

9.1 

19 

9.8 

9.5 

9.4 

19 

10.1 

9.8 

9.7 

20 

10.4 
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10.0 

21 

10.7 

10.3 

10.2 

22 

11.0 
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10.5 

23 
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10.8 

24 
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11.0 

25 
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11.4 
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Table  E3  (Continued) 


GRID  LOCATION  2,  3 lAT*45.89  LC&|s84.09 

GRID  POINT  NUMBER  24 


POINT  DOLOMITE 
Ml 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AK n WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

1 

3.8 

3.8 

3.7 

2 

4.8 

4.6 

4.8 

3 

5,3 
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5.3 

4 

5.7 

5,6 

5.8 

5 

6.1 

6.1 

6.2 

6 

6.4 

6.4 

6.5 

7 

6.6 

6.7 

6.7 

8 

6.9 

7.0 

7.0 

9 

7.1 

7.3 

7.2 

10 

7.4 

7.6 

7.5 

11 

7.6 

7.8 

7.8 

12 

7.9 

8.1 

8.0 

15 

8.1 

8.4 

8.3 

14 

8.4 

8.7 

8.5 

15 

8.6 

9.0 

8.8 

16 

8.9 

9.3 

9.1 

17 

9.1 

9.6 

9.3 

15 

9.4 

9.9 

9.6 

19 

9.6 

10.2 

9.8 

20 

9.9 

10.5 

10.1 

21 

10.1 

10.7 

10.4 

22 

10.4 

11.0 

10.6 

23 
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11.3 
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11.1 
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11.1 

(Continued) 
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Table  E3  (Continued) 


GRID  LOCATION  2.  4 LAT*45.89  L05j  = t3.;90  DETOUR  REEF 

MI 

GRID  POINT  NUMBER  25 

significant  period  by  angle  class  ami  wave  height 

WAVE  HEIGHT  (FT)  ANGLE  CLASS 


1 

2 

3 

1 

3.9 

3.4 

3.6 

2 

4.8 

4.6 

4.8 

3 

5.3 

5.3 

5.3 

4 

5.6 

5.7 

5.8 

5 

5.9 

6.2 

6.2 

A 

6.1 

6.5 

6.5 

7 

6.3 

6.8 

6.8 

3 

6.5 

7.0 

7.0 

9 

6.7 

7.3 

7.3 

11 

6.9 

7.6 

7.6 

11 

7.0 

7.9 

7.9 

12 

7.2 

8.2 

8.2 

13 

7.4 

8.4 

8.4 

14 

7.6 

8.7 

8.7 

15 

7.8 

9.0 

9.0 

16 

8.0 

9.3 

9.3 

17 

8.2 

9.6 

9.6 

19 

8.4 

9.8 

9.8 

19 

8.6 

10.1 

10.1 

21 

8.8 

10.4 

10.4 

21 

8.9 

10.7 

10.7 

22 

9.1 

11 . 0 

11.0 

23 

9.3 

11.2 

11.2 

24 

9.5 

11.5 

11.5 

25 

9.7 

11.8 

11.8 
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Table  E3  (Continued) 


GRID  LOCATION  3.  5 
GRID  POINT  NUMBER  26 


L ATs  45 , 75  L0fcl»f3.70  WEST  END  DRlMMOND  ISLAND 

MI 


significant  period  by  angle  class  ami  wave  height 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

1 

3.8 

3.6 

3.6 
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4.8 

4.8 

4.6 
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6 
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3 

6,6 
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7.0 

9 
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7.2 

in 
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7.4 
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11 

7.2 

7.6 

7.8 

12 

7.4 

7.8 

8.0 

13 

7.6 

8.0 

8.3 

14 

7.8 

8.2 

8.5 

15 

8.0 

8.4 

8.8 

16 

8.2 

8.6 

9.1 

17 

8.4 

8.8 

9.3 

18 

8.6 

9.0 

9.6 

19 

8.8 

9.2 

9.8 

20 

9.0 
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10.1 

21 

9.2 

9.7 

10.4 

22 

9.4 

9.9 
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23 

9.6 

10.1 

10.9 
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25 

10.0 

(Continued) 
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Table  E3  (Continued) 


GRID  LOCATION  3.  6 
GRID  POINT  NUMBER  27 


LAT345.75  L0&j  = f3.50  FALSE  DETOUR  CHANNEL 

MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AND  WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 
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9 
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7.0 

10 
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7.3 

7.3 

11 
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7.6 

7.6 

12 

7.9 

7.8 

7.9 

13 

8.1 

8.1 

8.1 

14 

8.4 

8.3 

8.4 

15 

8.6 

8.6 

8.7 

16 

8.8 

8.9 

9.0 

17 

9.1 

9.1 

9.3 
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9.9 

10.1 

21 

10.0 
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Table  E3  (Concluded) 


GRID  LOCATION  3,  7 |_AT  = 45.75  LCK|sfc3.30 
GRID  POINT  NUMBER  26 


COCKBURN  ISLAND 
MI 


SIGNIFICANT  PERIOD  BY  ANGLE  CLASS  AN  n WAVE  HEIGHT 


WAVE  HEIGHT  (FT) 

angle  class 

1 

2 

3 

1 

3.6 

3.8 
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4.7 
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5.5 

5 

6.0 

6 . 0 
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A 

6.3 

6.3 

6.2 

7 

6.5 

6.6 

6.4 

8 

6.8 

6.9 

6.7 

9 

7.0 

7.2 

6.9 

n 

7.3 

7.4 

7.2 

ii 

7.6 

7.7 

7.5 

12 

7.8 

8.0 

7.7 

ii 

8.1 

8.3 

8.0 

14 

8.3 

8.6 

8.2 

15 

8.6 

8.9 

8.5 

18 

8.9 
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8.8 

17 

9.1 
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9.0 
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9.4 
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9.6 
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10.9 

10.3 

21 

10.7 

11.2 

10.6 

24 

10.9 

11.5 
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Table  E4 

Azimuths  of  Normal  Shoreline  Vectors* 


Shoreline 

Point 

Azimuth 

Shoreline 

Point 

Azimuth 

1 

180 

15 

270 

2 

2U0 

16 

270 

3 

270 

17 

270 

4 

270 

18 

240 

5 

270 

19 

210 

6 

270 

20 

210 

7 

240 

21 

210 

8 

210 

22 

210 

9 

180 

23 

210 

10 

210 

24 

330 

11 

300 

25 

0 

12 

300 

26 

0 

13 

300 

27 

0 

ll; 

270 

28 

0 

* See  Figure  21 


APPENDIX  F:  NOTATION 


a 


Vas 


n 

b 


n 

B 


VBs 


G 

d 

d. 

J 

e 

E 

E(f) 

Ejr) 

f 


max 

frf2 

F 


F( f ,0 ) 
F(x) 
F(x,k,t ) 
F(H) 
F(L) 
F(io) 

g 

G 


Regression  constant  (Equation  2),  and  Phillips  equilibrium 
constant  (Equation  B9) 

Regression  coefficients  relating  wind  speeds  at  different 
elevations  ( S.  = above  20  ft,  s = at  20  ft) 

Parameter  in  asymptotic  distribution  of  extreme  values 

Regression  constant 

Parameter  in  asymptotic  distribution  of  extreme  values 

Linear  constant  for  log-linear  wind  profile  (Appendix  A), 
and  coefficient  of  energy  transfer  and  function  of  propaga- 
tion direction  relative  to  wind  direction  and  wind  speed 
(Appendix  B) 

Regression  coefficients  relating  wind  speeds  at  different 
elevations  ( l = above  20  ft,  s = at  20  ft) 

Coefficient  of  bottom  friction 

Group  velocity 

Water  depth 

Distance  between  wind  station  and  interpolation  point 
Napierian  base  = 2.71828... 

Total  wave  energy 

One-dimensional  representation  of  wave  spectrum 

Fully  developed  wave  spectrum 

Frequency 

Frequency  associated  with  maximum  spectral  density 

Lowest  and  highest  frequency,  respectively,  in  the  model 

Land-bounded  fetch  (Equation  10),  and  wave  component  (Equa- 
tion B15) 

Two-dimensional  representation  of  wave  spectrum 
Cumulative  probability  function  of  x 
Five-dimensional  representation  of  wave  spectrum 
Cumulative  probability  function  of  H 
Ice-limited  fetch 

One-dimensional  spectral  energy  density  and  function  of  w 

Function  (Equation  10)  and  gravitational  acceleration 

Geostrophic  wind  velocity  (Parts  I-IV)  and  function 
(Equation  A1 ) 


FI 


r l 


h(0)  Spreading  factor 

H Wave  height 

H'  Wave  height  transformed  by  fetch  reduction  (Equation  11 ) , 
and  wave  height  as  hindcast  by  model  (Equation  D9) 

H*  Reduced  wave  height  variate 

H Ice-limited  wave  height 

Li 

Significant  wave  height 

i Conventional  summation  over  all  dimensions 

k Parameter  denoting  curvature  in  extremal  distribution 

(Equation  5),  and  KEYPS  function  constant  (Equation  A13) 

— ► 

k Wave  number  vector 

K Von  Karman  constant 

l.  Series  product  of  distances  between  wind  stations  and  inter- 
polation point 

L'  Stability  length 

n Number  of  independent  observations  within  a season  or  year 

N Number  of  years  in  a sample  (Figure  23),  and  number  of  wind 
stations  (Appendix  C) 

NORM  Sum  of  multiple  products  Jl 

p(H)  Probability  of  a particular  wave  height  occurring  at  a site 

p'(H)  Probability  of  a particular  wave  height  occurring  estimated 

from  hindcasts  not  considering  ice  cover 

p(e)  Probability  of  no  ice  at  the  site 

q Ratio  of  energy  density  in  spectrum  to  that  in  fully 
developed  spectrum 

R Ratio  of  lake-wind  speed  to  land-wind  speed  (Equations  1 
and  2),  and  angular  spreading  function  in  Barnett  model 
(Equation  Bll) 

S Sum  of  all  exterior  sources  and  sinks  of  energy  transfer 
into  spectrum  (Equation  B7),  and  standard  deviation  of 
annual  maxima  (Equation  D8) 

SD  Sum  of  sources  in  Barnett  model 

S^  Slim  of  sources  in  Inoue  model 

S^  Sum  of  linear  source  terms 

S Standard  deviation  of  the  estimate 

x 

S „ Sum  of  nonlinear  source  terms 

nil 

S^-S^  Sources  involving  weak  interactions  between  winds  and  waves 


F2 


si-sh 


T' 

T 


w 

TrT2 


u 


# 


U 


u 

w 

V 


V 

0 

W. 

1 

x 


X 


X 

o 

X 


y 

z 


z 

o 


Sources  involving  weak  interactions  between  currents  and 
waves 

Source  involving  weak  interactions  among  different  spectral 
components  of  waves 

Sou  ce  involving  weak  interactions  between  bottom  and  waves 

Unknown  dissipation  function  due  to  breaking 

Time  (Equation  Bl),  and  wind  duration  (Equation  A3) 

Return  period  or  recurrence  interval  (Equation  6),  and  non- 
dimensional  time  factor  (Equation  Al) 

Return  period  as  modified  by  ice  cover 

Air  temperature 

Recurrence  interval  at  elevation  above  20  ft 

Lake  surface  temperature  (Figures  9 and  11 ),  and  recurrence 
interval  at  20-ft  elevation  (Equation  A 6) 

Water  temperature 

Transfer  functions  dependent  on  coupling  coefficients  be- 
tween spectral  components  involved  in  the  interactions 

Wind  speed  (Equation  AU ) , fluid  velocity  relative  to  the 
bottom  (Equation  B26),  and  parameter  of  the  double- 
exponential distribution  (Equation  D1 ) 

Friction  velocity 

East-west  wind  component 

Wind  speed  over  land 

Wind  speed  over  water 

Anemometer-level  wind  (Figure  11 ) , and  north-south  wind  com- 
ponent (Appendix  C) 

Wind  velocity  (Equation  A9) 

Weighting  of  interpolation  function  on  ith  wind  station 
Actual  fetch  (Equation  A2) 

Location  vector 

Parameter  of  extremal  distribution 

Nondimensional  fetcn  factor 

Distance  inland  from  lake-land  interface 

Reduced  variate  in  distribution  of  extremes 

Elevation  (Equation  Ak ) , and  parameter  of  the  combined 
double-exponential  and  Gaussian  distribution  (Equation  D15) 

Roughness  height 


F3 


_ 


r 


a 


8 


Y 

r 

6 


At 


E 


01,02 


n 


JJ 

o 


2 

a 


4> 


*(f,e) 

4> 

n 


Parameter  of  extremal  distribution  (Equation  5),  energy 
transfer  due  to  interaction  between  turbulent  pressure  fluc- 
tuations and  waves — Phillips'  mechanism  (Equation  BT ) , and 
parameter  of  double-exponential  distribution  (Equation  Dl) 

Energy  transfer  due  to  interaction  between  mean  wind  profile 
and  waves — Miles'  mechanism  (Equation  B9),  and  inverse  of 
a (Equation  D20) 

~-V  -f 

Function  of  k,kf  in  Equation  B7 
Parameter  of  wave-wave  interactions 
Integration  variable 
Land  boundary  layer  thickness 
Kronecker  delta 

Increment  of  time  in  numerical  equation 

Noiidimensional  wind  shear  (Equation  A12),  and  function  of 
k,k'  (Equation  B7 ) 

Direction  of  wave  propagation 

Direction  of  wind  heading 

Limiting  angles  for  wave  propagation  toward  the  shore 
Dynamic  viscosity  of  air 
Viscosity  tensor 
3.lUl6 

Nondimensional  wave  height 

Circular  frequency  associated  with  F (Equation  BU),  and 
standard  deviation  of  the  errors  in  hindcast  wave  heights 
(Equation  Dll) 

Variance  of  surface  in  a wave  spectrum  (related  to  total 
energy ) 

Parameter  of  wave-wave  interactions 

Modification  function  of  Z/L'  for  wind  profile  due  to 
thermal  effects  (Equation  A10),  and  limiting  function  for 
wind  input  in  Barnett  model  (Equation  BIO) 

Effect  of  ice  cover  at  site  on  return  periods  (Equation  8), 
and  probability  integral  (Equation  Dl6) 

Frictional  dissipation  function 

Normalized  function  dependent  on  air-sea  temperature 
difference  (Equation  2) 

Function  of  Z , which  relates  estimated  wave  height  dis- 
tribution to  actual  wave  height  distribution 


FU 


♦ Velocity-dependent  function  (Equation  2) 

0)  Angular  frequency 
wo  Wind  frequency 


In  accordance  with  letter  from  DAEN-RDC , DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Resio,  Donald  T 

Design  wave  information  for  the  Great  Lakes:  Report  4: 

Lake  Huron  / by  Donald  T.  Resio,  Charles  L.  Vincent.  Vicks- 
burg, Miss.  : U.  S.  Waterways  Experiment  Station,  1977. 

49,  c104j  p.  : ill.  ; 27  cm.  (Technical  report  - U.  S. 

Army  Engineer  Waterways  Experiment  Station  ; H-76-1,  Report  4) 
Prepared  for  U.  S.  Army  Engineer  Division,  North  Central, 
Chicago,  Illinois. 

References:  p.  45-49. 
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